分别用 欧拉法 和 四阶龙格-库塔法 解微分方程dy/dx=xy ,x=[0,3] ,x=0 时 y=1请写出matlab源程序,和运行结果,

来源：学生作业帮助网 编辑：六六作业网 时间：2024/11/18 01:47:58
分别用欧拉法和四阶龙格-库塔法解微分方程dy/dx=xy,x=[0,3],x=0时y=1请写出matlab源程序,和运行结果,分别用欧拉法和四阶龙格-库塔法解微分方程dy/dx=xy,x=[0,3],

分别用 欧拉法 和 四阶龙格-库塔法 解微分方程dy/dx=xy ,x=[0,3] ,x=0 时 y=1请写出matlab源程序,和运行结果,
分别用 欧拉法 和 四阶龙格-库塔法 解微分方程
dy/dx=xy ,x=[0,3] ,x=0 时 y=1
请写出matlab源程序,和运行结果,

分别用 欧拉法 和 四阶龙格-库塔法 解微分方程dy/dx=xy ,x=[0,3] ,x=0 时 y=1请写出matlab源程序,和运行结果,
f=inline('x*y','x','y'); %微分方程的右边项
dx=0.05; %x方向步长
xleft=0; %区域的左边界
xright=3; %区域的右边界
xx=xleft:dx:xright; %一系列离散的点
n=length(xx); %点的个数
y0=1;
%%(1)欧拉法
Euler=y0;
for i=2:n
Euler(i)=Euler(i-1)+dx*f(xx(i-1),Euler(i-1));
end
%%(2)龙格库塔法
RK=y0;
for i=2:n
k1=f(xx(i-1),RK(i-1));
k2=f(xx(i-1)+dx/2,RK(i-1)+k1*dx/2);
k3=f(xx(i-1)+dx/2,RK(i-1)+k2*dx/2);
k4=f(xx(i-1)+dx,RK(i-1)+k3*dx);
RK(i)=RK(i-1)+dx*(k1+2*k2+2*k3+k4)/6;
end
%%Euler和Rk法结果比较
plot(xx,Euler,xx,RK)
hold on
%精确解用作图
syms x
rightsolve=dsolve('Dy=x*y','y(0)=1','x');%求出解析解
rightdata=subs(rightsolve,xx);%将xx代入解析解,得到解析解对应的数值
plot(xx,rightdata,'r*')
legend('Euler','Runge-Kutta','analytic')

分别用 欧拉法 和 四阶龙格-库塔法 解微分方程dy/dx=xy ,x=[0,3] ,x=0 时 y=1请写出matlab源程序,和运行结果, 求这题的解集.分别用区间和数轴表示 韩国和朝鲜分别用英语怎么说? 什么时候和日期分别用英语表示 “初中”和“高中”分别用英语怎么说? “工资”和“薪水”用英文分别怎么说? “初中”和“高中”分别用英语怎么说? 字和词分别用英语怎么说 英语翻译分别用flight和plane or和and 分别怎么用?英语? 分别用一句话赞美父亲和母亲 price和prices分别在什么时候用? “单数”和“复数”分别用英语怎么说? 英语翻译分别用WHAT 和HOW 开头 分别用人迹罕至和人声鼎沸造句 请用“不足为奇”和“欣喜若狂”分别造句 上课和下课分别用英语怎么说? 用妒忌造句和捐赠分别造句