로렌츠 방정식의 상태공간표현(4차 런지-쿠타 방정식)

본문내용

Lorenz 방정식의 상태-공간 표현 (4th Runge Kutta Methode)
방정식 초기값은 t = 0
dx/dt = -10x + 10y x = 5, y = 5, z = 5
dy/dt = 28x - y - xz
dz/dt = -2.6666667z + xy
Table 작성
t x y z k11 k12 k13
0.0 5.000000 5.000000 5.000000 0.000000 110.000000 11.666665
0.1 9.781470 17.079464 10.439470 72.979943 154.688333 139.223675
0.2 17.702968 20.874102 35.896877 31.711335 -160.672272 273.808548
0.3 10.810880 -2.529237 39.307446 -133.401170 -119.714197 -132.163148
0.4 0.549578 -5.544186 28.074618 -60.937644 5.503178 -77.912621
0.5 -3.164605 -5.841285 22.368887 -26.766793 -11.978966 -41.165012
0.6 -5.575876 -8.420372 19.923124 -28.444962 -36.615284 -6.177390
0.7 -8.887190 -12.678892 22.141485 -37.917012 -39.386842 53.635756
0.8 -11.914211 -13.430034 29.800004 -15.158230 34.875658 80.541574
0.9 -10.666799 -7.217843 33.399029 34.489556 64.808197 -12.072805
1.0 -6.846781 -3.430182 29.307167 34.165995 12.380068 -54.666752

키워드

4차 런지-쿠타,   4th Runge Kutta Methode,   나비효과