clc;clear;n=[17 47 75 105 135 165 198 228 258 288 318 344];% day number for monthly averagephi=12.9833;%latitudeABA=[4.93 5.89 6.64 6.72 6.12 5.24 4.73 4.80 5.01 4.42 4.06 4.24];% Hg in kWh/m2-day Hg=ABA*3600; % in kJ/m2-dayB=23;%tilt anglep=0.2; %reflectivity of groundg=0;%surface azimuthal no.n1=1; %refractive index of airn2=1.42 ;%refractive index of materialm=1; %no. of glass coverK=15; %extinction cofficientd=0.004;%thickness of glass coveraf=0.94; % plate absortivityfor i=1:12 D(i)=23.45*(sind((360/365)*(n(i)+284)));%declination anglews(i)= acosd(-tand(D(i))*tand(12.9833));Isc=1.367; %in W/m2 Ho(i)=(24/pi)*Isc*3600*[1+(0.033*(cosd(n(i)*360/365)))]*[(ws(i)*pi/180)*(sind(phi))*(sind(D(i)))+((cosd(phi))*(cosd(D(i)))*(sind(ws(i))))];Hd(i)=(1.411-(1.686*Hg(i)/Ho(i)))*Hg(i); % Monthly average of daily diffuse radiationwst(i)=acosd(-(tand(phi-B))*(tand(D(i))));Rb(i)=(((sind(D(i)))*(sind(phi-B))*(wst(i)*pi/180))+((cosd(D(i)))*(cosd(phi-B))*(sind(wst(i)))))/(((ws(i)*pi/180)*(sind(phi))*(sind(D(i))))+((cosd(phi))*(cosd(D(i)))*(sind(ws(i)))));Rd(i)=(1+(cosd(B)))/2;Rr(i)=(1-(cosd(B)))*p/2;Ht(i)=(((1-(Hd(i)/Hg(i)))*Rb(i))+((Hd(i)/Hg(i))*Rd(i))+Rr(i))*Hg(i);a(i)=0.459+(0.5016*(sind(ws(i)-60)));b(i)=0.6609-(0.4767*(sind(ws(i)-60)));for j=1:12 WT(i,j)=j+6; EOT(i)=(9.87*sind(2*(360/364)*(n(i)-81)))-(7.53*cosd((360/364)*(n(i)-81)))-(1.5*sind((360/364)*(n(i)-81))); LAT(i,j)=WT(i,j)-(((4*(82.5-79.1833))+EOT(i))/60); w(i,j)=15*(LAT(i,j)-12);Io(i,j)=Isc*3600*(1+(0.033*(cosd(n(i)*360/365))))*(((sind(phi))*(sind(D(i))))+((cosd(phi))*(cosd(D(i)))*(cosd(w(i,j)))));Ig(i,j)=Io(i,j)*(a(i)+(b(i)*(cosd(w(i,j)))))*Hg(i)/Ho(i);if (Hd(i)/Hg(i))