Es fehlt noch s := s + DeltaS und t := t + DeltaT.
Wählen Sie als Startwert für DeltaT=0,01.