#include"stdio.h"
#define N 20
void ModEuler(float(*f)(float,float),float x0,float y0,float xn,int n)
{ int i;
float yp,yc ,x=x0,y=y0,h=(xn-x0)/n;
printf("x[0]=%f\ty[0]=%f\n",x,y);
for(i=1;i<=n;i++)
{
yp=y+h*(*f)(x,y);
x=x0+i*h;
yc=y+h*(*f)(x,yp);
y=(yp+yc)/2;
......................
阅读全部 | 2015年1月22日 15:45
#include"stdio.h"
#define N 20
void ModEuler(float(*f)(float,float),float x0,float y0,float xn,int n)
{ int i;
float yp,yc ,x=x0,y=y0,h=(xn-x0)/n;
printf("x[0]=%f\ty[0]=%f\n",x,y);
for(i=1;i<=n;i++)
{
yp=y+h*(*f)(x,y);
x=x0+i*h;
yc=y+h*(*f)(x,yp);
y=(yp+yc)/2;
......................
阅读全部 | 2015年1月22日 15:44
#include"stdio.h"
#define N 20
void ModEuler(float(*f)(float,float),float x0,float y0,float xn,int n)
{ int i;
float yp,yc ,x=x0,y=y0,h=(xn-x0)/n;
printf("x[0]=%f\ty[0]=%f\n",x,y);
for(i=1;i<=n;i++)
{
yp=y+h*(*f)(x,y);
x=x0+i*h;
yc=y+h*(*f)(x,yp);
y=(yp+yc)/2;
......................
阅读全部 | 2015年1月22日 15:44