Ngày tải lên :
24/12/2013, 12:16
... as
in the original Bulirsch-Stoer method.
The starting point is an implicit form of the midpoint rule:
y
n+1
− y
n−1
=2hf
y
n+1
+ y
n−1
2
(16.6.29)
738
Chapter 16. Integration of Ordinary Differential ... calculates
dydx
.
{
void lubksb(float **a, int n, int *indx, float b[]);
void ludcmp(float **a, int n, int *indx, float *d);
int i,j,nn,*indx;
float d,h,x,**a,*del,*ytemp;
indx=ivector(1,n);
a=matrix(1,n,1,n);
del=vector(1,n);
ytemp=vector(1,n);
h=htot/nstep; ... components of
y
.
{
void jacobn(float x, float y[], float dfdx[], float **dfdy, int n);
void lubksb(float **a, int n, int *indx, float b[]);
void ludcmp(float **a, int n, int *indx, float *d);
int i,j,jtry,*indx;
float...