Scribd Logo
Upload

Welcome to Scribd!

  • 4 views

    Q3 SQ7 Ma20011 Kuk2443 01G 1 PDF

    Uploaded by

    lim chuan yi
    This document presents data and calculations for solving the Van der Pol oscillator equation using Euler's method and the 4th order Runge-Kutta method. The data shows the values of y1, y2, and time t over 5 steps of 0.1 seconds each using Euler's method. It then shows the calculations using the 4th order Runge-Kutta method to solve the same system over the same time steps and interval.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd

    Q3 SQ7 Ma20011 Kuk2443 01G 1 PDF

    Uploaded by

    lim chuan yi
    4 views4 pages
    This document presents data and calculations for solving the Van der Pol oscillator equation using Euler's method and the 4th order Runge-Kutta method. The data shows the values of y1, y2, and time t over 5 steps of 0.1 seconds each using Euler's method. It then shows the calculations using the 4th order Runge-Kutta method to solve the same system over the same time steps and interval.

    Original Title

    Q3_SQ7_MA20011_KUK2443_01G__1_.pdf

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    This document presents data and calculations for solving the Van der Pol oscillator equation using Euler's method and the 4th order Runge-Kutta method. The data shows the values of y1, y2, and time t over 5 steps of 0.1 seconds each using Euler's method. It then shows the calculations using the 4th order Runge-Kutta method to solve the same system over the same time steps and interval.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    4 views4 pages

    Q3 SQ7 Ma20011 Kuk2443 01G 1 PDF

    Uploaded by

    lim chuan yi
    This document presents data and calculations for solving the Van der Pol oscillator equation using Euler's method and the 4th order Runge-Kutta method. The data shows the values of y1, y2, and time t over 5 steps of 0.1 seconds each using Euler's method. It then shows the calculations using the 4th order Runge-Kutta method to solve the same system over the same time steps and interval.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 4

    euler's method

    data given
    i t h y1 y2
    0 0 0.1 52.2900 83.8200 d/dt y1 = -5y1 + 3y2
    1 0.1 0.1 51.2910 63.9030 d/dty2 = y1 - 3y2
    2 0.2 0.1 44.8164 49.8612
    3 0.3 0.1 37.3666 39.3845
    4 0.4 0.1 30.4986 31.3058
    5 0.5 0.1 24.6410 24.9639

    RK 4TH ORDER METHOD

    i t y1 y2 k1,1 k1,2 y1 + (1/2)*k1,1*h y2 + (1/2)*k1,2*h k2,1 k2,2 y1 + (1/2)*k2,1*h y2 + (1/2)*k2,2*h k3,1 k3,2 y1 + k3,1*h y2 +k3,2*h k4,1 k4,2 h
    0 0 52.2900 83.8200 -9.9900 -199.1700 51.7905 73.8615 -37.3680 -169.7940 50.4216 75.3303 -26.1171 -175.5693 49.6783 66.2631 -49.6022 -149.1109 0.1
    1 0.1 49.1806 66.5032 -46.3935 -150.3290 46.8610 58.9868 -57.3445 -130.0993 46.3134 59.9982 -51.5723 -133.6813 44.0234 53.1351 -60.7118 -115.3818 0.1
    2 0.2 43.7650 53.2820 -58.9789 -116.0810 40.8160 47.4780 -61.6463 -101.6178 40.6827 48.2011 -58.8100 -103.9207 37.8840 42.8899 -60.7501 -90.7858 0.1
    3 0.3 37.7543 42.9829 -59.8226 -91.1945 34.7632 38.4232 -58.5461 -80.5065 34.8270 38.9576 -57.2620 -82.0459 32.0281 34.7784 -55.8053 -72.3070 0.1
    4 0.4 31.9669 34.8395 -55.3159 -72.5516 29.2011 31.2119 -52.3697 -64.4347 29.3484 31.6178 -51.8887 -65.5049 26.7780 28.2890 -49.0230 -58.0890 0.1
    5 0.5 26.7526 28.3308 -48.7706 -58.2399 24.3141 25.4188 -45.3139 -51.9424 24.4869 25.7337 -45.2335 -52.7142 22.2293 23.0594 -41.9681 -46.9490 0.1

    This study source was downloaded by 100000849462533 from CourseHero.com on 05-07-2023 12:17:58 GMT -05:00
    28/12/21, 12:25 PG Octave Online · Cloud IDE compatible with MATLAB

    octave:1> function yp =vanderpol(t,y)

    > > yp =[-5*y(1)+3*y(2);y(1)-3*y(2)];

    > > end

    octave:2> tspan=[0,0.5];

    octave:3> y0=[52.29,83.82];

    octave:4> h=0.1;

    octave:5> [t,y]=ode45("vanderpol";tspan,y0,h)

    error: parse error:

    syntax error

    >>> [t,y]=ode45("vanderpol";tspan,y0,h)

    octave:5> [t,y]=ode45("vanderpol",tspan,y0,h)

    t =

    0.0500

    0.1000

    0.1500

    0.2000

    0.2500

    0.3000

    0.3500

    0.4000

    0.4500

    0.5000

    y =

    52.290 83.820

    51.193 74.551

    49.194 66.498

    46.642 59.461

    43.780 53.277

    40.782 47.817

    37.766 42.978

    34.814 38.675

    31.976 34.836

    29.285 31.404

    26.759 28.328

    octave:6> h=0.1

    h = 0.1000

    octave:7> plot(t,y(:,1))

    This study source was downloaded by 100000849462533 from CourseHero.com on 05-07-2023 12:17:58 GMT -05:00

    https://octave-online.net/# 1/3
    https://www.coursehero.com/file/127545968/Q3-SQ7-MA20011-KUK2443-01G-1pdf/
    28/12/21, 12:25 PG Octave Online · Cloud IDE compatible with MATLAB

    55

    50

    45

    40

    35

    30

    25
    0 0.1 0.2 0.3 0.4 0.5

    octave:8> plot(t,y(:,1),t,y(:,2))

    This study source was downloaded by 100000849462533 from CourseHero.com on 05-07-2023 12:17:58 GMT -05:00

    https://octave-online.net/# 2/3
    https://www.coursehero.com/file/127545968/Q3-SQ7-MA20011-KUK2443-01G-1pdf/
    28/12/21, 12:25 PG Octave Online · Cloud IDE compatible with MATLAB

    90

    80

    70

    60

    50

    40

    30

    20
    0 0.1 0.2 0.3 0.4 0.5

    This study source was downloaded by 100000849462533 from CourseHero.com on 05-07-2023 12:17:58 GMT -05:00

    https://octave-online.net/# 3/3
    https://www.coursehero.com/file/127545968/Q3-SQ7-MA20011-KUK2443-01G-1pdf/
    Powered by TCPDF (www.tcpdf.org)

    You might also like