|
|
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
mingdong(金币+3,VIP+0):感谢交流！ 11-11 15:35
mltgood(金币+10,VIP+0):非常感谢您的帮助 11-11 19:28
function sdf
global a
a=66;
[t,y]=ode45(@myfun,[0,100-a],0.1);
[t y]
plot(t,y)
function f=myfun(t,x)
global a
f=(1/(x+a)+log(x+a)/(x+a)+log(x+a)^2/(x+a)+log(x+a)^3/(x+a))*x^3; |
|
回复此楼
