clear,clc
clf
g=9.81;
Dure = 60;
N = 51;
for i= 1: 10
    eval(['L',num2str(i),'=g*((Dure/(2*pi*(N+',num2str(i),')))^2);']);
end
    
b=.005;
M = 1;

 AngleInit =pi/3;

Theta0_1 = AngleInit;
Vit0_1 = 0;
Theta0_2 = AngleInit;
Vit0_2 = 0;
Theta0_3 = AngleInit;
Vit0_3 = 0;
Theta0_4 = AngleInit;
Vit0_4 = 0;
Theta0_5 = AngleInit;
Vit0_5 = 0;
Theta0_6 = AngleInit;
Vit0_6 = 0;
Theta0_7 = AngleInit;
Vit0_7 = 0;
Theta0_8 = AngleInit;
Vit0_8 = 0;
Theta0_9 = AngleInit;
Vit0_9 = 0;
Theta0_10 = AngleInit;
Vit0_10 = 0;
axis([-L1 L1 -L1 0])
axis square
hold on
B1 = plot(L1*sin(Theta0_1),-L1*cos(Theta0_1),'ok');
B2 = plot(L2*sin(Theta0_2),-L2*cos(Theta0_2),'ok');
B3 = plot(L3*sin(Theta0_3),-L3*cos(Theta0_3),'ok');
B4 = plot(L4*sin(Theta0_4),-L4*cos(Theta0_4),'ok');
B5 = plot(L5*sin(Theta0_5),-L5*cos(Theta0_5),'ok');
B6 = plot(L6*sin(Theta0_6),-L6*cos(Theta0_6),'ok');
B7 = plot(L7*sin(Theta0_7),-L7*cos(Theta0_7),'ok');
B8 = plot(L8*sin(Theta0_8),-L8*cos(Theta0_8),'ok');
B9 = plot(L9*sin(Theta0_9),-L9*cos(Theta0_9),'ok');
B10 = plot(L10*sin(Theta0_10),-L10*cos(Theta0_10),'ok');

set(B1,'markersize',10,'Linewidth',4);
set(B2,'markersize',10,'Linewidth',4);
set(B3,'markersize',10,'Linewidth',4);
set(B4,'markersize',10,'Linewidth',4);
set(B5,'markersize',10,'Linewidth',4);
set(B6,'markersize',10,'Linewidth',4);
set(B7,'markersize',10,'Linewidth',4);
set(B8,'markersize',10,'Linewidth',4);
set(B9,'markersize',10,'Linewidth',4);
set(B10,'markersize',10,'Linewidth',4);

pause(2)
Tr = 0;
dt = 0.025;

while (Tr < Dure)
    sim('plume',[Tr Tr+dt]);
    Tr = Tr + dt;
    fprintf('Temps écoulé %2.2f\n',Tr);
    Xend1 = x1(length(x1));
    Yend1 = y1(length(y1));
    set(B1,'Xdata',Xend1,'Ydata',Yend1);
    Xend2 = x2(length(x2));
    Yend2 = y2(length(y2));
    set(B2,'Xdata',Xend2,'Ydata',Yend2)
    Xend3 = x3(length(x3));
    Yend3 = y3(length(y3));
    set(B3,'Xdata',Xend3,'Ydata',Yend3)
    Xend4 = x4(length(x4));
    Yend4 = y4(length(y4));
    set(B4,'Xdata',Xend4,'Ydata',Yend4)
    Xend5 = x5(length(x5));
    Yend5 = y5(length(y5));
    set(B5,'Xdata',Xend5,'Ydata',Yend5)
    Xend6 = x6(length(x6));
    Yend6 = y6(length(y6));
    set(B6,'Xdata',Xend6,'Ydata',Yend6)
    Xend7 = x7(length(x7));
    Yend7 = y7(length(y7));
    set(B7,'Xdata',Xend7,'Ydata',Yend7)
    Xend8 = x8(length(x8));
    Yend8 = y8(length(y8));
    set(B8,'Xdata',Xend8,'Ydata',Yend8)
    Xend9 = x9(length(x9));
    Yend9 = y9(length(y9));
    set(B9,'Xdata',Xend9,'Ydata',Yend9)
    Xend10 = x10(length(x10));
    Yend10 = y10(length(y10));
    set(B10,'Xdata',Xend10,'Ydata',Yend10)
    drawnow;
    Theta0_1 = Theta1;
    Vit0_1 = Vit1;
    Theta0_2 = Theta2;
    Vit0_2 = Vit2;
    Theta0_3 = Theta3;
    Vit0_3 = Vit3;
    Theta0_4 = Theta4;
    Vit0_4 = Vit4;
    Theta0_5 = Theta5;
    Vit0_5 = Vit5;
    Theta0_6 = Theta6;
    Vit0_6 = Vit6;
    Theta0_7 = Theta7;
    Vit0_7 = Vit7;
    Theta0_8 = Theta8;
    Vit0_8 = Vit8;
    Theta0_9 = Theta9;
    Vit0_9 = Vit9;
    Theta0_10 = Theta10;
    Vit0_10 = Vit10;
   
end