% Got some more info via fax from Pasco guy

l=[2325.4,1970.1,1529.6,1060.0,1014.0,852.1,706.5,656.3,643.8,632.8,589.3,587.6,546.1,486.1,480.0,435.8,404.7];
n=[1.68218,1.68750,1.69378,1.70227,1.70345,1.70887,1.71661,1.72085,1.72200,1.72309,1.72803,1.72825,1.73430,1.74648,1.74805,1.76198,1.77579];

p=matmin('chisq',[1.7,0,0,0,0],[],'poly',n,ones(size(n)),l);

A=6015783.2978286
B=-13717512.2572366
C=11725406.1052594
D=-4452989.0254548
E=633985.0106661
pl=3000./(A+B*n +C*n.^2+D*n.^3+E*n.^4).^0.5

plot(l,n,'.-')
hold on
plot(pl,n,'*r-');
%l=400:20:3000;
%plot(l,cauchy(p,l),'r');
%plot(l,poly(p,l),'r');
hold off
