因为你提供的微分方程用dsolve()得不到其解析解,只能用数值分析来求解。可以用matlab提供了微分方程数值求解函数ode45()。求解过程如下:
y0=30;
[x,y] = ode45(@func,[0 10],y0);
plot(x,y)
xlabel('x'),ylabel('y(x)')
legend('y(x)数值解')
func.m 自定义函数文件
function f= func(x,y)
f=(-7/6*1.4*30^(-3)-64*0.0001*30^(-3)+128*0.0001*y^(-3))/(7/6*y^(-3));
end
运行结果