read("/p/edu/mat-1.414/maple/v201.mpl"):
Tässä on toimiva SD-askel komento komennolta. Hiukan piti kokeilla ja sen tuloksena lisätä mm. evalm. Tästä pitäisi olla helppo jatkaa. (Pari yllättävää kaatumista tapahtui. Muistakaa, että CTR-S tallettaa, kannattaa viljellä usein.)
Tässä pidetään rinnakkain
jonoa x0 ja Vectoria x0v x1 x1v ....Funktiota f voidaan soveltaa suoraan jonoon, kuten x0, x1 ...
Jos vektorilasku ei mene komponenteille, täytyy komentaa vanhakantaiseen tyyliin evalm();
f:=(x,y)->4*x^2-4*x*y+2*y^2; g:=[D[1](f),D[2](f)]; x0:=1,1: x0v:=; u:=-Vector(g(x0)); jana:=v2l(evalm(x0v+t*u)); phi:=simplify(f(op(jana))); plot(phi,t=0..1); t1:=solve(diff(phi,t)=0,t); x1v:=x0v+t1*u;x1:=op(v2l(x1v)); f(x1),subs(t=t1,phi);