%load amp_data1
%load hdata
h=figure(1);
set(h,'Position', [200   300   250   400])

%a=[0.6793    0.2113
%    1.0000   -0.2146];

a=[ 0.6180    1.0000
    1.0000    -0.6180];



c1=(10)*[-1:.4:1];
c1=[c1 c1(5:-1:2)];
c1=[c1 c1 c1 c1 c1 c1 c1 c1 c1 c1 c1 c1 c1];
c2a=c1;
c2=(10)*[-1:(2/13):1];
c2=[c2 c2(12:-1:2)];
c2=[c2 c2 c2 c2 c2];
c1=c2;
c2=c2a;

b=a;
a=(1/max(b(:,1)))*a;
rdata=(1/sqrt(2))*[1 sqrt(2) 1 0 1];
tdata=[0 pi/4 2*pi/4 0 pi/4];
kl=length(c1);
   
axis image
axis off
view(25,3);
axis((3/4)*[-.5 4.5 -2 2 -8.5 1])
title('Even Mode Shapes');
Meven = moviein(kl);
j=0;
for i=1:kl
   cla

fprintf('%.0f out of %.0f frames for EVEN modes...\n',i,kl);
   tdata2a=a(1,1)*c1(i)*(pi/180)+tdata;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'k');
   hold on;
   tdata2a=a(1,1)*c1(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'b');
   tdata2a=a(1,1)*c1(i)*(pi/180)+tdata+pi;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'k');
   tdata2a=a(1,1)*c1(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'b');
   
   tdata2b=a(1,2)*c1(i)*(pi/180)+tdata;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'k'); 
   tdata2b=a(1,2)*c1(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'b'); 
   tdata2b=a(1,2)*c1(i)*(pi/180)+tdata+pi;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'k'); 
   tdata2b=a(1,2)*c1(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'b'); 


   tdata2a=a(1,2)*c1(i)*(pi/180)+tdata;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'k');
   tdata2a=a(1,2)*c1(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'b');
tdata2a=a(1,2)*c1(i)*(pi/180)+tdata+pi;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'k');
tdata2a=a(1,2)*c1(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'b');

   tdata2b=a(1,1)*c1(i)*(pi/180)+tdata;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'k');
   tdata2b=a(1,1)*c1(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'b');
   tdata2b=a(1,1)*c1(i)*(pi/180)+tdata+pi;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'k');
   tdata2b=a(1,1)*c1(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'b');




   tdata2a=a(2,1)*c2(i)*(pi/180)+tdata;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'k');
   tdata2a=a(2,1)*c2(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'r');
   tdata2a=a(2,1)*c2(i)*(pi/180)+tdata+pi;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'k');
   tdata2a=a(2,1)*c2(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'r');

   tdata2b=a(2,2)*c2(i)*(pi/180)+tdata;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'k'); 
   tdata2b=a(2,2)*c2(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'b');
   tdata2b=a(2,2)*c2(i)*(pi/180)+tdata+pi;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'k');
   tdata2b=a(2,2)*c2(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'b');
   
   tdata2a=a(2,2)*c2(i)*(pi/180)+tdata;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'k');
   tdata2a=a(2,2)*c2(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'b');
   tdata2a=a(2,2)*c2(i)*(pi/180)+tdata+pi;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'k');
   tdata2a=a(2,2)*c2(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'b');
   
   
   tdata2b=a(2,1)*c2(i)*(pi/180)+tdata;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'k');
   tdata2b=a(2,1)*c2(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'r');
   tdata2b=a(2,1)*c2(i)*(pi/180)+tdata+pi;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'k');
   tdata2b=a(2,1)*c2(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'r');

 plot3(-1+[0 1],[0 0],[-3.5 -3.5],'k')
  plot3(-1+[0 1],[0 0],[0 0],'k')
   plot3([0 1],[0 0],[-3.5 -3.5],'k')
   plot3([0 1],[0 0],[0 0],'k')
   plot3(1+[0 1],[0 0],[-3.5 -3.5],'k')
plot3(1+[0 1],[0 0],[0 0],'k')
plot3(2+[0 1],[0 0],[-3.5 -3.5],'k')
plot3(2+[0 1],[0 0],[0 0],'k')
plot3(3+[0 1],[0 0],[-3.5 -3.5],'k')
plot3(3+[0 1],[0 0],[0 0],'k')
plot3([-1 4 -1 4],[0 0 0 0],[0 0 -3.5 -3.5],'k.');

a=0-a;
plot3([0 0 1 2 3 3 0],[0 0 0 0 0 0 0],.2-2+.07*c1(i)*[0 a(1,1) a(1,2) a(1,2) a(1,1) 0  0],'b');
plot3([ 0 1 2 3 ],[ 0 0 0 0],.2-2+.07*c1(i)*[ a(1,1) a(1,2) a(1,2) a(1,1) ],'b.');


plot3([0  0 1 2 3 3 0],[0 0 0 0 0 0 0],-3.5-2+.07*c2(i)*[0 a(2,1) a(2,2) a(2,2) a(2,1)  0  0],'r');
plot3([0 1 2  3 ],[ 0 0 0 0],-3.5-2+.07*c2(i)*[a(2,1) a(2,2) a(2,2) a(2,1) ],'r.');
a=-a;

   
   axis image
   axis off
   view(25,3);
   axis((3/4)*[-.5 4.5 -2 2 -8.5 1])
   j=j+1;
title('Even Mode Shapes');
   Meven(:,j) = getframe;
end