Guten Tag zusammen,
Wir haben in der Schule folgende Aufgabe bekommen:
"Erdbebenwellen pflanzen sich längs der Erdoberfläche
mit 3,3 km/s, zugleich aber gradlinig durch das Erdinnere mit 10 km/s fort.
Die erste Welle kommt im Ort A t1 Sekunden, im Ort B t2 Sekunden früher an als die zweite.
Aus diesen Angaben ist der Ort des Bebens zu bestimmen.
(von der Tiefe des Bebens unter der Erdoberfläche wird abgesehen)"
Wir mussten diese Aufgabe mithilfe von Mathematica lösen, natürlich kam ich auch auf die Lösung
-------------------------------------------------------------------------------------------------
Beben1[t1_, t2_, lan1_, phi1_, lan2_, phi2_, lan4_, phi4_] :=
Block[{},
a = FindRoot[10*(((6371*e)/3.3) - t1) - (2*6371*Sin[(e/2)]) == 0, {e, 1}];
e1 = e /. a[[1]];
b = FindRoot[10*(((6371*d)/3.3) - t2) - (2*6371*Sin[(
d/2)]) == 0, {d, 1}];
d1 = d /. b[[1]];
c = Table[FindRoot[{Cos[e1] == Sin[ phi3]*Sin[phi1] +
Cos[phi3]*Cos[phi1]*Cos[ lan3 - lan1], Cos[d1] ==
Sin[phi3]*Sin[phi2] + Cos[phi3]*Cos[phi2]*Cos[lan3 - lan2]},
{phi3, phi4}, {lan3, lan4 + phi2}]];
Print[c]]
----------------------------------------------------------------------------------------------
zum Beispiel:
Beben1[32*60 + 23, 32*60 + 9, (-18 - 3.5/60)°, (59 +
20.5/60)°, (-12 - 29/60)°, (41 + 54/60)°, -1.3, -0.5]
{phi3 -> -0.0346761, lan3 -> -1.57049}
(phi4 und lan4 sind nährungswerte, da findroot zu wenig gut ist :schimpf)
--------------------------------------------------------------------------------------------
Wüsste jemand eine Methode dieses Progrämmchen ein bischen schöner darzustellen??
gibt es evtl. eine Möglichkeit eine schönere Oberfläche zu erstellen?
Vielen Dank an die lieben Leute da draussen und noch schöne Festtage.
McClane