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:

Code1 alt.nb