function [dd] = elektron( t, Y)
% ======================================================================
% Beispiel einer rechten Seite einer Differentialgleichung
% trudelndes Elektron
%
% Erstellt fuer die Veranstaltung
% Weiterbildenes Studium fuer Ingenieure WS 98/99
%
% S.F.
% ----------------------------------------------------------------------

x1  = Y(1,:); % x
dx1 = Y(2,:); % x'
x2  = Y(3,:); % y
dx2 = Y(4,:); % y'
x3  = Y(5,:); % z
dx3 = Y(6,:); % z'

r=sqrt(x1.^2+x2.^2+x3.^2);

B1 = -3.*x1.*x3 ./ r.^5;
B2 = -3.*x2.*x3 ./ r.^5;
B3 = -(3.*x3.^2 - r.^2) ./ r.^5;

ddx1 = dx3.*B2 - dx2.*B3;
ddx2 = dx1.*B3 - dx3.*B1;
ddx3 = dx2.*B1 - dx1.*B2;

dd = [dx1' ddx1' dx2' ddx2' dx3' ddx3']';

return

%
%	Beispiel für eine Aufruf von startawa:
%
%     startawa('elektron',[0,50],[0.7,0,0,0.8,0,0.6],'nurdrei',[1,3,5])





