Hello, I have an oscillatory decaying equation f(h), when it intersects with zero, it should give me multiple values, but vpasolve or fsolve is only giving one value close to the initial value I provided, but I prefer to get all the solution in the range from [25, 50]. The equation is below: syms h P0=4.181323907591596e+03; d1=13.5; d=9.734259414551886; d2=4.706483519517566; F=(2*pi*P0*sin(2*pi*(h-d1)/d1)*exp(d^3/(d1^2*d2)-(h-d1)/d2)/(d1*10e-9)+… P0*cos(2*pi*(h-d1)/d1)*exp(d^3/(d1^2*d2)-(h-d1)/d2)/(d2*10e-9))==0; S=vpasolve(F,h,[25,50])
If you can tell me where I went wrong, that would be great, thanks. PS: You can try to change h from [10,50] to see the intersected values, but I am not able to obtain each one of them numerically.
Best Answer