如何使用matlab用梯形法解常微分方程dy/dx=y/x-y^2/x,y(1)=0.5,求出其数值解

jimmyol 1年前 已收到1个回答 举报

泡沫小子 幼苗

共回答了14个问题采纳率:78.6% 举报

y=dsolve('Dy=y/x-y^2/x','y(1)=0.5','x')
y =
1/(1/x + 1)
梯形公式:
y(1)=0.5;x(1)=1;h=0.1;
for n=1:99
x(n+1)=x(n)+h;
y(n+1)=y(n)+h*( y(n)/x(n)-y(n)^2/x(n));
y(n+1)=y(n)+h/2*( y(n)/x(n)-y(n)^2/x(n) +y(n+1)/x(n+1)-y(n+1)^2/x(n+1) );
end
plot(x,y,'y')

1年前

8
可能相似的问题
Copyright © 2024 YULUCN.COM - 雨露学习互助 - 17 q. 0.033 s. - webmaster@yulucn.com