share/calculus/qual.dem

   1 ("A demonstration of the qualitative analysis functions described in file qual.usg.")$
   2
   3 load("qual");
   4
   5 ("Consider the following expression:")$
   6 temp: tanh(cos(x/3)**2-cos(x/3-%pi/2)**2) + sin(x/5)*erf(tan(x)) + atan(y);
   7 ("To do a qualitative analysis of this expression, we write:")$
   8 qual(temp);
   9 ("We may use the ASSUME function to establish non-default bounds on indeterminates.  For example:")$
  10 assume(p>0, q>=-1/4, q<=1/4, b<=bmax)$
  11 ("We may also individually invoke the various property analysis functions invoked by QUAL.  For example:")$
  12 temp: -asin(q+1/2) - q**2 + p*q - 4*q
  13       + log(p) - p**2 + 3*p + sqrt(p) - b**2 + b - 18;
  14 curvature(temp);
  15 temp: expand((x-1)**4 + (x-1)**2 + (y-2)**2 + 5);
  16 stationarypoints(temp);
  17 temp: expand((log(x-3)-5)*(x-y+8)/tan(x-y+2));
  18 zerosandsingularities(temp);
  19 ("QUAL also automatically reveals successive levels of the structure of lengthier expressions.  For example, the second derivative of TEMP with respect to X is about 2 pages long, but:")$
  20 qual(diff(temp,x,2))$