Using 3 dimensional to plot 2-d plot? Urgent Help Required.

Hello

Using this equation

. to get this plot .

i did a approach bt output isnt coming. my codes is provided.

pmr = 180 which is the pi value.

This is my codes im working on.

clear all; close all; clf; clc

%Fractional Order PID
i=1;
j=1;
k=1;
pmr=180;  %Value of PI

i=1;
for w=0:1:100; % Omega Value Set Range

      j=1;
      for mu4=0:0.1:1; % Mu Value Set Range

          k=1;
          for lam4=0:0.1:1;  % Lambda Value Set Range

   % Short Hand calcalution for Theta values of cos and sin
              th1=cos(lam4*(pmr/2));
              th2=cos((1+lam4)*(pmr/2));
              th3=cos((lam4+mu4)*(pmr/2));
              th4=sin(mu4*(pmr/2));
              th5=sin(lam4*(pmr/2));

              ki(i,j,k)=((-w^(3+lam4))+(42.46*(w^(1+lam4)))+(100*w^(mu4+lam4)*th4))/(250*th5)

              kp(i,j,k)=(-1/(250*(w^lam4)*th1))*(((-w^(3+lam4))*th2)+((42.46*w^(1+lam4))*th2)...
              -(15.88*(w^(2+lam4))*th1)+(106.2*(w^lam4)*th1)+(100*(w^(lam4+mu4))*th3)...
              +250*ki)   

              i=i+1; % for next update

          end % End loop for Lambda

          j=j+1; % for next update

      end % End loop for mu

      k=k+1; % for next update

end % End loop for Omega

figure(1)
plot(ki,kp,'r', 'Linewidth',0.5)

Im using the above coding using 3 dimensional array for loop. However, im having a error which i cant resolve..

Assignment has more non-singleton  rhs dimensions than non-singleton  subscripts

Your help will be appreciated. Thank you.

Nitesh

[Walter Roberson]

http://www.mathworks.com/matlabcentral/answers/29922-why-your-question-is-not-urgent-or-an-emergency