Dear Maple help me to to plot the graph please see and rectify. thanks in advance
i am attaching the codes
inf:=5:
pdes:= R(X,R,t)*diff(U(X,R,t),X)+U(X,R,t)*diff(R(X,R,t),X)+R(X,R,t)*diff(V(X,R,t),R)+V(X,R,t),
diff(U(X,R,t),t)+U(X,R,t)*diff(U(X,R,t),X)+V(X,R,t)*diff(U(X,R,t),R)=Gr*T(X,R,t)+Gc*C(X,R,t)+(1/R(X,R,t))*diff(R*diff(U(X,R,t),R),R),
diff(T(X,R,t),t)+U(X,R,t)*diff(T(X,R,t),X)+V(X,R,t)*diff(T(X,R,t),R)+(1/(Pr*R(X,R,t)))*diff(R*diff(T(X,R,t),R),R),diff(C(X,R,t),t)+U* diff(C(X,R,t),X)+V(X,R,t)*diff(C(X,R,t),R)+(1/(Sc*R(X,R,t)))*diff(R*diff(C(X,R,t),R),R):
conds:= U(X,R,0)=0, V(X,R,0)=0, T(X,R,0)=0, C(X,R,0)=0,
U(X,1,t)=1, V(X,1,t)=0, T(X,1,t)=1, C(X,1,t)=1, U(0,R,t)=0, T(0,R,t)=0, C(0,R,t)=0,
U(X,int,t)=0,T(X,int,t)=0,C(X,int,t)=0:
pars:= { Gr=5, Gc=10,Sc=2.0}
pars := {Gc = 10, Gr = 5, Sc = 2.0}
PrVals:=[0.71, 1.00, 1.25, 2.00]:
colors:=[red, green, blue, black]:
for j from 1 to numelems(PrVals) do
pars1:=`union`( pars, {Pr=PrVals[j]}):
pdSol:= pdsolve( eval([pdes], pars1),
eval([conds], pars1),
numeric
);
plt[j]:=pdSol:-plot( U(X,R,t),X=1, t=2, R=0..inf, numpoints=200, color=colors[j]);
od:
plots:-display( [seq(plt[j], j=1..numelems(PrVals))]);