A Chaotic Circuit

Analysis of a simple electric circuit that displays chaotic behavior

eqns=
	{q'[t]==i[t],
		(Subscript[L, 1]) i'[t]==-(Subscript[R, 1]) i[t] - 
                ((((Subscript[C, 2]) - (Subscript[C, 1]))/(2 (Subscript[C, 1]) (Subscript[C, 2]))) Abs[q[t]] + 
                    (((Subscript[C, 1]) + (Subscript[C, 2]))/(2 (Subscript[C, 1]) (Subscript[C, 2]))) q[t] + (Subscript[E, 0])) + 
                𝒱  Sin[2 π 𝒻 t],
		q[0] == 0,
		i[0] == 0};

eqns = 
      eqns /. {(Subscript[C, 1]) -> 0.1 ×(10^(-6)), 
          (Subscript[C, 2]) -> 400. ×(10^(-12)), (Subscript[R, 1]) -> 60., 
          (Subscript[L, 1]) -> 100.×(10^(-6)), 𝒻 -> 700000., 
          (Subscript[E, 0]) -> 0.1};

sol[v_] := 
      NDSolve[eqns /. 𝒱 -> v, {q, i}, {t, 0, 
          1.4×(10^(-4))}, MaxSteps -> (10^4)];

Do[
      Plot[Evaluate[q[t] /. sol[𝒱]], {t,  
          1.2×(10^(-4)), 1.4×(10^(-4))}, 
        PlotLabel -> 𝒱], {𝒱, 0.1, 0.3, 
        0.1}];

Do[
      ParametricPlot[
        Evaluate[{q[t], i[t]} /. sol[𝒱]], {t,  
          1.2×(10^(-4)), 1.4×(10^(-4))}, 
        PlotLabel -> 𝒱], {𝒱, 0.1, 0.3, 
        0.1}];

ListPlot[
      Flatten[Table[
          With[{s = sol[𝒱]}, 
            Table[{𝒱, First[q[t] /. s]}, {t, (10^(-4)), 
                1.4×(10^(-4)), 1/700000.}]], {𝒱, 0, 
            0.3, 0.01}], 1], Frame -> True];