Question

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--')

Answer 1

好复杂

Answer 2

kgfiospf