Advanced Search

Question: How to solve algebraic equations?

Question:How to solve algebraic equations?

torabi 90 Maple 2015
solve fsolve
April 06 2016
0

hi.please help me for solve algebric equations below?

solve or fsolve dos not any answer?

i think use Newton  Raphson, because solve or fsolve not work.

thankscomparision.mw

restart; mu1 := .1; mu2 := .1; sigma1 := -40; F := 25; upsilon1 := (1/2)*sigma2; upsilon2 := (3/2)*sigma2-sigma1; sigma2 := 100

100

 (1)

gamma11 := 2.686901; -1; gamma12 := 7.175339; -1; gamma21 := 2.436735; -1; gamma22 := 12.94855; -1; gamma22 := 12.94855; -1; delta1 := .928207; -1; delta2 := .105073; -1; s11 := .629894; -1; s12 := .217601; -1; s21 := 0.73897e-1; -1; s22 := .805815

.805815

 (2)

 

 

 

 
Q1 := -mu1*p1-upsilon1*q1+gamma11*q1*(p1^2+q1^2)+gamma12*q1*(p2^2+q2^2)-delta1*(2*p1*q1*p2-q2*(p1^2-q1^2))-s11*F*q1+s12*F*q2 = 0

-.1*p1-65.747350*q1+2.686901*q1*(p1^2+q1^2)+7.175339*q1*(p2^2+q2^2)-1.856414*p1*q1*p2+.928207*q2*(p1^2-q1^2)+5.440025*q2 = 0

 (3)

Q2 := -mu1*q1+upsilon1*p1-gamma11*p1*(p1^2+q1^2)-gamma12*p1*(p2^2+q2^2)-delta1*(2*p1*q1*q2+p2*(p1^2-q1^2))-s11*F*p1-s12*F*p2 = 0

-.1*q1+34.252650*p1-2.686901*p1*(p1^2+q1^2)-7.175339*p1*(p2^2+q2^2)-1.856414*p1*q1*q2-.928207*p2*(p1^2-q1^2)-5.440025*p2 = 0

 (4)

Q3 := -mu2*p2-upsilon2*q2+gamma21*q2*(p1^2+q1^2)+gamma22*q2*(p2^2+q2^2)+delta2*(3*p1^2*q1-q1^3)+s21*F*q1 = 0

-.1*p2-190*q2+2.436735*q2*(p1^2+q1^2)+12.94855*q2*(p2^2+q2^2)+.315219*p1^2*q1-.105073*q1^3+1.847425*q1 = 0

 (5)

Q4 := -mu2*q2+upsilon2*p2-gamma21*p2*(p1^2+q1^2)-gamma22*p2*(p2^2+q2^2)+delta2*(-p1^3+3*p1*q1^2)-s21*F*p1 = 0

-.1*q2+190*p2-2.436735*p2*(p1^2+q1^2)-12.94855*p2*(p2^2+q2^2)-.105073*p1^3+.315219*q1^2*p1-1.847425*p1 = 0

 (6)

 

-.1*q2+((3/2)*sigma2-40)*p2-2.436735*p2*(p1^2+q1^2)-12.94855*p2*(p2^2+q2^2)-.105073*p1^3+.315219*q1^2*p1-1.847425*p1

 (7)

NULL

fsolve({Q1, Q2, Q3, Q4}, {p1, p2, q1, q2})

fsolve({-.1*p1-65.747350*q1+2.686901*q1*(p1^2+q1^2)+7.175339*q1*(p2^2+q2^2)-1.856414*p1*q1*p2+.928207*q2*(p1^2-q1^2)+5.440025*q2 = 0, -.1*p2-190*q2+2.436735*q2*(p1^2+q1^2)+12.94855*q2*(p2^2+q2^2)+.315219*p1^2*q1-.105073*q1^3+1.847425*q1 = 0, -.1*q1+34.252650*p1-2.686901*p1*(p1^2+q1^2)-7.175339*p1*(p2^2+q2^2)-1.856414*p1*q1*q2-.928207*p2*(p1^2-q1^2)-5.440025*p2 = 0, -.1*q2+190*p2-2.436735*p2*(p1^2+q1^2)-12.94855*p2*(p2^2+q2^2)-.105073*p1^3+.315219*q1^2*p1-1.847425*p1 = 0}, {p1, p2, q1, q2})

 (8)

solve(Q1, Q2, Q3, Q4)

Error, invalid input: too many and/or wrong type of arguments passed to solve; first unused argument is -.1*q1+34.252650*p1-2.686901*p1*(p1^2+q1^2)-7.175339*p1*(p2^2+q2^2)-1.856414*p1*q1*q2-.928207*p2*(p1^2-q1^2)-5.440025*p2 = 0

  

Sol := [fsolve(ZZ, omega)]



Download comparision.mw
Please Wait...