Question: explicate plotting


restart:with(MultiSeries):with(plots):

parvalues:={C=2.02,H_liq=8.74,R=22.11,R1=0.0006};

(1)

k:=27:

G:=R3*k*(1-R1)/(C*Pec_i);

(2)

F1:=1/(R*S+1-S);

(3)

F:=F1*(1-R*G*S);

(4)

F2:=Pec_i*F;

(5)

A1:=Theta0/C=(exp(F2*(1-S))-1)*(H_liq+1/(1-exp(-C/k*F2*(S))));

(6)

p1:=plots:-implicitplot([subs(parvalues,Theta0=4,R3=0,Pec_i=-20,A1),subs(parvalues,Theta0=4,R3=-1,Pec_i=-20,A1),subs(parvalues,Theta0=4,R3=-3,Pec_i=-20,A1),subs(parvalues,R3=-4,Pec_i=-20,Theta0=4,A1)],k=0..20,S=0..1,color=[red,blue,green,black],numpoints=10000):plots:-display(p1);

 

I want to plot F2 vs k using the same data used for the above plot, but F2 is dependent on S (range of S is (0, 1)) and, S has to be found as apart of the solution (equation A1). Any idea

 


Download section_.mws

Please Wait...