oxy_fit_poly3

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