*** Begin of history file: Fri Apr 11 07:50:30 2003
*** Begin of history file: Thu Jan 30 12:46:15 2003
* Ceate a vector of length 5 for pressure values
vec/cre P(5) r 1.036 .9104 .3010 .1919 .0953
* Use this vector to make a vector of errors
sigma ep=p*.05
* Create a vector of temperatures
vec/cre T(5) r .01 -38.8344 -189.3442 -218.7916 -248.5939
* fit the vector with a polynomial of degree 1 (a straight line)
vec/fit T P ep p1
* create a vector of temperature errors (need for fancy plotting)
sigma et=t*0
* put fancy error bars on our plot
hplot/err t p et ep 5
* put fancy labels on axes
atitle 'Temperature in Kelvin' 'Pressure in Atomospheres'
*Want to print this out, save it to a file this way
*lets make a fancy histogram first
hi/cre/1d 100 'Absolute Zero Determination' 10 -300 20
hi/set/max 100 1.1
hi/set/min 100 -0.1
fort/fil 66 myplot.ps
meta 66 -111
hi/pl 100
fun/pl 1.01315+0.00370454*x -300 10 s
atitle 'Temperature in Kelvin' 'Pressure in Atomospheres'
hplot/err t p et ep 5
text -275 .9 'Fit to,  Pressure = aTemperature +b   gives,' 0.3
text -275 .8 'a=3.70E-3 +/- 1.3E-4' 0.3
text -275 .75 'b=1.013 +/- 0.031' 0.3
close 66

*** End   of history file: Fri Apr 11 07:50:54 2003