function result = windChill(temperatureF,velocityMPH)
%windchill  compute the windchill with NWS formula (www.weather.gov)
% inputs:  temperatureF (temperature in Fahrenheit)
%          velocityMPH (wind speed in miles per hour)
% output:
%          windchill factor
%
% Examples:
%         >>  windChill(20,30)
%         ans = 1
%         >>  windChill(30,20)
%         ans = 17
%
% P. Conrad for CISC106, 10/03/2007

   result = 35.74 + (0.6215 * temperatureF) - 35.75*(velocityMPH^0.16) ...
	    + 0.4275*temperatureF*(velocityMPH^0.16);

  return;
end % function windChill