微分方程-第一章习题5实验报告.docVIP

第一章习题5实验报告 一、实验目的 通过用Euler法、中点格式、预报－校正格式、Adams预报—修正格式等算法求解一阶常微分方程初值问题 的数值解，经典四级四阶R-K法是高精度单步法，四阶Adams预估－校正法是线性多步法，这两个算法较实验一算法精度要高。通过本次实验，掌握R-K法的编程要领，掌握多步法初值的准备方法，深刻领会微分方程数值解的实质，体会单步法和线性多步法各自的优缺点，熟练掌握各算法的计算机实现过程，并能从理论及实验结果分析各种算法的优缺点。 二、实验内容 用Euler法、中点格式、预报－校正格式、经典四级四阶R-K格式、Adams预报—修正格式算法求解一阶常微分方程初值问题 的数值解，其精确解为。 三、算法 1 Euler法 其中，为步长。 2　中点格式 3　预报－校正格式 4 四级四阶格式： 5 Adams预报—修正格式 其中为步长。 四 程序 一、Euler法程序如下： #include stdio.h #define LEN 100 int Euler(int a,int b,double h) { double y[LEN]={0}; double x[LEN]={0}; int i=0; double num=0; num=((b-a)/h)+1;//节点个数 //循环计算出各个节点的坐标 for(i=1;inum;++i) { x[0]=a; x[i]=a+i*h; } //带公式计算各个节点的y的值 for(i=0;inum;i++) { y[0]=1; y[i+1]=(y[i]/(2*x[i]))+(x[i]*x[i])/(2*y[i]); } //输出x和y printf(\n ************************* Euler法求解方程 ******************************* \n); printf(\n x[n] y[n]\n\n); for(i=0;inum;i++) { printf( %.4f %.4f\n,x[i],y[i]); } printf(\n ************************* Euler法求解方程 ******************************** \n\n); return 0; } } 二、中点格式程序如下： int Mid(int a,int b,double h) { double y[LEN]={0}; double x[LEN]={0}; int i=0; double num=0; num=((b-a)/h)+1;//节点个数 //循环计算出各个节点的坐标 for(i=0;inum;++i) { x[0]=a; x[i]=a+i*h; } //带公式计算各个节点的y的值 for(i=0;inum;i++) { y[0]=1; y[i+1]=y[i]+h*f((x[i]+h/2),(y[i]+(h/2)*(f(x[i],y[i])))); } //输出x和y printf(\n ************************* 中点格式求解方程 ******************************* \n); printf( x[n] y[n]\n\n); for(i=0;inum;i++) { printf( %.4f %.4f\n,x[i],y[i]); } printf(\n ************************* 中点格式求解方程 ******************************* \n\n); return 0; } 三、预报——校正格式程序如下： int PreCorr(int a,int b,double h) { double y1[LEN]={0}; double y[LEN]={0}; double x[LEN]={0}; int i=0; double num=0; num=((b-a)/h)+1;//节点个数 //循环计算出各个节点的坐标 for(i=0;inum;++i) { x[0]=a