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)]))$
fr1/(%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'"])$

The attached file shows two versions, one merely literal translation and one more "paper and pencil" like version.
Both versions exploit the still experimental Maxima Draw2D region.

---

Section Graphics Maxima Draw region.sm (112 KiB) downloaded 52 time(s).

Excellent on window with maxima plugin.
As I work on linux I am trying to convert them
I can't make resonance plot (f|Vout|)
include example

---

rlcsmath.sm (9 KiB) downloaded 48 time(s).