Ploting a function and its roots

f=@(x)x.^2-5*sin(x);
[p1,p2]=fzero(f,1,[-1,3])

fplot(f1,[-1,3])

Can I "draw" the roots p1,p2 in the same function fplot?

Answer 1

[Star Strider]

The fzero function will only return one value for each initial estimate. You need to provide different initial estimates in different calls to it, if you want other roots. (For several roots, this usually requires a for loop.) Plotting the results is then straightforward.

Try this:

f=@(x)x.^2-5*sin(x);
p1 = fzero(f,1);
p2 = fzero(f,5);

fplot(f,[-1,3])
hold on
plot([p1,p2], [0 0], 'p')
hold off

[Harel Harel Shattenstein]

I can use fsolve too?

[Star Strider]

Yes. However it is subject to the same initial estimate constraints, so that different initial estimates can result in different zeros and different parameter estimates, at least in my experience. The difference is that fsolve (link) finds parameters defining roots of systems of functions, and is not limited to single functions, as is fzero.