如何使用matlab用梯形法解常微分方程dy/dx=y/x-y^2/x，y(1)=0.5，求出其数值解
1个回答
展开全部
解析解：
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')
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')
本回答由提问者推荐
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
