i=n+2;
    h = i;
    k=0;
    z=2;
        for j= 1:n-1
            while (i<=h+n-k-1)
                if (i==h)
                    if (i==node_number-1)
                        y(i)=0;
                        theta(i)=0;
                        M(i) = Me;
                        nu(i) =PrandtlMeyer(M(i),gamma);
                        mu(i)= asin(1/M(i))*180/pi;
                        mI = tand((mu(i-n+k)-theta(i-n+k)+mu(i))/2);
                        end
                    else
                        y(i)=0;
                        theta(i)=0;
                        nu(i) = theta(i-n+k)+nu(i-n+k)+(1/(sqrt((M(i-n+k)^2)-1)-cotd(theta(i-n+k))));
                        M(i) = InversePrandtlMeyer(nu(i));
                        mu(i)= asin(1/M(i))*180/pi;
                        mI = tand((mu(i-n+k)-theta(i-n+k)+mu(i))/2);
                        disp(mI);

                        end
                    end
                end
                i = i+1;
            end
            k=k+1;
            h=i;
        end

2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844
2.1844