A:matrix([3,1,1,0],[1,2,1,1],[1,1,2,1],[0,1,1,3]);

eigenvectors(A);

P:transpose(matrix([1,0,0,-1],[1,1,1,1],[1,0,-2,1],[0,1,-1,0]));

An=P.matrix([3^n,0,0,0],[0,5^n,0,0],[0,0,1^n,0],[0,0,0,1^n]).invert(P);

ode2('diff(y,t)+(1/(2*t))*y=2,y,t);

atvalue(x(t),t=0,1);atvalue(y(t),t=0,3);desolve(['diff(x(t),t)=x(t)+2*y(t)-1,'diff(y(t),t)=-x(t)+4*y(t)+t],[x(t),y(t)]);

wxplot2d([(50*%e^(3*t))/9-(11*%e^(2*t))/2+t/3+17/18,(50*%e^(3*t))/9-(11*%e^(2*t))/4-t/6+7/36],[t,0,1]);

taylor(asin(%e^x-1)-%e^(asin(x)),x,0,6);

ds(a,b):=block([c,d],
if (a<b) then (c:a,a:b,b:c),d:1,while (not(d=0)) do (d:a-b, a:b,b:d,if (a<b) then (c:a,a:b,b:c)),
a
);

ds(45,60);