# set terminal pngcairo  transparent enhanced font "arial,10" fontscale 1.0 size 600, 400 
# set output 'electron.3.png'
set dummy jw, y
set format x "%g" 
set format y "%g" 
set format x2 "%g" 
set format y2 "%g" 
set format z "%g" 
set format cb "%g" 
set format r "%g" 
set grid nopolar
set grid xtics nomxtics noytics nomytics noztics nomztics nortics nomrtics \
 nox2tics nomx2tics y2tics nomy2tics nocbtics nomcbtics
set grid layerdefault   lt 0 linecolor 0 linewidth 0.500,  lt 0 linecolor 0 linewidth 0.500
set key fixed center bottom vertical Right noreverse enhanced autotitle box lt black linewidth 1.000 dashtype solid
set samples 160, 160
set style data lines
set xtics border out scale 1,0.5 mirror norotate  autojustify
set xtics  norangelimit logscale autofreq 
set ytics border out scale 1,0.5 nomirror norotate  autojustify
set ytics  norangelimit logscale autofreq  textcolor rgb "dark-violet" 
set ztics border out scale 1,0.5 nomirror norotate  autojustify
set y2tics border out scale 1,0.5 nomirror norotate  autojustify
set y2tics  norangelimit autofreq  textcolor rgb "#56b4e9" 
set cbtics border out scale 1,0.5 mirror norotate  autojustify
set rtics axis out scale 1,0.5 nomirror norotate  autojustify
set title "Amplitude and Phase Frequency Response" 
set xlabel "jw (radians)" 
set xrange [ 1.10000 : 90000.0 ] noreverse nowriteback
set x2range [ * : * ] noreverse writeback
set ylabel "magnitude of A(jw)" 
set y2label "Phase of A(jw) (degrees)" 
set yrange [ * : * ] noreverse writeback
set y2range [ * : * ] noreverse writeback
set zrange [ * : * ] noreverse writeback
set cbrange [ * : * ] noreverse writeback
set rrange [ * : * ] noreverse writeback
set logscale y 10
set logscale x 10
set logscale x2 10
set colorbox vertical origin screen 0.9, 0.2 size screen 0.05, 0.6 front  noinvert bdefault
Ie(Vbe)=Ies*exp(Vbe/kT_q)
Ic(Vbe)=alpha*Ie(Vbe)+Ico
Ida(Vd)=Ido*(1-Vg/Vp)**2
Idb(Vd)=Ido*(2*Vd*(Vg-Vp)-Vd*Vd)/(Vp*Vp)
Id(Vd)= (Vd>Vg-Vp) ? Ida(Vd) : Idb(Vd)
A(jw) = ({0,1}*jw/({0,1}*jw+p1)) * (1/(1+{0,1}*jw/p2))
NO_ANIMATION = 1
Ies = 4e-14
kT_q = 0.025
alpha = 0.99
Ico = 1e-09
Ido = 2.5
Vg = 0.3125
Vp = -1.25
p1 = 10
p2 = 10000
plot abs(A(jw)) lt 1, 180/pi*arg(A(jw)) axes x1y2 lt 3