function delta=rksteg(t,U,dt,funktionen)
k1=dt*funktionen(t,U);
k2=dt*funktionen(t+dt/2,U+k1/2);
k3=dt*funktionen(t+dt/2,U+k2/2);
k4=dt*funktionen(t+dt,U+k3);
delta=(k1+2*k2+2*k3+k4)/6;