function incidencecal
%msgbox('solar software','calculating')
%input w n gama(surface azimuth angle slop of plate,phe
h1=findobj(0,'tag','now');
m=findobj(h1,'tag','EditText2');
phe=get(m,'string');
phe=str2num(phe)
h2=findobj(0,'tag','incidencef');
t=findobj(h2,'tag','EditText1hour');
wk=get(t,'string');
w=str2num(wk)
t=findobj(h2,'tag','EditText2day');
nk=get(t,'string');
n=str2num(nk)
t=findobj(h2,'tag','EditText3azimuth');
gamak=get(t,'string');
gama=str2num(gamak)
t=findobj(h2,'tag','EditText4slop');
bk=get(t,'string');
b=str2num(bk)
if strcmp(wk,'')|strcmp(nk,'')|strcmp(gamak,'')|strcmp(bk,'')
msgbox('you should enter all the required number in balanks')
else
del=23.45*sin((360*(n+284)/365)*3.14/180);
delrad=del*pi/180;
wdegree=(w-12)*15;
wrad=wdegree*(pi/180);
pherad=phe*pi/180;
brad=b*pi/180;
gamarad=gama*pi/180;
coosincidenceangle=sin(delrad)*sin(pherad)*cos(b)-sin(delrad)*cos(pherad)*sin(brad)*cos(gamarad)+cos(delrad)*cos(pherad)*cos(brad)*cos(wrad)+cos(delrad)*sin(pherad)*sin(brad)*cos(gamarad)*cos(wrad)+cos(delrad)*sin(brad)*sin(gamarad)*sin(wrad)
incidence=(180/pi)*acos(coosincidenceangle);
msgbox(num2str(incidence),'incidence angle')
end