# set terminal pngcairo  transparent enhanced font "arial,10" fontscale 1.0 size 600, 400 
# set output 'multiaxis.1.png'
set bar 1.000000 front
set style circle radius graph 0.02, first 0.00000, 0.00000 
set style ellipse size graph 0.05, 0.03, first 0.00000 angle 0 units xy
set dummy jw, y
set grid nopolar
set grid xtics nomxtics noytics nomytics noztics nomztics \
 nox2tics nomx2tics y2tics nomy2tics nocbtics nomcbtics
set grid layerdefault   lt black linewidth 0.200 dashtype solid,  lt black linewidth 0.200 dashtype solid
set key inside center top vertical Right noreverse enhanced autotitle nobox
set style textbox transparent margins  1.0,  1.0 border
unset logscale
set logscale x 10
set logscale y 10
set logscale x2 10
set style data lines
set xtics border out scale 1,0.5 mirror norotate  autojustify
set ytics border out scale 1,0.5 nomirror norotate  autojustify
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 
set cbtics border out scale 1,0.5 mirror norotate  autojustify
set rtics axis out scale 1,0.5 nomirror norotate  autojustify
set paxis 1 tics border out scale 1,0.5 nomirror norotate  autojustify
set paxis 1 tics  rangelimit autofreq 
set paxis 2 tics border out scale 1,0.5 nomirror norotate  autojustify
set paxis 2 tics  rangelimit autofreq 
set paxis 3 tics border out scale 1,0.5 nomirror norotate  autojustify
set paxis 3 tics  rangelimit autofreq 
set paxis 4 tics border out scale 1,0.5 nomirror norotate  autojustify
set paxis 4 tics  rangelimit autofreq 
set paxis 5 tics border out scale 1,0.5 nomirror norotate  autojustify
set paxis 5 tics  rangelimit autofreq 
set paxis 6 tics border out scale 1,0.5 nomirror norotate  autojustify
set paxis 6 tics  rangelimit autofreq 
set paxis 7 tics border out scale 1,0.5 nomirror norotate  autojustify
set paxis 7 tics  rangelimit autofreq 
set title "Transistor Amplitude and Phase Frequency Response" 
set xlabel "jw (radians)" 
set xrange [ 1.10000 : 90000.0 ] noreverse nowriteback
set x2range [ 1.10000 : 90000.0 ] noreverse nowriteback
set ylabel "magnitude of A(jw)" 
set y2label "Phase of A(jw) (degrees)" 
set paxis 1 range [ * : * ] noreverse nowriteback
set paxis 2 range [ * : * ] noreverse nowriteback
set paxis 3 range [ * : * ] noreverse nowriteback
set paxis 4 range [ * : * ] noreverse nowriteback
set paxis 5 range [ * : * ] noreverse nowriteback
set paxis 6 range [ * : * ] noreverse nowriteback
set paxis 7 range [ * : * ] noreverse nowriteback
set colorbox vertical origin screen 0.9, 0.2, 0 size screen 0.05, 0.6, 0 front  noinvert bdefault
A(jw) = ({0,1}*jw/({0,1}*jw+p1)) * (1/(1+{0,1}*jw/p2))
p1 = 10
p2 = 10000
GPFUN_A = "A(jw) = ({0,1}*jw/({0,1}*jw+p1)) * (1/(1+{0,1}*jw/p2))"
plot abs(A(jw)) axes x1y1, 180./pi*arg(A(jw)) axes x2y2