# set terminal pngcairo  transparent enhanced font "arial,10" fontscale 1.0 size 600, 400 
# set output 'electron.2.png'
set dummy Vd, 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 key at 0.200000, 0.00450000 right top vertical Right noreverse enhanced autotitle box lt black linewidth 1.000 dashtype solid
unset key
set label 1 "-0.5 Vp" at 4.10000, 0.625000, 0.00000 left norotate back nopoint
set label 2 "-0.25 Vp" at 4.10000, 1.40000, 0.00000 left norotate back nopoint
set label 3 "0" at 4.10000, 2.50000, 0.00000 left norotate back nopoint
set label 4 "Vg = 0.5 Vp" at 4.10000, 3.90000, 0.00000 left norotate back nopoint
set samples 160, 160
set style data lines
set title "JFET Mutual Characteristic" 
set xlabel "Drain voltage Vd (V)" 
set xrange [ 0.00000 : 4.00000 ] noreverse nowriteback
set x2range [ * : * ] noreverse writeback
set ylabel "Drain current Id (mA)" 
set yrange [ 0.00000 : 5.00000 ] noreverse nowriteback
set y2range [ * : * ] noreverse writeback
set zrange [ * : * ] noreverse writeback
set cbrange [ * : * ] noreverse writeback
set rrange [ * : * ] noreverse writeback
set rmargin  12
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)
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
plot Vg=0.5*Vp,Id(Vd),Vg=0.25*Vp,Id(Vd),Vg=0,Id(Vd),Vg=-0.25*Vp,Id(Vd)