function traject(u0, u1, teta, omega, N, ns, nf, aRf)
n=N;
fi= 0.01: pi/6 : n*pi;
m=9.1*10^(-31);
e=1.6*10^(-19);
nse=1; if nf>(n-1), nf=(n-1); end;
ve=sqrt(2*e*u0/m);
m1=(sin(0.5*teta)/(0.5*teta));
hold off;
plot(fi(ns:nf),0.01*sin(fi(ns:nf)),'g-');
hold on
for i = 1:48
 v(i)=ve*(1+m1*u1/(2*u0)*sin(fi(i)));
 for j = 1:(n-1)
  z(j)=v(i).*(fi(j)/omega-fi(i)/omega)-0.5*aRf*(fi(j)/omega-fi(i)/omega).^2;
 end;
 if ns>i, nse=ns; else nse=i;end;
 plot(fi(nse:nf),z(nse:nf),'g-');
end;
grid on;
hold off;