Tmax=100; XY0=[10^7 0;0 0;10^7/2 10^7*sqrt(3)/2];
[A D E mymovie] = generator(Tmax,XY0,'ran');