function out = ptheta(theta,mu,sigma)
% Function to give the density value for a given angle theta in degrees 
% in a Gaussian distribution given mean mu and standard 
% deviation sigma in degrees.
% Tim Brown 9/6/00.

th = theta .* pi ./ 180;
m = mu .* pi ./ 180;
si = sigma .* pi / 180;
exponent = ((th - ((pi/2) - m)).^2) * -1 / (2 * si.^2);
out = exp(exponent);