diff(x*f(x),x);
f(x):=sin(x);
ev(%th(2),diff)

x;
x:3;
x;
'x

f(x):=x^2;
'f(2);
ev(%,f);
'(f(2));
f(2);
sum(i!,i,1,4);
'sum(i!,i,1,4);
remvalue(x);
'integrate(f(x),x,a,b)

for i thru 5 do s:i^2+s;
s;
ev(%,s:0);
ev(%th(2))

'sum(g(i),i,0,n);
z*%e^z;
ev(%,z:x^2);
subst(x^2,z,%th(3));
a:%;
1+a;
kill(a,y);
a

integrate(y^2,y)

f(y):=diff(y*log(y),y,2)

f(y):=1/y

(y+x)^3;
diff(%,x)

1+x = y^2;
x-1 = 1+2*y;
%+%th(2);
%th(3)/y;
1/%

(sqrt(2.25)+sqrt(-4))^2;
expand(%)

expand(sqrt(2*%i))

a[n]:=n*a[n-1];
a[0]:1;
a[5];
a[n]:=n;
a[6];
a[4]

f(x):=y+x^2;
f(2);
ev(f(2),y:7)

f(x):=1+sin(x)^2;
f(1+x);
g(y,z):=3*y+f(z);
ev(g(z+2*y,-0.5),y:7);
functions

h(n):=sum(i*x^i,i,0,n)

t[n](x):=ratexpand(2*x*t[n-1](x)-t[n-2](x));
t[0](x):=1;
t[1](x):=x;
t[4](y)

g[n](x):=sum(ev(x),i,n,2+n);
h(n,x):=sum(ev(x),i,n,2+n);
g[2](i^2);
h(2,i^2)

p[n](x):=ratsimp(diff((x^2-1)^n,x,n)/(2^n*n!));
q(n,x):=ratsimp(diff((x^2-1)^n,x,n)/(2^n*n!));
p[2];
p[2](1+y);
q(2,y);
p[2](5)

f[i,j](x,y):=y^j+x^i;
g(fun,a,b):=print(fun,"
applied to ",a,"
and ",b,"
is ",fun(a,b));
g(f[2,1],sin(%pi),2*c)