function out=oxyfunc(zsr,p);
Co=156.5852;
Cm=186.3462;
Zm=2023.6;
%Zstar=497.4128;
%Zstar=589.6686;
Zstar=543.5407;
%out=Co+(Cm-p(1))*[(exp(zsr/Zstar)-1)/(exp(Zm/Zstar)-1)]+p(2)*(zsr-Zm*[(exp(zsr
% /Zstar)-1)/(exp(Zm/Zstar)-1)]);

out=Co+(Cm-Co)*[(exp(zsr/Zstar)-1)/(exp(Zm/Zstar)-1)]+p(1)*(zsr-Zm*[(exp(zsr/Zstar)-1)/(exp(Zm/Zstar)-1)]);