//Definição dos parâmetros
E0 = 100;
beta = 15;
w = 10;

lambda = 2*%pi/beta;
printf("Lambda: %.2f \n\n", lambda);

r0 = 1.5;
Dr = 0.00001;

printf("R0: %.2f", r0);

t0 = 7;

r = [-r0:Dr:r0];

E_minus = E0*exp(%i*(beta*r-w*t0));

Dt = 0.1;
t1 = t0 + Dt;

E_minus_t1 = E0*exp(%i*(beta*r-w*t1));


figure(1);
plot(r,E_minus,r,E_minus_t1);



E_plus = E0*exp(%i*(beta*r+w*t0));
E_plus_t1 = E0*exp(%i*(beta*r+w*t1));

figure(2);
plot(r,E_plus,r,E_plus_t1);