oxy_fit_poly5

0001 function [new_oxy,Soc,Voffset,Boc,tcor,pcor,p1,p2,p3,p4,p5]=...
0002     oxy_fit_poly5(volt,OXSAT,temp,press,sta,bot_oxy,Soc,Voffset,Boc,tcor,pcor,p1,p2,p3,p4,p5);
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 beta(1,10)=p5;
0014 beta0=non_lin_fit(x,bot_oxy,'oxygen_fit_poly5',beta);
0015 Soc=beta0(1);
0016 Voffset=beta0(2);
0017 Boc=beta0(3);
0018 tcor=beta0(4);
0019 pcor=beta0(5);
0020 p1=beta0(6);
0021 p2=beta0(7);
0022 p3=beta0(8);
0023 p4=beta0(9);
0024 p5=beta0(10);
0025 new_oxy=((Soc.*(volt+Voffset))+Boc.*exp(-0.03.*temp)+(p1.*sta+(p2.*sta.^2)+(p3.*sta.^3)+(p4.*sta.^4)+(p5.*sta.^5))).*exp(tcor.*temp).*OXSAT.*exp(pcor.*press);
0026 return