I have this code and I want to calculate the maximum value of alpha but it's giving me the minimum value:
phi = 53;h = 1/24:1/24:1;d = 100:100:365;for i = 1:1:length(d) delta(i) = -23.45 * cosd((360/365) * (d(i)+10)); for j = 1:1:length(h) h_a(j) = 360 * (h(j)-0.5); alpha(i,j) = asind(((sind(phi) * sind(delta(i))) + (cosd(phi) * cosd(delta(i)) * cosd(h_a(j))))); theta(i,j) = 90 - alpha(j); if h <= 0.5 beta(i,j) = -acosd((((sind(delta(i)) * cosd(phi)) - (cosd(h_a(j)) * cosd(delta(i)) * sind(phi)))/sind(theta(j)))); else beta(i,j) = acosd((((sind(delta(i)) * cosd(phi)) - (cosd(h_a(j)) * cosd(delta(i)) * sind(phi)))/sind(theta(j)))); end endenda = max(alpha(i,j));fprintf('%.1f\n',a);disp(a)
Best Answer