solinit = bvpinit(linspace(0,1,10),[-1 1]);
f1=@(x,y)[y(2);1100*(y(2)-0.5*(1-y(1)).^2)];
f2=@(ya,yb)[ya(1);yb(2)];
sol = bvp4c(f1,f2,solinit);
xint = linspace(0,1,50);
yint = deval(sol,xint);
plot(xint,yint(1,:),'-ro');

2楼: Originally posted by 凌风7388 at 2012-06-30 11:23:02
给出matlab代码，谢谢
13/72/759985_1341026557_410.png
求方程的解，给出matlab代码。