1+1;

%o1;

%;

2+3/4*5;

2^1000;

100!;

set_display('ascii)$

3^1000;

2.0^1000;

set_display('xml)$

sqrt(2);

sqrt(2.0);

float(sqrt(2));

if numer#false then numer:false else numer:true;

sqrt(2);

if numer#false then numer:false else numer:true;

fpprec : 32;

bfloat(sqrt(2));

sin(4);

%pi;

sin(%pi/4);

sqrt(4);

f(x):=1/(1+x^2);

f(2);

f(2.0);

f(d);

plot2d([f(x)], [x,-5,5],
 [plot_format, gnuplot],
 [gnuplot_preamble, "set grid;"])$

plot3d(1/(1+x^2+y^2), [x,-5,5], [y,-5,5], [plot_format,gnuplot],
 [gnuplot_preamble, "set pm3d map; unset surf"])$

plot3d(1/(1+x^2+y^2), [x,-5,5], [y,-5,5], [plot_format,gnuplot],
 [gnuplot_preamble, "set pm3d at s; unset surf; unset colorbox"])$

100!;

factor(100!);

factor(1-x^4);

ratsimp(%);

expr:(1-x^4)/(1-x^3);

ratsimp(expr);

gcd(1-x^4,1-x^3);

f(x);

integrate(f(x), x, 0, 1);

romberg(f(x), x, 0, 1);

integrate(f(x),x);

diff(f(x),x);

diff(f(x),x,3);

taylor(f(x), x, 0, 8);

p(x):=1-x^2+x^4;

plot2d([f(x),p(x)], [x,-1,1],
 [plot_format, gnuplot],
 [gnuplot_preamble, "set zeroaxis;"])$

niceindices(powerseries(f(x), x, 0));

sum(1/k^2, k, 1, inf), simpsum;

sum(1/k^2, k, 1, inf);

sum(1/k^3, k, 1, inf), simpsum;

sum(k,k,1,n),simpsum;

factor(%);

sum(k^2,k,1,n),simpsum;

v:[1,2,3];

A:matrix([1,2,3],[4,5,6]);

v.A;

A.v;

transpose(A);

v.transpose(A);

w:makelist(k^2,k,1,10);

M: genmatrix(lambda([i,j], i^2+j^3), 4, 4);

h[i,j]:=i^2+j^3;

H:genmatrix(h,4,4);

determinant(H);

invert(H);

H^^-1;

H^-1;

2+H;

H.H;

H^^2;

H*H;

H^2;

A;

A.A;

A*A;

eigenvalues(H);

expand(charpoly(H,x));

eigenvectors(H);

%[1];

%o86[2];

linsolve([2*x+y=3, 2*y+z=3, x+2*z=3], [x,y,z]);

echelon(H);

triangularize(H);

set_display('none)$

H;

for k:1 thru 20 do print(k^2);

for k:1 step 2 while k<20 do print(k^2);