clear all;
close all;
clc;

Zeit = [0 4200];              %70 min
sigma= [0];                 %Anfangswert
[Tout,SigmaOut] = ode45(@dgl,Zeit,[0])
plot(Tout,SigmaOut)