Numerische Lösung des mathematischen Pendels mit ...

46 Export["hamExakt.plot", N[h[r]], "List"]
47
48 (*** Jetzt wollen wir noch die exakte Lösung *)
49 (*** mit den Näherungen vergleichen. *)
50
51 (** Unsere Differenz-Funktion *)
52 d[t_, x_] := { t, x - f[t] }
53
54 dateien := FileNames[{"euler*.plot", "banach*.plot"}]
55
56 (*** Jetzt eine Schleife über alle Dateien ... *)
57 Do[
58
59 name = dateien[[i]];
60 liste = Import[name, "List"];
61
62 (* Das ist jetzt eine Liste mit abwechselnd
63 Zeit und Näherung dazu *)
64
65 liste = Partition[liste, 2];
66
67 (* Jetzt haben wir eine Liste aus
68 zweielementigen Listen, jeweils Zeit und
69 Näherung *)
70
71 (*** Ersetze die Einzel-Listen (Tiefe: 1)
72 durch eine Anwendung von d. *)
73 liste = Apply[d, liste, {1}];
74
75 (* Jetzt haben wir eine Liste aus
76 zweielementigen Listen, jeweils Zeit
77 und Differenz *)
78
79 (* Beim Dateinamen wird einfach vorne
80 ein D angehängt. **)
81
82 Export["D" name, liste, "List"]
83
84 (* Das ganze machen wir für alle i
85 von 1 bis Anzahl der Dateien. *)
86
87 , {i, Length[dateien] }]
88
89 (* Wieder zurückwechseln *)
21