Also, you might find ParametricNDSolveValue useful:

solPar2 = 
  ParametricNDSolveValue[Flatten@{Eqns1, init1}, 
   vars, {t, 0, 0.0001}, {\[ScriptD]}];

Plot[solPar2[.5] // Evaluate, {t, 0, .0001}, PlotRange -> All, 
 PlotLegends -> vars]

Hello, Jesus.

I don't know how the parametric feature works. And I also find it puzzling. I find it strange that we are not required to provide a range for the parameters. When we don't, I wonder what the implied range is, since it seems unreasonable to expect a parameter in such cases to be valid over the extent of R or C. I do see that we can provide a range in the form a {p, pmin, pmax}. That may be a safer alternative than just p.

It would be good for someone from Wolfram to comment! Maybe we can wake up the Moderators.

Kind regards, David

Hi Jesus,

In your note book you try to plot with

Plot[iL[0.5][t] /. solPar, {t, 0, 0.0001}, PlotRange -> All]

But really iL[t] is the returned parametric function, so you should plot

Plot[iL[t][0.5] /. solPar, {t, 0, 0.0001}, PlotRange -> All]

Best regards,

David