9

program
.model raf res r=1 dev/gauss=15%
.model caf cap c=1 dev/uniform=10%
**         we+ we- wy
.subckt wo 1 2 3
r1 1 4 .5meg
rw 4 2 RAF .5meg
E1 5 0 4 2 200k
r3 5 3 1
c1 3 0 caf 0.032
.ends
x1 1 0 3 wo
v1 1 0 ac 1u dc 1u
.dc v1 -1u 1u .1u
.ac DEC 100 1 10meg
.mc 600 DC V(3) ymax list output all
.probe
.end
**cutoff_lowpass_3db(V(2)) w performance analysis
**SQRT(R(V(1)/I(r1))^2 +IMG(V(1)/I(r1))^2)