Deleted The Following Issue was solved.

I have the following code , i'm trying to convert from Cartesian coordinates to polar but so far when i plot , the plot is empty.

    Latitude=31.13;
    for n=1:30:365/2;
    for h=-12:0.1:12
        w=15*h;
        delta=23.45*sind((360/365)*(284+n)); %Declination Angle
        alpha=asind(sind(delta)*sind(Latitude)+cosd(delta)*cosd(w)*cosd(Latitude));%Solar Altitude Angle
        gamma1=asind(cosd(delta)*sind(w)/(cosd(alpha)));%Azimuth Angle Dashed
        we=acosd(tand(delta)/tand(Latitude));
        if abs(w)<we
            C1=1;
        else 
            C1=-1;
        end 
        if Latitude*(Latitude-delta)>=0
            C2=1;
        else
            C2=-1;
        end
        if w>=0
           C3=1;
        else
            C3=-1;
        end
        Gamma=C1*C2*gamma1+C3*((1-C1*C2)/2)*180;
        [phi,r]=cart2pol(Gamma,alpha);
        polar(phi,r)
    end
    end

[Star Strider]

You have to create vectors (or matrices) out of ‘Gamma’ and ‘alpha’, but since you have two loops, and I have no idea what result you want, I have no idea what to suggest. It’s straightforward to reference them by ‘n’, but if you want to reference them by ‘h’, you need to use a counter that increments with each value of ‘h’, because ‘h’ begins with a negative value and has fractional increments.

The assignments would be:

Gamma_array(n) = Gamma;
alpha_array(n) = alpha;

or:

Gamma_array(n,k) = Gamma;
alpha_array(n,k) = alpha;

where ‘k’ is the ‘h’ counter.

Assign your vectors (or matrices) for ‘Gamma’ and ‘alpha’ at the end of the loop, then put your car2pol and polar calls outside and after both loops. That will also make your code considerably faster.

[Star Strider]

I apologise for repeating myself, but as I wrote in my first Comment: ‘but if you want to reference them by ‘h’, you need to use a counter that increments with each value of ‘h’, because ‘h’ begins with a negative value and has fractional increments’

So I would create the counter as:

      for n=1:30:365/2;
      k1 = 0;
      for h=-12:0.1:12
          k1 = k1 + 1;
          ... CODE ...
          Gamma_array(n,k1) = Gamma;
          alpha_array(n,k1) = alpha;
      end
      end

See if that does what you want.

[Rena Berman]

(Answers Dev) Restored question.