作业帮 > 综合 > 作业

MATLAB求方程y’= y - 2x/y,y(0) = 1 (0≤x≤1,h = 0.1) 的数值解,要求编写程序

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/07/17 23:02:03
MATLAB求方程y’= y - 2x/y,y(0) = 1 (0≤x≤1,h = 0.1) 的数值解,要求编写程序
用向前欧拉公式和改进的欧拉公式、龙格那什么,三种方法(要是嫌麻烦的话,只给出一种也可以) 求方程y’= y - 2x/y,y(0) = 1 (0≤x≤1,h = 0.1) 的数值解,要求编写程序
是这样吗
function f=rigid(x,y)
f=y-2*x/y;
[x,y]=ode23(‘rigid’,[0,1],0.1);
x(1)=0;y1(0)=1;y2(0)=1;h=0.1;
for k=1:10
x(k+1)=x(k)+h;
y1(k+1)=y1(k)+h*(y(k)-2*x(k)/y(k));
y2(k+1)=y2(k)+h*(y(k+1)-2*x(k+1)/y(k+1));
end
x=0:0.1:1;
plot(x,y1,'k:',x,y2,'r--')
好复杂