I have my code set up as follows:
thmaxn = ;
for ang=[0, pi/2, pi, 3*pi/2, 2*pi]
thmax = phasedist(ang,N,rhoss);
Don't worry what the function does, it just gives a value for each angle. I just wondered how I could put each value into the thmaxn vector successively?