oxy_fit_poly4

0001 function [new_oxy,Soc,Voffset,Boc,tcor,pcor,p1,p2,p3,p4]=...
0002     oxy_fit_poly4(volt,OXSAT,temp,press,sta,bot_oxy,Soc,Voffset,Boc,tcor,pcor,p1,p2,p3,p4);
0003 x=[volt OXSAT temp press sta];
0004 beta(1,1)=Soc;
0005 beta(1,2)=Voffset;
0006 beta(1,3)=Boc;
0007 beta(1,4)=tcor;
0008 beta(1,5)=pcor;
0009 beta(1,6)=p1;
0010 beta(1,7)=p2;
0011 beta(1,8)=p3;
0012 beta(1,9)=p4;
0013 beta0=non_lin_fit(x,bot_oxy,'oxygen_fit_poly4',beta);
0014 Soc=beta0(1);
0015 Voffset=beta0(2);
0016 Boc=beta0(3);
0017 tcor=beta0(4);
0018 pcor=beta0(5);
0019 p1=beta0(6);
0020 p2=beta0(7);
0021 p3=beta0(8);
0022 p4=beta0(9);
0023 new_oxy=((Soc.*(volt+Voffset))+Boc.*exp(-0.03.*temp)+(p1.*sta+(p2.*sta.^2)+(p3.*sta.^3)+(p4.*sta.^4))).*exp(tcor.*temp).*OXSAT.*exp(pcor.*press);
0024 return