function gr05m012
% Fuer MATLAB und OCTAVE gleichermassen geeignet!
% In Anlehnung an die Loesung der Gruppe 05 (Trabs/Friedrich) 

clear; clc; clf;

dt=0.01;

k1=4;
k2=3;
c10 = 7;

t=0:dt:3;
c=f(t,k1,k2,c10);

plot(t,c(1,:),'b');
hold on
plot(t,c(2,:),'r');
hold on
plot(t,c(3,:),'m');
hold off
title('Analytische Loesung des Reaktionsystems')

% legend (Funktion haengt von der OCTAVE-Version ab!)

% ============================================================

function c = f(t,k1,k2,c10)
% Funktionswerte der Loesung

n=numel(t);
for i = 1:n
  c(1,i) = c10*exp(-k1*t(i));
  c(2,i) = (-1/(k1-k2)*k1*c10*exp(-t(i)*(k1-k2))+1/ ...
                                         (k1-k2)*k1*c10)*exp(-k2*t(i));
  c(3,i) = k2*k1*c10/(k1-k2)*(1/k1*exp(-k1*t(i))-1/k2*exp(-k2*t(i)))- ...
                                          k2*k1*c10/(k1-k2)*(1/k1-1/k2);
end

% =============================================================
