Hello everyone,

I am trying to solve numerically int( f(t,z) , t=0..T ) = 0 , in z for a cumbersome f.

I tried z1=fsolve( int( f(t,z) , t=0..T ) = 0 , z). But then I tried int( f(t,z1) , t=0..T ) and the result is clearly not zero nor anything small.

It looks like Maple evaluates analytically the integral, and does it wrong (check this for more details) so fsolve uses the wrong equations.

Anyone knows how I can force Maple to evaluate numerically the integral at each step of the fsolve function?

Thank you!

Hello everyone

I have a problem with the code below in maple:

> restart;
> f := signum(a*cos(t+b)+sin(t+c));
> res := int(f, t = 0 .. 2*Pi);
                       2 Pi signum(a cos(b) + sin(c))
> a := 0.01; b := 0; c := 0;
> res; int(f, t = 0 .. 2*Pi);
                                    2 Pi
                                     0.

As you can see, results are not consistent with one another...

If I set a := 0.01; b := 1.35; c := 4.4; It is even worse, the computation of int(f, t = 0 .. 2*Pi) lasts more than 10 minutes (when I interrupt computation)

I tested it in Maple 12 and Maple 2015 and I get exactly the same results. I also replaced signum with csgn and I get the same problem.

I really need your help!!