function res=sundist(day);
%Earth-sun distance ratio calculation (d0/d)^2
% input : 
%       day - Julian day (DAY of Year)

res=1+(0.034.*cos(2*pi.*day/365));