PlanetaryAtmospheres: phi.m

function y=phi(x,lambda)
  y=exp(-lambda/x*(x-1)-mypsi(x,lambda)/2);
end