Can I set the initial value of the solution in "solve" function? Like "fsolve" function?

Can I set the initial value of the solution in "solve" function? Like "fsolve" function?

My code is:

syms p z

b=150;beta=10;s=10;s0=5;alpha=0.5;c1=30;cv=10;N=150;

f2=@(x)((x-z).*normpdf(x,10,40));

eqns=[p-(b+beta.*(s-s0)+10+alpha.*c1)./(2.*alpha)+(int(@(x) f2(x),z,N))./(2.*alpha)==0,p-c1-(p-cv).*cdf('Normal',z,10,40)==0];

S=solve(eqns,[p,z]);

S.p

S.z

The result is [p z]=[30,-185], but I know that there exists the other solution [223, 63].

In order to find the latter solution, can I set the initial value of the solution in "solve" function? Like "fsolve" function? If the answer is yes,how?

Thank you very much!

Answer 1

[Star Strider]

The vpasolve function allows initial estimates and ranges, See: Find Multiple Solutions by Specifying Initial Guesses and related sections for more information.

[qingwei wang]

Thank you!

[Star Strider]

As always, my pleasure!