[Maxima] problems in storing values into an array made with make_array

Robert Gloeckner Robert.Gloeckner at Web.de
Tue Mar 20 11:56:20 CDT 2007


Hello,

/*i would like to put results of a function into a file for later
gnuplotting it.
first step is to evaluate the function without quitting if an error
occurs this is done with "try1d" - the result is either a value or "?".
second step is to put all the values (or the "?") into an array. this
should be done with "testfunc1d". i need this step, because if i would
try to put the values directly into a file using with_stdout, then all
error messages appear in this file, too.
last step would be to put those values into a file.

i fail in doing second step.
below you can see the function definitions and a test run. as you can
see, maxima does not store the value inside the array i created via
make_array, it does store a lisp expression. how can i force evaluation?
*/

/**********************************************************/
try1D(func1D, x) := block([result],
  result : ev( errcatch( func1D(x)), numer),
  if( length(result) < 1) then( "?") else (first(result)));



testfunc1d(func, xmin, xmax, xstep) := block([
  n : ev((xmax - xmin)/xstep, numer), x, result_array],
  print( "n = ", n),
  result_array : make_array( any, ?truncate( 2 * (n + 1))),
  (for i: 0 thru n do (
    x : xmin + xstep * i,
    result_array[ 2 * i ] : x,
    result_array[ 2 * i + 1] : ev(try1D(func, x)))),
result_array);
/**********************************************************/

(%i72) testfunc1d(sin, 2.0, 9.0, 2.0);

n =  3.5
(%o72) {Array:  #(2.0 (($try1D SIMP) %SIN 2.0) 4.0 (($try1D SIMP) %SIN
4.0) 6.#
0
  (($try1D SIMP) %SIN 6.0) 8.0 (($try1D SIMP) %SIN 8.0) NIL)}




thanks for any hint
Robert


More information about the Maxima mailing list