Question

matlab解微分方程，无解析解怎么做图

-l*(l^2+x^2+h^2+y^2)-2*(-l^4-h^4-y^4+3*x^4+2*h^2*y^2-4*l^2*h^2-2*x^2*l^2+2*h^2*x^2+2*y^2*x^2)*Dy=0
这个微分方程在matlab中好像没有解析解（Explicit solution could not be found），怎么用数值方法求出方程的图形
其中 l=0.3 h=20 谢谢
y(0)=20 能在用初始值时假定一个么，主要是方法。。。。

Answer 1

初始值很重要，方法是这样的：

function hh
x0=0;%t0=20,即y(0)=20
x_span=0:10;
[x,y] = ode45(@myodefun,x_span,x0)
plot(x,y)
function dy=myodefun(y,t)%调用Dy
%将Dy变成显式函数
Dy=solve('-l*(l^2+x^2+h^2+y^2)-2*(-l^4-h^4-y^4+3*x^4+2*h^2*y^2-4*l^2*h^2-2*x^2*l^2+2*h^2*x^2+2*y^2*x^2)*Dy','Dy');
l=0.3; h=20;
dys=subs(Dy);
dy=subs(dys,'x','t');%t是matlab的内部自变量
dy=eval(dy);%转换