! Physics within FMS since 1999

year=[1999:2012]
FMStotals={1073;1207;1157;1157;1401;1397;1385;1377;1379;1279;1405;1436;1597;1733}
PHYStotals={18;24;30;26;41;51;47;48;48;53;70;72;76;73}
frac=PHYStotals/FMStotals*100

orientation portrait
defaults
clear
wind 3

set %yuaxis 95
!set histyp 2
set lintyp 9
set font triumf.2
set %ylabsz 2.5
set %xlabsz 2.5
scales 0 0 0 0 500 5
set nsyinc 2
set nsxinc 2

label\y `Total Head Count, 1999=100%'
label\x `Academic year beginning September of ...'
legend on
legend frame 25 65 50 95
legend frame off
legend nsymb 3

col 1 1
set pchar -11 
set %charsz 1.5
FMStotals=FMStotals/FMStotals[1]*100
graph `FMS' year,FMStotals
scalar\vary a
fit FMStotals=100+a*(year-1999)
fit\update trend
set pchar 0
graph\noax rchar(a,`G7.2')//`%/year' year,trend

col 2 2 
set pchar -9
set %charsz 2
Phystotals=Phystotals/PHYStotals[1]*100
graph\noax `Physics' year,Phystotals
fit Phystotals=100+a*(year-1999)
fit\update trend
set pchar 0
graph\noax rchar(a,`G8.3')//`%/year' year,trend
replot

terminal
wind 4

label\y `Physics Head Counts, % of FMS'
!clear
col 1 1
set pchar -9
scales 0 0 0 0 0 0
graph\ax year,frac
col 2 2
graph\noax `Physics in FMS' year,frac
scalar\vary a,b
fit frac=b+a*(year-1999)
fit\update trend
set pchar 0
set lintyp 1
graph\noax rchar(a,`G8.2')//`%/year' year,trend