oxygen_fit_exp

 Function to calculate the dissolved oxygen

0001 function yhat=oxygen_fit_exp(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 x1=x(:,1);
0012 x2=x(:,2);
0013 x3=x(:,3);
0014 x4=x(:,4);
0015 x5=x(:,5);
0016 yhat=((b1*(x1+b2))+(b3*exp(-0.03*x3))+b6*exp(b7*x5)).*exp(b4.*x3).*x2.*exp(b5.*x4);