SetOptions[Plot,PlotStyle->{Black,AbsoluteThickness[2]},BaseStyle->{16}];

Input appropriate parameters for the LRC circuit lab below:

width:=300

R=800;
L=.540;
C0=10^(-8);
ω0=1/Sqrt[L*C0];
β=R/2/L;
Y[ω_]:=(ω/L)/Sqrt[(4*β^2*ω^2+(ω0^2-ω^2)^2)]
ϕ[ω_]:=π/2+ArcTan[(ω0^2-ω^2),-R*ω/L]
fs:=16
Plot[Y[ω], {ω,0, 3 ω0},
PlotRange->All,
BaseStyle->{FontSize->fs,FontFamily->"Times"},
AxesLabel->{"ω (rad/s)","mod Y"},
PlotStyle->{Purple,AbsoluteThickness[3]},
Ticks->{{{ω0,"(Subscript[ω, 0])"},{2ω0,"2(Subscript[ω, 0])"},{3ω0,"3(Subscript[ω, 0])"}},None},
ImageSize->width]
Plot[ϕ[ω], {ω,0, 3 ω0},
PlotRange->{-π/2,π/2},
BaseStyle->{FontSize->fs,FontFamily->"Times"},
AxesLabel->{"ω (rad/s)","phase(Y)"},
PlotStyle->{Magenta,AbsoluteThickness[3]},
Ticks->{{{ω0,"(Subscript[ω, 0])"},{2ω0,"2(Subscript[ω, 0])"},{3ω0,"3(Subscript[ω, 0])"}},None},
ImageSize->width]

Here apply sinusoidal voltages one at a time.  n determines freq in units of the fundamental.  n1=1 => fund freq = re freq.  n2=2 => fund freq = 2* res freq.  Width detrmines size of image.

width:=300

n1=.7;
T=2 π/(ω0);
V1[t_]:=Exp[I*n1*2 π t/T];
curr1[t_]:=R*Y[n1 ω0]*Exp[I*ϕ[n1 ω0]]*V1[t];
Plot[{Re[V1[t]],1*Re[curr1[t]]} ,{t,0, 3T},PlotRange->{-2,2},AxesLabel->{"t (s)",None},Ticks->{{{T,"(Subscript[T, 0])"},{2T,"2(Subscript[T, 0])"},{3T,"3(Subscript[T, 0])"}},{-1,1}}, PlotStyle->{{Red,AbsoluteThickness[2]},{Blue,AbsoluteThickness[2]}},
ImageSize->width]

n2=1.4;
T=2 π/(ω0);
V2[t_]:=-1*Exp[I*n2*2 π t/T];
curr2[t_]:=1*R*Y[n2 ω0]*Exp[I*ϕ[n2 ω0]]*V2[t];
Plot[{Re[V2[t]],Re[curr2[t]]} ,{t,0, 3T},
PlotRange->{-2,2},AxesLabel->{"t (s)",None},
Ticks->{{{T,"(Subscript[T, 0])"},{2T,"2(Subscript[T, 0])"},{3T,"3(Subscript[T, 0])"}},{-1,1}}, PlotStyle->{{Red,AbsoluteThickness[2]},{Blue,AbsoluteThickness[2]}},
ImageSize->width]