function [dd] = dreikoerper( t, Y)
% ======================================================================
% Beispiel einer rechten Seite einer Differentialgleichung
% Drei-Körper-Problem
%
% Erstellt fuer die Veranstaltung
% Weiterbildenes Studium fuer Ingenieure WS 98/99
%
% S.F.
% ----------------------------------------------------------------------

x  = Y(1,:); % x
dx = Y(2,:); % x'
y  = Y(3,:); % y
dy = Y(4,:); % y'

mu=1/82.45;
mus=1-mu;

A=((x+mu ).^2+y.^2).^(1.5);
B=((x-mus).^2+y.^2).^(1.5);

ddx = x + 2*dy - mus*(x+mu)./A -mu*(x-mus)./B; %  x''
ddy = y - 2*dx - mus*(y)./A -mu*(y)./B;        %  y''

dd = [dx' ddx' dy' ddy']';

return

%
%	Beispiel fuer einen Aufruf der Datei startawa
%
%  startawa('dreikoerper',[0,6.19216933140],[1.2,0,0,-1.04935750983],'nur',[1,3])








