Numerische Integration partieller Differentialgleichungen, die physikalische Systeme mit endlicher Ausbreitungsgeschwindigkeit beschreiben, kann dadurch erfolgen, daß das ursprüngliche System mit Hilfe eines diskreten dynamischen Systems modelliert wird. Wenn das ursprüngliche System im eigentlichen physi- kalischen Sinn passiv ist, so läßt es sich durch eine Zeit-Raum-Koordinatentrans- formation in ein System transformieren, das mehrdimensional passiv ist, also passiv in einem verallgemeinerten, nämlich mehrdimensionalen Sinn. Entspre- chend kann dann auch das zugehörige diskrete System mehrdimensional passiv gestaltet werden. Dadurch gelingt es insbesondere, eine geeignete mehrdimensio- nale vektorielle Ljapunow-Funktion verfügbar zu machen. Die wichtigsten Vorteile, die das Verfahren für den sich ergebenden Algorith- mus liefert, sind: massiver Parallelismus, volle Lokalität aller Operationen, leichte Beherrschbarkeit der numerischen Stabilität, hohe Robustheit gegenüber den unvermeidbaren Rechenfehlern (Rundungs- bzw. Schneidefehler, Überlauf- korrekturen), die durch die Beschränktheit der auf einem Rechner zur Verfügung stehenden Wortlängen entstehen, sinnvolle Interpretationsmöglichkeit von Frequenzbereichs-Betrachtungen, Eignung als Grundlage für den Bau massiv paral- leler Spezialrechner. Die Anwendbarkeit des Verfahrens ist für die Akustik, Elektrodynamik, Elastizität und Fluiddynamik nachgewiesen worden.