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´úÂë¡£