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

	% 2*pi/365 is radial frequency [1/year]
	% p(1) - year long amplitude
	% p(2) - shift along the year 
	% p(3) - temperature offset or median yearly temperature
	% note that if p(2)=0 temerature will be minimal i.e. it is winter
	t=-p(1)*cos(2*pi/365*d + p(2)) + p(3);
end



