微分方程数值解法课程设计报告--许.docVIP

- .-- 微分方程数值解法课程设计----A组 计算022班3号 许** 用复化梯形计算积分 解： 算法思想： 根据复化梯形公式:取步长， 有，写出程序，如下所示： N=100; a=0;b=10; h=(b-a)/N; T=zeros(0,N); for j=1:N T(j)=0.5*h*[exp(-(a+h*(j-1)).^2)+exp(-(a+h*j).^2)] t=sum(T) end 运行结果：t = 0.8862 用Euler、改进的Euler法、梯形法、R-K法解 解： 根据算法思想： （1）Euler法 （2）改进的Euler法 （3）梯形法 （4）R-K法 我们可以写出程序如下所示： format long; N=100; h=1/N;%步长 U=zeros(1,N+1);%euler T=zeros(1,N+1);%真解 X=zeros(1,N+1); R=zeros(4,N+1);%误差 G=zeros(1,N+1);%改进eule O=zeros(1,N+1);%梯形法 K=zeros(1,N+1)%R-K法 c=zeros(1,4); a=0; for j=2:N+1; X(j)=h*(j-1); end U(1)=1 T(1)=1; G(1)=1; X