My problem was that it takes long time to guess the initial value for fsolve.
B ut, now that you have given the idea , i tried the following code :
options = optimset('Display','off','TolFun',eps);
>> [op,val] = fsolve(fh, [ 1,1] ,options )
Yes. It gave the correct result as : 445.5183 & 845.9731
Thank you very much for your advice (taking eps as the tolerance ). I hope it works for other problems also.
Thank you very much for your advice.
With best regards,
kalpana
B ut, now that you have given the idea , i tried the following code :
options = optimset('Display','off','TolFun',eps);
>> [op,val] = fsolve(fh, [ 1,1] ,options )
Yes. It gave the correct result as : 445.5183 & 845.9731
Thank you very much for your advice (taking eps as the tolerance ). I hope it works for other problems also.
Thank you very much for your advice.
With best regards,
kalpana