Re:1;mu:0.1,numer;
epsilon:0.001,numer;
z:sqrt((mu/epsilon)),numer;
%omega:float(sqrt(1/(mu*epsilon)));
R:z$freq:[]$H:[]$
for f:2 step 1 thru 100
do (w:2*%pi*f,Zs:R+%i*w*mu-%i/(w*epsilon),
Vout:R/Zs,freq:append(freq,[f]),
H:append(H,[abs(Vout)]))$
fr

1/(%pi*2*sqrt(mu*epsilon))),numer;
wxplot2d([discrete,float(freq),float(H)],
[logx],[xlabel,"Fréquence (MHz)"],
[ylabel,"Impédance"], [style,[lines,1,11]],
[gnuplot_preamble,
"unset key; set title 'Fréquence de résonance MHz'"])$