function out=tempfunc(zsr,p);
Co=1.7175;
Cm=3.6799;
Zm=2023.6;
out = Co+(Cm-Co)*[(exp(zsr/p(1))-1)/(exp(Zm/p(1))-1)];
%out = Co+(Cm-p(1))*[(exp(zsr/p(2))-1)/(exp(Zm/p(2))-1)];