实验五 使用Euler方法和改进的Euler方法求解初值问题

实验五 使用Euler方法和改进的Euler方法求解初值问题
实验目的:用Euler方法和改进的Euler方法求解初值问题。
实验内容:用改进欧拉方法解初值问题

取步长h=0.1计算,并与准确值 y = -x-1+2ex相比较。第一行输出0.1这个节点的准确值、Euler值、Euler值的误差、改进的Euler值、改进的Euler值的误差。分析一下Euler方法和改进的Euler方法哪个误差比较小,说明一下原因。
实验源程序:
下面这一段是参考程序,使用Euler方法求解初值问题,同学们可以分开编写Euler方法和改进的Euler方法。

#include 
#include 
#define f(x,y) (x+y)
#define y(x) (-x-1+2*exp(x))
void main()
{
	float a=0,b=1.0,h=0.1,y0=1.0,ye,yx,x;
	printf("\n分点   准确解     Euler解     Euler误差\n");
	x=a;
	ye=y0;//ye表示使用Euler法求解的数值解
	yx=y(0);//yx表示每个节点的准确值
	printf("%3.1f   %8.6f    %8.6f    %8.6f\n", x,yx,ye,fabs(ye-yx));
	while(x<b)
	{
		ye=(第一个空);
		x=(第二个空);
		yx=(第三个空);
        printf("%3.1f   %8.6f    %8.6f    %8.6f\n", x,yx,ye,fabs(ye-yx));
	}
}

你可能感兴趣的:(c语言)