改进欧拉法.txt

数值计算方法所用的上机实验C++源程序,个人所编.

#include <stdio.h>
float func(float x,float y)
{  return(y-x);
}

float euler(float x0,float xn,float y0,int N)
{  float x,y,yp,yc,h;
   int i;
   x=x0;
   y=y0;
   h=(xn-x0)/(float)N; /* 计算步长 */
   for(i=1;i<=N;i++) /* 改进的欧拉公式 */
   {  yp=y+h*func(x,y);
      x=x0+i*h;
      yc=y+h*func(x,yp);
      y=(yp+yc)/2.0;
   }
   return(y);
}
main()
{  float x0,xn,y0,e;
   int N;
   clrscr();
   printf("\ninput n:\n    "); /* 输入区间等分数 */
   scanf("%d",&N);
   printf("input x0,xn:\n    "); /* 输入x的区间 */
   scanf("%f%f",&x0,&xn);
   printf("input y0:\n    "); /* 输入x0处的y的值 */
   scanf("%f",&y0);
   e=euler(x0,xn,y0,N);
   printf("y(%f)=%6.4f",y0,e);
}