% Berechnung der Werte für die Hauptfunktion
% Rückgabe der Werte für die Ausgabe

function dc = berechnung( t,c,p )

% Übergabe der Werte
cA=c(1);
cB=c(2);
cC=c(3);

r1=p.k1*cA*cB^2;

dcA=p.FAzu/p.VR*p.CAzu-p.Fab/p.VR*cA-r1;

dcB=p.FBzu/p.VR*p.CBzu-p.Fab/r.VR*cB-2*r1;

dcC=(-p.Fab/p.VR)+cC+r1;


dc=[dcA; dcB; dcC];

dcdt=dcdt'; % transponiert den lösungsvektor da 
            % Rückgabewert ein Spaltenvektor sein muss

end