-3
$\begingroup$

This question already has an answer here:

I have combined three plots with "Show[p1,p2,p3]. Can I extract the data points from my plot and want to plot with origin?? enter image description here`

F = 0.025;
kv = 0.0043;
kw = 0.046;
A = (F + kv);
B = (F + kw);
Du = 290;
Dv = 100.443;
Dw = 7.3;
G = (F - Sqrt[F^2 - 4*F*(A^2 + B^2)])^2;
j11 = (-((4*F^2*(A^2 + B^2))*G^(-1)) - F);
j12 = -(2*A);
j13 = -(2*B);
j21 = (4*A^2*F^2)/(G);
j22 = (A);
j31 = (4*B^2*F^2)/G;
j33 = B;

C1 = (k^2 (Du + Dv + Dw) - j11 - j22 - j33);

C2 = ((-j11*(Dv + Dw) - j22*(Du + Dw) - j33*(Du + Dv))*
     k^2 + (Du*Dv + Du*Dw + Dv*Dw)*k^4 - j12*j21 - j13*j31 + j11*j22 +
     j11*j33 + j22*j33);

C3 = ((Du*Dv*Dw)*k^6 + (-j11*Dv*Dw - j22*Du*Dw - j33*Du*Dv)*
     k^4 + (Dw*j11*j22 + Dv*j11*j33 + Du*j22*j33 - Dw*j12*j21 - 
       Dv*j13*j31)*k^2 - j11*j22*j33 + j12*j21*j33 + j13*j31*j22);
p1 = ContourPlot[
  x^3 + (C1) *x^2 + (C2)*x + C3 == 0, {k, 0.005, 0.020}, {x, -0.01, 
   0.008}, ContourStyle -> Magenta]

` Then repeated the same code with Different Dv=149 and Dv=174, combine the three plots with Show[p1,p2, p3]

$\endgroup$

marked as duplicate by xzczd, Henrik Schumacher, gwr, Anton Antonov, b3m2a1 Jan 17 at 21:17

This question has been asked before and already has an answer. If those answers do not fully address your question, please ask a new question.

  • $\begingroup$ Strongly related: mathematica.stackexchange.com/q/14222/1871 mathematica.stackexchange.com/q/31858/1871 $\endgroup$ – xzczd Jan 14 at 7:39
  • $\begingroup$ This question might be a duplicate, but also might not be, depending how the mentioned combined plot looks like. Can you provide an image of your combined plot? $\endgroup$ – Anton Antonov Jan 14 at 14:39
  • $\begingroup$ Anton AntonovThank you so much for your reply I have added the image file here. $\endgroup$ – Albert Jan 17 at 1:33
  • 3
    $\begingroup$ All this fuss and to-do over the lack of an uppercase O on origin in the original post ! $\endgroup$ – High Performance Mark Jan 17 at 15:11
  • 1
    $\begingroup$ @HighPerformanceMark Yeah, a complete waste of time... $\endgroup$ – Anton Antonov Jan 17 at 15:16
9
$\begingroup$

@AntonAntonov I want to plot with origin. – Albert

I think OP wants this (it is not entirely clear):

F = 0.025;
kv = 0.0043;
kw = 0.046;
A = (F + kv);
B = (F + kw);
Du = 290;
(*Dv=100.443;*)
Dw = 7.3;
G = (F - Sqrt[F^2 - 4*F*(A^2 + B^2)])^2;
j11 = (-((4*F^2*(A^2 + B^2))*G^(-1)) - F);
j12 = -(2*A);
j13 = -(2*B);
j21 = (4*A^2*F^2)/(G);
j22 = (A);
j31 = (4*B^2*F^2)/G;
j33 = B;

{kmin, kmax} = {0.005, 0.020};

grs =
  MapThread[
   Function[{Dv, color},
    C1 = (k^2 (Du + Dv + Dw) - j11 - j22 - j33);
    C2 = ((-j11*(Dv + Dw) - j22*(Du + Dw) - j33*(Du + Dv))*
        k^2 + (Du*Dv + Du*Dw + Dv*Dw)*k^4 - j12*j21 - j13*j31 + 
       j11*j22 + j11*j33 + j22*j33);
    C3 = ((Du*Dv*Dw)*k^6 + (-j11*Dv*Dw - j22*Du*Dw - j33*Du*Dv)*
        k^4 + (Dw*j11*j22 + Dv*j11*j33 + Du*j22*j33 - Dw*j12*j21 - 
          Dv*j13*j31)*k^2 - j11*j22*j33 + j12*j21*j33 + j13*j31*j22);
    ContourPlot[
     x^3 + (C1)*x^2 + (C2)*x + C3 == 0, {k, kmin, kmax}, {x, -0.01, 
      0.008}, ContourStyle -> color, PlotRange -> All]
    ], {{100.443, 149, 174}, {Red, Blue, Green}}];

Show[grs, PlotRange -> {{-kmin, kmax}, All}, GridLines -> {{0}, {0}}, 
 Epilog -> {Red, PointSize[0.03], Point[{0, 0}], Black, 
   Text[Style["Origin", Bold], {0, 0}, {-1, -1}]}]

enter image description here

$\endgroup$
  • $\begingroup$ +1 for banter related arguments. $\endgroup$ – QuantumPenguin Jan 17 at 15:32
  • $\begingroup$ @QuantumPenguin Thanks, I spent too much time coding in order to get OP to say from "origin" I mean to Plot in "origin" software. $\endgroup$ – Anton Antonov Jan 17 at 15:35
  • 2
    $\begingroup$ Totally worth it, this really tickled me. $\endgroup$ – QuantumPenguin Jan 17 at 15:40
  • 1
    $\begingroup$ Totally worth an upvote. $\endgroup$ – halirutan Jan 17 at 19:05

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