IMAG0151 (12)

Rys. 4. Trajektoria fazowa oscylatora harmonicznego uzyskana po wykonaniu

Ma ona kształt odpowiadający punktowi osobliwemu typu „środek co iwtadczy st ibilm sci rozpatrywanego układu przy zadanych parametrach.

Przedstawione tutaj m~pliki zostały zredagowane w możliwie najprostszy i ogólny sposób, celem łatwiejszego zrozumienia zasad ich budowania. Mogą one być przez analogię wykorzystywane do innych równań, wystarczy zmienić wartość parametrów inicjujących i opis równania różniczkowego. W dość łatwy sposób można spowodować, że rozwiązywanie równania będzie przebiegać interaktywnie. Wystarczy bowiem zamiast wpisywać na sztywno parametry w skrypt, zastosować funkcję w/prowadzania wartości zmiennej z klawiatury

nazwa zmiennej = input [Textpojawiając}’ się na ekranie:']

Takie rozwiązanie zostało zastosowane w drugim przykładzie, opisującym ruch ciężarka wiszącego na sprężynie. Współczynnik tłumienia b został zadeklarowany jako zmienna globalna i jest wprowadzony z klawiatury.

|% ciężar. in

I'- T. Salacinski & G. Sycz, Warszawa 1997

ciężar.m rozwiązuje numerycznie r. różniczkowe 11 rzędu

7