>

norm3d:=proc(a,b,c) sqrt(a^2+b^2+c^2) end;
norm3d(3,4,0);

>

i:=`i`:total:=0:
for i from 0 to 100 do
total:=i+total:
end do:
total;

>

p:=4:
if p<2 then printf("p is less than 2");
elif p>=2 then printf("p is not less than 2");
end if;

p is not less than 2

>

F:=array(1..10000):
F[1]:=1:F[2]:=1:N:=15:
for i from 3 to N do
F[i]:=F[i-1]+F[i-2]:
end do:
seq(F[i],i=1..N);

>

mu:=3.2:x[0]:=0.2:
for n from 0 to 99 do
x[n+1]:=mu*x[n]*(1-x[n]):
end do:
for n from 90 to 99 do
nprintf("x[%-d]=%g",n+1,x[n+1]);
end do;

>

with(plots):
p1:=plot(x^2,x=-4..4,color=blue):p2:=plot(4-x^2,x=-4..4,color=red):
t1:=textplot([2.6,6,`y=x^2`],align=RIGHT):
t2:=textplot([-3,-6,`y=4-x^2`],align=RIGHT):
display({p1,p2,t1,t2},font=[TIMES,ROMAN,15],tickmarks=[3,3]);

>

restart:with(DEtools):with(plots):
deqn1:=diff(x(t),t$2)=-2*diff(x(t),t)-25*x(t):
p3:=DEplot(deqn1,x(t),t=0..10,[[x(0)=1,D(x)(0)=0]],stepsize=0.1,linecolor=blue,thickness=1):
deqn2:=diff(x(t),t$2)=-25*x(t):
p4:=DEplot(deqn2,x(t),t=0..10,[[x(0)=1,D(x)(0)=0]],stepsize=0.1,linecolor=red,thickness=1):
t3:=textplot([6,1,`Harmonic motion`],align=RIGHT,color=red):
t4:=textplot([1.8,0.2,`Damped motion`],align=RIGHT,color=blue):
display({p3,p4,t3,t4},labels=[`t`,`x`]);

>

restart:with(DEtools):with(plots):
interactiveparams(plot,[{1-a*x,(x-1)/a},x=0..2],a=0..1);

>

sys1:=diff(x(t),t)=-a*x(t),diff(y(t),t)=a*x(t)-b*y(t),diff(z(t),t)=b*y(t);

>

dsolve([sys1,x(0)=10,y(0)=0,z(0)=0]);

>

interactiveparams(plot,[{10*exp(-a*t),10*a*exp(-b*t)/(a-b)-10*a*exp(-a*t)/(a-b),(10*exp(-a*t)*b-10*exp(-b*t)*a^2/(a-b)+10*exp(-b*t)*a*b/(a-b)+10*a-10*b)/(a-b)},t=0..30],a=0.1..0.9,b=0.1..0.8);