Add a numpy version of the beta function

This takes a numpy array of temperatures - easer
for plotting and other manipulation.
Automated script to do the mass substitution
Now calculate beta and fit to function
Add tool to plot beta