Figure 2.1.1
# number sign # introduces a comment
# statements must be ended by a semi-colon
# or by a colon (surpresses printing)
# but can span multiple lines
p1:=plot([C,(9/5)*C+32,C=0..100],-10..100,-30..220,tickmarks=[5,0]):
p2:=textplot({[-5,32,`32`],[-5,212,`212`]}):
p3:=textplot({[-5,120,`Temperature F`],[60,-30,`Temperature C`]}, align=LEFT):
with(plots):
interface(plotdevice=postscript,plotoutput=`fig2.1.1.ps`);
display({p1,p2,p3});



      Page 13
assume(a>1); assume(k>0);
limit(x^k/a^x,x=infinity);

sol:=solve(x^3=2^x,x);
evalf({sol[1],sol[2]});

     Page 14
plot({[x,x^3,x=0..12],[x,2^x,x=0..12],
  [x,100*x^3/2^x,x=0..14]},x=0..14,y=0..4000);




       Exercises/Experiments

1A.
# Change of units is built-in
# type: ?convert
convert(x*inches,metric);
convert(y*pounds/gallon,metric,US);
convert(z*miles/hour,metric);

2B.
plot({Pi*r^2,4/3*Pi*r^3},r=0..1);
plots[loglogplot]({ Pi*r^2,4/3*Pi*r^3},r=0.1..1);
plot({ 3*5^x,5*3^x},x=0..1);
plots[logplot]({ 3*5^x,5*3^x},x=0..1);



3.
plot({ 3*x^2+5*x+7,2^x},x=0..7);
plot((3*x^2+5*x+7)/2^x,x=0..10,y=0..10);
limit((3*x^2+5*x+7)/2^x,x=infinity);

4.
eq:=diff(y(t),t)=3*y(t);
sol:=dsolve({ eq,y(0)=2},y(t));
y:=unapply(rhs(sol),t); plot(y(t),t=0..1);
y:='y';