function [MC,MK,uc,ub,nb] = Mod1_C1C2O3O4

% vorangestelltes Mod_ für Ausgabe von 6 Variablen (neu)

% stellt die Matrizen MC, MK und MV her
% MC enthält alle faktoriellen Ratekonstanten, Hauptdiagonalen noch Null
% MK gibt mit 1 an, wo die Konzentration (in M) als Faktor auftritt
% MV enthält Konstanten die mit der Spannung (in V) im Exponenten zu
% multiplizieren ist
% uc ist der Vektor der Offenzustände
% ub ist der Vektor der Bindungsgrade
% nb ist die maximala Anzahl an Bindungsstellen


ns = 4;           % Anzahl Zustände gleich 4
nb = 1;           % Anzahl Bindungsstellen gleich 1
uc = [0 0 1 1];   % Vektor der offenen Zustände - Strom fließt (current)
ub = [0 1 0 1];   % Vektor der gebundenen Zustände - Fluoreszenz 


% Vordefinierte Konstanten

Po0 = 0.0;          % wird hier nicht benötigt 

% Parameter, die i.a. aus dem Fitprogramm stammen

% p(1) = 1.0e-3;
% p(2) = 0.10e3;


% Ratekonstanten, werden aus den Parametern und/oder vordefinierten
% Konstanten berechnet, oder gleich als Zahl eingegeben

k12 = 3.0e7;
k21 = 4.22e2 ;
k34 = 3.0e7;
k43 = 22.4*10^-3;
k13 = 1.97e-7;
k31 = 6.32*10^-3;
k24 = 3.74*10^-3;
k42 = 11.7*10^-3;


% MC - Matrix, in welcher die Ratekonstanten stehen - sonst Nullen  
% MK - Matrix, mit 1 wenn zugehörige Ratekonstante kozentrationsabhängig
% - sonst Nullen
% MV - Matrix, mit Faktor bei Spannungsabhängigkeit im Exponenten
% - sonst Nullen

MC = zeros(ns,ns);
MK = MC;

MC(1,2)=k21;    % bekommt 1 von 2
MC(2,1)=k12;    % bekommt 2 von 1

MC(3,4)=k43;    % bekommt 3 von 4
MC(4,3)=k34;    % bekommt 4 von 3

MC(2,4)=k42;    % bekommt 2  von 4
MC(4,2)=k24;    % bekommt 4 von 2

MC(1,3)=k31;    % bekommt 1 von 3
MC(3,1)=k13;    % bekommt 3 von 1



MK(2,1) = 1;    % Übergang "bekommt 1 von 2" ist konzentrationsabhängig
MK(4,3) = 1;    % Übergang "bekommt 3 von 4" ist konzentrationsabhängig