activityFunction(timezoneOffset, UTC)

function [activity] = activityFunction(timezoneOffset, UTC) 

localTime=mod(UTC+timezoneOffset+24,24);

if localTime < 6,
    activity=0.1;
else
    activity=1-0.9*(cos(((localTime-6)*pi)/18))^10;
end