Please, what is the right value that I am supposed to iterate in the "for" vectors "a" and "b" of line 12 and 13 in the picture attached? When I make the max limit to be more than 20, there is no output after some time waiting for calculation but I get an output when "a" is less than 30 and "b" is less than 40.

Note that vectors "As" and "as" all has 8760 X 1 entries.

location_filename = 'Delft.mat'; %You can replace 'Delft.mat' by: 'Buenos_Aires.mat'; 'New_Dehli.mat'; or 'Vancouver.mat' load(location_filename, 'As','as', 'DHI', 'DNI', 'GHI')% Loading vectors 'As', 'as', 'DHI', 'DNI' and 'GHI' 4 6 7 8 9 alpha = 0.2; Am = 0:2:360; theta = 6:2:90; am = 90-theta; % Albedo coefficient % Initializing PV Module Azimuth range from 0-360 degrees in steps of 2 degrees as a row vector % Initializing PV Module tilt angle range from 0-90 degrees in steps of 2 degrees as a row vector % PV Module altitude range as a row vector le 11 % Create 12 for a = 1:1: ... %'a' is the index to iterate through vector 'am' for b = 1:1: ... %'b' is the index to iterate through vector 'Am' + (sind(a). *sind(As)); % Calculate the cosine of the angle of incidence for e % Direct Irradiance COS_AOI = (cosd(a). *cosd(as). *cosd(b-As)) disp(cos_AOI); %direct = DNI .* COS_AOI; , ) estoas) + (santo), sindas)