doc/info/figures/colnew-ex2.mac

   1 /* colnew-ex2.mac
   2
   3   Generate figure for colnew example 1
   4
   5 */
   6
   7 load("colnew.mac")$
   8 load("prob2.mac")$
   9
  10 /* Function to generate list of n+1 equally spaced values */
  11 xlist(xmin,xmax,n):=block([dx:(xmax-xmin)/n],makelist(i,i,0,n)*dx+xmin)$
  12 X: xlist(aleft,aright,500)$
  13
  14 /* First run: No initial approx provided */
  15 ipar[9] : 0;
  16
  17 [iflag, fspace, ispace] :
  18 colnew_expert(ncomp, m, aleft, aright, zeta, ipar, ltol, tol,
  19  fixpnt, ispace, fspace,
  20  0, fsub, dfsub, gsub, dgsub, dummy)$
  21
  22 /* Generate solution values */
  23 ans:colnew_appsln(X,4,fspace,ispace)$
  24
  25 /* List of first variable */
  26 z1:maplist(first,ans)$
  27
  28 /* Second run: initial approx provided */
  29 ipar[9] : 1;
  30
  31 [iflag, fspace, ispace] :
  32 colnew_expert(ncomp, m, aleft, aright, zeta, ipar, ltol, tol,
  33  fixpnt, ispace, fspace,
  34  0, fsub, dfsub, gsub, dgsub, solutn)$
  35
  36 /* Generate solution values */
  37 ans:colnew_appsln(X,4,fspace,ispace)$
  38
  39 /* List of first variable */
  40 z2:maplist(first,ans)$
  41
  42 plot2d([[discrete,X,z1],[discrete,X,z2]],
  43   [legend,"Solution I","Solution II"],
  44   [xlabel, "x"],
  45   [ylabel, "φ(x)"],
  46   [pdf_file,"./figures/colnew-ex2.pdf"]);
  47
  48 plot2d([[discrete,X,z1],[discrete,X,z2]],
  49   [legend,"Solution I","Solution II"],
  50   [xlabel, "x"],
  51   [ylabel, "φ(x)"],
  52   [png_file,"./figures/colnew-ex2.png"]);