Solving system of nonlinear euations

In mathematica we use the command NSolve[{x^2 + y^3 == 1, 2 x + 3 y == 4}, {x, y}, Reals] without needing any initial guess to solve a system of non-linear equations. I need the same command in Matlab to solve system of non-linear equations BUT without needing the initial guess. THANKS

[Sara]

What's the problem with assigning initial guesses? Can't you make them up with some knowledge of the system you're trying to solve?

[Kashif]

I dont have any initial guess as i am trying to solve a system of 4 equations and 4 four variable, and each equation is containing 24 nonlinear terms. Can you please suggest me some way in Matlab, as i am new in matlab. Thanks

[Sara]

I don't think you can. You'll have to use fsolve, which requires an initial guess. If it's a mathematical exercise, I'd go with an array of zeros; otherwise, if the eqn represent a physical system, you can use some reasonable values (e.g. mass is usually >= 0).

[Kashif]

Actually to use fsolve we need to make another m file for the function, in my case ceofficients of variables in equations depend on many other values in my program . Can we make function which is called by fsolve within the program.

Answer 1

[Shoaibur Rahman]

syms x y
[x,y] = vpasolve(x^2 + y^3 == 1, 2*x + 3*y == 4, x,y)

[Kashif]

I have used the above command but i am getting this error Expression or statement is incorrect--possibly unbalanced (, {, or [.

My Codes were the following

syms x y z s [x,y,z,s] = vpasolve(Aphi(n/8).*x+Bphi(n/8).*y+Aphi(n/8).*z+Bphi(n/8).*s^2+Ephi(n/8).*x^2+Ephi(n/8).*y^2+Ephi(n/8).*z^2+Ephi(n/8).*s^2+Iphi(n/8).*x*z+Jphi(n/8).*y*s+Kphi(n/8).*x*s+Kphi(n/8).*y*s+Mphi(n/8).*y*z^2 +Mphi(n/8).*y*s^2+Mphi(n/8).*s*x^2+Mphi(n/8).*z*y^2+Ophi(n/8).*x*z^2+Ophi(n/8).*x*s^2+Ophi(n/8).*z*x^2+Ophi(n/8).*z*y^2-Uphi(n/8).*x^2*z^2-Uphi(n/8).*x^2*s^2-Uphi(n/8).*y^2*z^2-Uphi(n/8).*y^2*s^2 ==Vphi(n/8), Aphi(n+n/8).*x+Bphi(n+n/8).*y+Aphi(n+n/8).*z+Bphi(n+n/8).*s^2+Ephi(n+n/8).*x^2+Ephi(n+n/8).*y^2+Ephi(n+n/8).*z^2+Ephi(n+n/8).*s^2+Iphi(n+n/8).*x*z+Jphi(n+n/8).*y*s+Kphi(n+n/8).*x*s+Kphi(n+n/8).*y*s+Mphi(n+n/8).*y*z^2 +Mphi(n+n/8).*y*s^2+Mphi(n+n/8).*s*x^2+Mphi(n+n/8).*z*y^2+Ophi(n+n/8).*x*z^2+Ophi(n+n/8).*x*s^2+Ophi(n+n/8).*z*x^2+Ophi(n+n/8).*z*y^2-Uphi(n+n/8).*x^2*z^2-Uphi(n+n/8).*x^2*s^2-Uphi(n+n/8).*y^2*z^2-Uphi(n+n/8).*y^2*s^2 ==Vphi(n+n/8), Aphi(2*n+n/8).*x+Bphi(2*n+n/8).*y+Aphi(2*n+n/8).*z+Bphi(2*n+n/8).*s^2+Ephi(2*n+n/8).*x^2+Ephi(2*n+n/8).*y^2+Ephi(2*n+n/8).*z^2+Ephi(2*n+n/8).*s^2+Iphi(2*n+n/8).*x*z+Jphi(2*n+n/8).*y*s+Kphi(2*n+n/8).*x*s+Kphi(2*n+n/8).*y*s+Mphi(2*n+n/8).*y*z^2 +Mphi(2*n+n/8).*y*s^2+Mphi(2*n+n/8).*s*x^2+Mphi(2*n+n/8).*z*y^2+Ophi(2*n+n/8).*x*z^2+Ophi(2*n+n/8).*x*s^2+Ophi(2*n+n/8).*z*x^2+Ophi(2*n+n/8).*z*y^2-Uphi(2*n+n/8).*x^2*z^2-Uphi(2*n+n/8).*x^2*s^2-Uphi(2*n+n/8).*y^2*z^2-Uphi(2*n+n/8).*y^2*s^2 ==Vphi(2*n+n/8), Aphi(n/16).*x+Bphi(n/16).*y+Aphi(n/16).*z+Bphi(n/16).*s^2+Ephi(n/16).*x^2+Ephi(n/16).*y^2+Ephi(n/16).*z^2+Ephi(n/16).*s^2+Iphi(n/16).*x*z+Jphi(n/16).*y*s+Kphi(n/16).*x*s+Kphi(n/16).*y*s+Mphi(n/16).*y*z^2 +Mphi(n/16).*y*s^2+Mphi(n/16).*s*x^2+Mphi(n/16).*z*y^2+Ophi(n/16).*x*z^2+Ophi(n/16).*x*s^2+Ophi(n/16).*z*x^2+Ophi(n/16).*z*y^2-Uphi(n/16).*x^2*z^2-Uphi(n/16).*x^2*s^2-Uphi(n/16).*y^2*z^2-Uphi(n/16).*y^2*s^2 ==Vphi(n/16), , x,y,z,s)

[Shoaibur Rahman]

I see an extra comma (,) in the very last part of the code. ,,x,y,z,s --- remove one comma. It also seems that you have to define Aphi, Bphi........ with n