oxygen_fit_poly6

 Function to calculate the dissolved oxygen

0001 function yhat=oxygen_fit_poly6(beta,x)
0002 % Function to calculate the dissolved oxygen
0003 %
0004 b1=beta(1);
0005 b2=beta(2);
0006 b3=beta(3);
0007 b4=beta(4);
0008 b5=beta(5);
0009 b6=beta(6);
0010 b7=beta(7);
0011 b8=beta(8);
0012 b9=beta(9);
0013 b10=beta(10);
0014 b11=beta(11);
0015 x1=x(:,1);
0016 x2=x(:,2);
0017 x3=x(:,3);
0018 x4=x(:,4);
0019 x5=x(:,5);
0020 yhat=((b1*(x1+b2))+(b3*exp(-0.03*x3))+(b6*x5+(b7*x5.^2)+(b8*x5.^3)+(b9*x5.^4)+(b10*x5.^5)+(b11*x5.^6))).*exp(b4.*x3).*x2.*exp(b5.*x4);