piecewise a possibility?

Preben Alsholm — Sat, 05 Nov 2011 19:48:35 Z

You don't give all the code for norma, but the problem is exhibited in this simple example:

f:=x->if x<Pi/2 then sin(x)^2 else cos(x) end if;
int(f(x),x=0..2*Pi);

#In this simple case using piecewise instead is a good idea

fp:=x->piecewise(x<Pi/2,sin(x)^2,cos(x));
int(fp(x),x=0..2*Pi);

Delay evaluation

Joe Riel — Sat, 05 Nov 2011 19:52:42 Z

One way is to delay the evaluation of the function call norma by quoting it with single forward quotes:

int( 'norma'(...), ..., numeric )