function planet_2()
clf;

% So stehen die Massen von Sonne, Erde und Mars in beiden Funktionen zur 
% Verfuegung
global m_S m_E m_M G
m_S = 1.998e30;     % Sonnenmasse
m_E = 5.974e24;     % Erdmasse
m_M = 6.419e23;     % Marsmasse
G = 6.673e-11;      % Gravitationskonstante
AE = ??;            % Mittlerer Abstand Sonne-Erde (Astronomische Einheit)
r_E_a = ??;         % Größter Abstand Sonne - Erde (Aphel)
v_E_a = ??;         % Geschwindigkeit der Erde im Aphel
r_M_a = ??;         % Größter Abstand Sonne - Mars (Aphel)
v_M_a = ??;         % Geschwindigkeit der Mars im Aphel

% Langweilige Anfangsbedingungen - finden Sie bessere
x_E_0 = 0.0;
y_E_0 = 0.0;
vy_E_0 = 0.0;
vy_E_0 = 0.0;
x_M_0 = 0.0;
y_M_0 = 0.0;
vx_M_0 = 0.0;
vy_M_0 = 0.0;

% Die Rechnung soll über 1000 Tage laufen. Die Genauigkeit soll
% +- 1 km bzw. +- 1 m/s sein
options = odeset ( 'RelTol', ??, 'AbsTol', [ ??, ??, ??, ??, ??, ??, ??, ?? ]);
[ T, Q ] = ode45( @pmotion, [ 0 ?? ], ...
                  [ x_E_0, y_E_0, vx_E_0, vy_E_0, ...
                    x_M_0, y_M_0, vx_M_0, vy_M_0 ], options );

r_E = ??        % Bahnradius der Erde in Abhängigkeit von der Zeit
                % gemessen in AE
r_M = ??        % Bahnradius des Mars in Abhängigkeit von der Zeit
                % gemessen in AE
phi_E = ?       % Phase der Erdposition in Abhängigkeit von der Zeit
phi_M = ?       % Phase der Marsposition in Abhängigkeit von der Zeit

% Was wird hier wohl geplottet?
subplot ( 2, 1, 1);
plot ( td, rE, '-r', td, rM, '--b', 'LineWidth', 2 );
set( gca(), 'FontSize', 20 );
xlabel ( 'Zeit [d]', 'FontSize', 20 ); 
ylabel ( 'Bahnradius [AE]', 'FontSize', 20 ); 
box on;
grid on;

subplot ( 2, 1, 2 );
plot ( td, phiE, '-r', td, phiM, '--b', 'LineWidth', 2 );
set( gca(), 'FontSize', 20 );
xlabel ( 'Zeit [d]', 'FontSize', 20 ); 
ylabel ( 'Phase [rad]', 'FontSize', 20 ); 
box on;
grid on;


end

function dq = pmotion ( t, q )
global m_S m_E m_M G

dq = zeros ( 8, 1 );

% Hier fehlen die Differentialgleichungen
dq(1) = 
dq(2) = 
dq(3) = 
dq(4) = 
dq(5) = 
dq(6) = 
dq(7) = 
dq(8) = 

end