clear; clf
n=288;
u0=1000;  u1=50;  f=2e+9; omega=2*pi*f;
teta=pi/2;
uRf=2000; sRf=0.05;
aRf=1.6*10^(-19)*(u0-u1)/sRf/9.1*10^(-31)
ns=1; nf=ns+172;
figure(gcf);
traject(u0, u1, teta, omega, n, ns, nf, aRf);