solve nonlinear equations
Show older comments
Im trying to solve 5 nonlinear equations with 5 unknowns, but im getting an error?
x0 = [-0.0018; 2.8914; 30.6782; 35.3729; 0.2481]; % Make a starting guess at the solution
options=optimset('Display','iter'); % Option to display output
[x,fval] = fsolve(@moments,x0,options) % Call solver
??? Error using ==> vertcat
CAT arguments dimensions are not consistent.
Error in ==> moments at 2
F = [2*x(4)*x(2)*x(5)*(x(3)^((x(5)-1)/x(5))) + x(1) + 0.000216678;
Error in ==> fsolve at 254
fuser = feval(funfcn{3},x,varargin{:});
Caused by:
Failure in initial user-supplied objective function evaluation. FSOLVE cannot continue.
What is the error in moments.m at line 2? what did i do wrong?
The m file, moments.m is
function F = moments(x)
F = [2*x(4)*x(2)*x(5)*(x(3)^((x(5)-1)/x(5))) + x(1) + 0.000216678;
2*x(2)*x(5)*(x(3)^((x(5)-1)/x(5)))-4*(x(4)^2)*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1)+(2*x(4)*x(2)*x(5)*(x(3)^((x(5)-1)/x(5))) + x(1))^2 - 0.000045035;
-12*x(4)*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1)+8*(x(4)^3)*(x(3)^((x(5)-3)/x(5)))*x(2)*x(5)*(x(5)-1)*(x(5)-2)+3*(2*x(2)*x(5)*(x(3)^((x(5)-1)/x(5)))-4*(x(4)^2)*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1))+(2*x(4)*x(2)*x(5)*(x(3)^((x(5)-1)/x(5))) + x(1))^3 +0.00000020584664 ;
-12*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1)+48*x(4)^2*(x(3)^((x(5)-3)/x(5)))*x(2)*x(5)*(x(5)-1)*(x(5)-2)-16*x(4)^4*(x(3)^((x(5)-4)/x(5)))*x(2)*x(5)*(x(5)-1)*(x(5)-2)*(x(5)-3)+3*(2*x(2)*x(5)*(x(3)^((x(5)-1)/x(5)))-4*(x(4)^2)*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1))^3+4*(-12*x(4)*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1)+8*(x(4)^3)*(x(3)^((x(5)-3)/x(5)))*x(2)*x(5)*(x(5)-1)*(x(5)-2))*(2*x(4)*x(2)*x(5)*(x(3)^((x(5)-1)/x(5))) + x(1))+6*(2*x(2)*x(5)*(x(3)^((x(5)-1)/x(5)))-4*(x(4)^2)*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1))*(2*x(4)*x(2)*x(5)*(x(3)^((x(5)-1)/x(5))) + x(1))^2+(2*x(4)*x(2)*x(5)*(x(3)^((x(5)-1)/x(5))) + x(1))^4 -0.000000017773077721;
(120*x(4)*(x(3)^((x(5)-3)/x(5)))*x(2)*x(5)*(x(5)-1)*(x(5)-2)-160*x(4)^3*(x(3)^((x(5)-4)/x(5)))*x(2)*x(5)*(x(5)-1)*(x(5)-2)*(x(5)-3)+32*x(4)^5*(x(3)^((x(5)-5)/x(5)))*x(2)*x(5)*(x(5)-1)*(x(5)-2)*(x(5)-3)*(x(5)-4))+10*(2*x(2)*x(5)*(x(3)^((x(5)-1)/x(5)))-4*(x(4)^2)*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1))*(-12*x(4)*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1)+8*(x(4)^3)*(x(3)^((x(5)-3)/x(5)))*x(2)*x(5)*(x(5)-1)*(x(5)-2))+5*(-12*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1)+48*x(4)^2*(x(3)^((x(5)-3)/x(5)))*x(2)*x(5)*(x(5)-1)*(x(5)-2)-16*x(4)^4*(x(3)^((x(5)-4)/x(5)))*x(2)*x(5)*(x(5)-1)*(x(5)-2)*(x(5)-3))*(2*x(4)*x(2)*x(5)*(x(3)^((x(5)-1)/x(5))) + x(1))+15*(2*x(2)*x(5)*(x(3)^((x(5)-1)/x(5)))-4*(x(4)^2)*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1))^2*(2*x(4)*x(2)*x(5)*(x(3)^((x(5)-1)/x(5))) + x(1))+10*(-12*x(4)*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1)+8*(x(4)^3)*(x(3)^((x(5)-3)/x(5)))*x(2)*x(5)*(x(5)-1)*(x(5)-2))*(2*x(4)*x(2)*x(5)*(x(3)^((x(5)-1)/x(5))) + x(1))^2+10*(2*x(2)*x(5)*(x(3)^((x(5)-1)/x(5)))-4*(x(4)^2)*(x(3)^((x(5)-2)/x(5)))*x(2)*x(5)*(x(5)-1))*(2*x(4)*x(2)*x(5)*(x(3)^((x(5)-1)/x(5))) + x(1))^3+(2*x(4)*x(2)*x(5)*(x(3)^((x(5)-1)/x(5))) + x(1))^5 + 0.000000000396254743];
Accepted Answer
More Answers (0)
Categories
Find more on Numerical Integration and Differential Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!