fsolve fails when eqs are 2D integrations

> restart;

> infolevel[all]:=0;

> infolevel[evalf]:=0;

> infolevel[fsolve]:=3;

> e1 := b=evalf(Int(exp(sin(a*x)),x=0..1));

> e2 := a=evalf(Int(exp(sin(2*b*x)),x=0..1));

> fsolve({e1,e2}); #should be 2 variables, why maple recognized as three?Error, (in fsolve) number of equations, 2, does not match number of variables, 3

> fsolve({e1,e2},{a,b}); #works ok!fsolve/sysnewton: trying multivariate Newton iterationfsolve/sysnewton: guess vector [6.5647247568434774, -9.1035796650957585]fsolve/sysnewton: norm of errors: 15.667670001840765fsolve/sysnewton: new norm: 1.6754775047232522 fsolve/sysnewton: iter = 1 |incr| = 16.247 new values a = .62914 b = 1.2079fsolve/sysnewton: new norm: .59028056918433893 fsolve/sysnewton: iter = 2 |incr| = 2.4059 new values a = 1.9642 b = 2.2787fsolve/sysnewton: new norm: 0.42572858304347783e-1 fsolve/sysnewton: iter = 3 |incr| = .53511 new values a = 1.6596 b = 2.0482fsolve/sysnewton: new norm: 0.266378329912618e-3 fsolve/sysnewton: iter = 4 |incr| = 0.46156e-1 new values a = 1.6833 b = 2.0257fsolve/sysnewton: new norm: 0.10369039257e-7 fsolve/sysnewton: iter = 5 |incr| = 0.27633e-3 new values a = 1.6835 b = 2.0256fsolve/sysnewton: new norm: 0.57e-16 fsolve/sysnewton: iter = 6 |incr| = 0.11512e-7 new values a = 1.6835 b = 2.0256

>

>

> e3 := b=evalf(Int(exp(sin(a*x+3*b*y)),[x=0..1,y=0..1],method = _cuhre));

e4 := a=evalf(Int(exp(sin(2*b*x+4*a*y)),[x=0..1,y=0..1],method = _cuhre));

fsolve({e3,e4},{a,b}); #*********maple connot work out in long times.***********fsolve/sysnewton: trying multivariate Newton iterationfsolve/sysnewton: guess vector [6.5647247568434774, -9.1035796650957585]Warning, computation interrupted

evalf(Int(exp(sin(6.5647247568434774*x+3*(-9.1035796650957585)*y)),[x=0..1,y=0..1],method = _cuhre)); # a single integration can be evaluated quickly

This post generated using the online HTML conversion tool
Download the original worksheet