function t=temp_vs_day_time(d,p)
	% d - day time
	% t - temperature outside
	% p input parameters
	% special care is taken to make it arrys data compatible

	% 2*pi/1is radial frequency [1/day]
	% p(1) - daily amplitude
	% p(2) - shift with respect to midnight (minimum temperature lag)
	% note that if p(2)=0 temerature will be minimal i.e. it is midnight
	t=-p(1)*cos(2*pi/1*d - p(2));
end
