Luento ti 17.10.00
EHY vakioiden variointikaavalla
> with(linalg):
> A:=matrix([[-2,1],[1,-2]]);g:=t->vector([2*exp(-t),3*exp(t)]);
> psi:=t->matrix([[exp(-3*t),exp(-t)],[-exp(-3*t),exp(-t)]]);
> inverse(psi(s));evalm(% &* g(s));
> map(int,%,s);
> evalm(psi(s)&*%);
> map(simplify,%);
> x:=evalm(psi(s)&*vector([c1,c2])+%);
> map(simplify,%);
> vasen:=map(diff,x,s);
> oikea:=evalm((A&*x)+g(s));
> evalm(vasen-oikea);
>