I have a set of data points giving me the values for the second virial coefficient, for various values of T, of the virial expansion which is an equation that corrects the ideal gas law for empirical deviations: I'm trying to do a least squares fit to determine how well the van der Waals equation predicts using MATLAB.
The equation you derive for B using van der Waals' equation ends up being:
where a and b are unknown constants I need to determine and R is just the ideal gas constant.
Now my data points for true values are here:
B = [-160 -35 -4.2 9 16.9 21.3]*10^-6 ;
T = [100 200 300 400 500 600]; (kelvin)
How can I write some code which will use the least squares method to generate estimates for a and b using the given data points?
My problem is that this is essentially a equation which I'm not sure how to represent in MATLAB.
I can only find options for quadratic, cubic, polynomial fit etc.