function PH=HGfun(x,g)
% Henyey Greenstein Phase function
% This function isn't normalized to 1 only to 4pi
% x-scattering angle in deg
% g-asymetric parameter

x=x*pi/180;

PH=(1-g^2)./(1+g^2-2*g*cos(x)).^(3/2);


