fx=(b(1)+b(2)*x+b(3)*x.^2+b(4)*x.*y+b(5)*x.^3)./(1+b(6)*x+b(7)*x.*y+b(8)*x.^2.*y);
SSy = 29114
b =[3.536419787 1.508853635 -0.03445725593 0.06125633835 0.01355562355 -0.0338963446 0.0007619172333 -2.277850887e-05]
RSS = 10.2281528862
MSe = 0.10654
R^2 = 0.99965
