1
$\begingroup$ 
vpw = 3.3
Mp = 80
rp = 0.014*Mp
b0 = 26.31016273 - (1.7168*10^3/t) - 3.5519*Log[t]
b1 = 24.65681838 - (1.547*10^3/t) - 3.4314*Log[t]
b2 = 9.080370819 - (6.9445*10^2/t) - 1.2222*Log[t]
g = b0 + b1*phip + b2*phip^2
dg = b1 + 2*b2*phip
ddg = 2*b2
dddg = 0

S1 = (1/(phip*rp*vpw)) + (1/(1 - phip)) - (2*g - 2*(1 - 2*phip)*dg - phip*(1 - phip)*ddg)
S2 = -(1/(phip^2*rp*vpw)) + (1/(1 - phip)^2) - (6*dg - 3*(1 - 2*phip)*ddg - 
      phip*(1 - phip)*dddg)

I have tried FindRoot to solve S1==0 and S2==0 but it gives only one root which may be real or imaginary depending upon initial guess.

In order to get all real roots, I am trying to solve S1==0 and S2==0 by using NSolve. But it takes too long time and does not give any solution.

In order to get all possible real roots of above non-linear equations, please suggest me the possible solution of my problem. Thanks in advance and Regards

Improve this question
$\endgroup$
2
  • $\begingroup$ Just to clarify: You want to solve for t and phip? $\endgroup$
    – sebhofer
    Commented Jul 28, 2015 at 12:36
  • $\begingroup$ Here goes a numerical solution for you {phip -> -1.2873402582429818, t -> 414.62463948474345}. As far as I can see there is no real root for t<0. Also using Plot3D I doubt if there is any solution to {S1==0, S2==0} for t>0 and phip>1. Function is also unbounded in that region. $\endgroup$
    – PlatoManiac
    Commented Jul 28, 2015 at 12:51

1 Answer 1

Reset to default
1
$\begingroup$ 
vpw = 3.3; 
Mp = 80; 
rp = 0.014*Mp; 
b0 = 26.31016273 - 1.7168*(10^3/t) - 3.5519*Log[t]; 
b1 = 24.65681838 - 1.547*(10^3/t) - 3.4314*Log[t]; 
b2 = 9.080370819 - 6.9445*(10^2/t) - 1.2222*Log[t]; 
g = b0 + b1*phip + b2*phip^2; 
dg = b1 + 2*b2*phip; 
ddg = 2*b2; 
dddg = 0; 

S1 = 1/(phip*rp*vpw) + 1/(1 - phip) - (2*g - 2*(1 - 2*phip)*dg - 
     phip*(1 - phip)*ddg); 
S2 = -(1/(phip^2*rp*vpw)) + 1/(1 - phip)^2 - 
    (6*dg - 3*(1 - 2*phip)*ddg - phip*(1 - phip)*dddg); 

   p = Apply[{phip, t} /. FindRoot[{S1 == 0, S2 == 0}, {phip, #1}, {t, #2}] & , 
   {{-2, 400}, {-0.05, 20}, {3, 800}}, {1}]

{{-1.28734, 414.625}, {-0.0436006, 23.3922}, {2.6102, 820.852}}

ContourPlot[{S1 == 0, S2 == 0}, {phip, -2, 3}, {t, 0, 1000}, 
  PlotTheme -> "Detailed", PlotPoints -> 100, 
  Epilog -> {Red, PointSize[Large], Point[p]}]

Contourplot3

   ContourPlot[{S1 == 0, S2 == 0}, {phip, -0.2, 0.3}, {t, 0, 100}, 
  PlotTheme -> "Detailed", PlotPoints -> 100, 
  Epilog -> {Red, PointSize[Large], Point[p]}]

Contourplot4

A question, how can I insert plots here?

Improve this answer
$\endgroup$
1
  • 1
    $\begingroup$ Use "Save Selection As" in Mathematica, save the plot as a gif or other format, use the image selector to place it in your post (sixth icon from the left). $\endgroup$
    – Daniel Lichtblau
    Commented Jul 28, 2015 at 17:49

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.