%==========================================================Airplane===== f(z) = 1/z 
f0=figure(3);
clf;
dx=1.5;
dy=1.5;

du=10.5;
dv=10.5;
%---------main body
t=[0:600]/600*2*pi;
x=0.1*cos(t);
y=1*sin(t);
u=x./(x.^2+y.^2);
v=-y./(x.^2+y.^2);
f1=subplot(1,2,1);
plot(x,y,'r');
xlabel('x');
ylabel('y');
axis([-dx dx -dy dy]);
hold on;

f2=subplot(1,2,2);
plot(u,v,'r');
xlabel('u');
ylabel('v');
axis([-du du -dv dv]);
hold on;
pause;

%---------front wing
t=[0:600]/600*2*pi;
x=1*cos(t);
y=0.1*sin(t);
u=x./(x.^2+y.^2);
v=-y./(x.^2+y.^2);
subplot(1,2,1);
plot(x,y,'b');
xlabel('x');
ylabel('y');

subplot(1,2,2);
plot(u,v,'b');
xlabel('u');
ylabel('v');
pause;

%---------rear wing
t=[0:600]/600*2*pi;
x=0.4*cos(t);
y=-1+0.1*sin(t);
u=x./(x.^2+y.^2);
v=-y./(x.^2+y.^2);
subplot(1,2,1);
plot(x,y,'k');
xlabel('x');
ylabel('y');

subplot(1,2,2);
plot(u,v,'k');
xlabel('u');
ylabel('v');
pause;

%---------paddle 1
x=[-0.5:0.01:0.5];
y=x/2+1;
u=x./(x.^2+y.^2);
v=-y./(x.^2+y.^2);

subplot(1,2,1);
plot(x,y,'m');
plot(-x,y,'m');
xlabel('x');
ylabel('y');

m1(1)=getframe;

subplot(1,2,2);
plot(u,v,'m');
plot(-u,v,'m');
xlabel('u');
ylabel('v');

m2(1)=getframe(f0);

pause;

%---------paddle 2
subplot(1,2,1);
plot(x,y,'w');
plot(-x,y,'w');
subplot(1,2,2);
plot(u,v,'w');
plot(-u,v,'w');

x=[-0.7:0.01:0.7];
y=x./x;
u=x./(x.^2+y.^2);
v=-y./(x.^2+y.^2);
subplot(1,2,1);
plot(x,y,'m');
xlabel('x');
ylabel('y');
subplot(1,2,2);
plot(u,v,'m');
xlabel('u');
ylabel('v');
%pause;

y=[0.7:0.01:1.3];
x=1-y./y;
u=x./(x.^2+y.^2);
v=-y./(x.^2+y.^2);
subplot(1,2,1);
plot(x,y,'m');
xlabel('x');
ylabel('y');

m1(2)=getframe;

subplot(1,2,2);
plot(u,v,'m');
xlabel('u');
ylabel('v');

m2(2)=getframe(f0);

pause;

%-------------movie

%subplot(1,2,1);
%movie(m1,100);

%subplot(1,2,2);
movie(f0,m2,10,4,[0 0 0 0]);
