Hi,
I'm using fsolve in a loop, as an example can we find the first 10 roots of sin(x) and store then in an array.
restart;
x:=array(0..10):
x[0]:=0:
for i from 1 to 10 do
x[i]:= fsolve(sin(y)=0, y=x[i-1]..10 , avoid={y=x[i-1]});
end do;
Firstly this solve only 0,2*Pi, 3*Pi...Why does it skip Pi? also can I change the range to x[i-1]..infinity?
Basically how can I force fsolve to find a solution that is GREATER than the previous solution???
Thanks!
Matt