Hello, I am having a problem to plot a system of two differential equations while one initial value is assigned values from a Table. This is my code
g = 5/3;
c4 = 6;
smax = 10;
Rinit = 1.05;
Linit = 0.5;
t = Flatten@Table[Linit + i, {i, 0, 0.3, 0.1}];
Sol1 = NDSolveValue[{D1[s] ==
2 (-1 + R[s]) R[
s] (-(-1 + g) L[s] +
c4 ((-1 + R[s]) R[s]^3)^((1 - g)/2) (-2 + g) g L[s]^(2 - g) +
c4 ((-1 + R[s]) R[s]^3)^((1 - g)/2) (-1 + g) g L[s]^-g),
N1[s] == -(1 + L[s]^2) (1 - g +
c4 (-1 + R[s])^(1/2 - g/2)
R[s]^(-(3/2) (-1 + g)) (2 + 4 R[s] (-1 + g) - 3 g) g L[s]^(
1 - g)), R'[s] == D1[s]/(\[Sqrt](D1[s]^2 + N1[s]^2)),
Derivative[1][L][s] == N1[s]/(\[Sqrt](D1[s]^2 + N1[s]^2)),
R[0] == Rinit, L[0] == #,
WhenEvent[
D1[s] == 0, {smax = s, Print[N1[s]/D1[s]], "StopIntegration"},
"DetectionMethod" -> "DerivativeSign"]}, {R, L}, {s, 0,
smax}] & /@ t
(*s1=ParametricPlot[Evaluate@Through[Sol1[{R[s],L[s]}]],{s,0,smax},\
PlotStyle\[Rule]Red,AxesLabel\[Rule]{R,L},AspectRatio\[Rule]Full,\
PlotRange\[Rule]All]*)
(*ParametricPlot[{Rsol[s],Lsol[s]},{s,0,smax},PlotStyle\[Rule]Red,\
AxesLabel\[Rule]{R,L},AspectRatio\[Rule]Full,PlotRange\[Rule]All]*)
I think there is a mistake on the ParametricPlot, but I am not sure. I have included the file. Thank you.
Attachments: