Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Upload

Welcome to Scribd!

  • 46 views

    Numerical Analysis: MATLAB Practical (Autumn 2020) B.E. III Semester Thapar Institute of Engineering & Technology Patiala

    Uploaded by

    Aarohan Verma
    This document presents the numerical analysis of solving ordinary differential equations using the Modified Euler's Method and the Runge-Kutta Fourth Order Method. It provides the algorithms, MATLAB code, and an input/output table comparing the results of applying each method to equations for problems y'=-y+2cos(x) and y'=(x+y) with initial conditions.

    Copyright:

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

    Numerical Analysis: MATLAB Practical (Autumn 2020) B.E. III Semester Thapar Institute of Engineering & Technology Patiala

    Uploaded by

    Aarohan Verma
    46 views6 pages
    This document presents the numerical analysis of solving ordinary differential equations using the Modified Euler's Method and the Runge-Kutta Fourth Order Method. It provides the algorithms, MATLAB code, and an input/output table comparing the results of applying each method to equations for problems y'=-y+2cos(x) and y'=(x+y) with initial conditions.

    Original Title

    TIET_MATLAB4

    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 the numerical analysis of solving ordinary differential equations using the Modified Euler's Method and the Runge-Kutta Fourth Order Method. It provides the algorithms, MATLAB code, and an input/output table comparing the results of applying each method to equations for problems y'=-y+2cos(x) and y'=(x+y) with initial conditions.

    Copyright:

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

    Numerical Analysis: MATLAB Practical (Autumn 2020) B.E. III Semester Thapar Institute of Engineering & Technology Patiala

    Uploaded by

    Aarohan Verma
    This document presents the numerical analysis of solving ordinary differential equations using the Modified Euler's Method and the Runge-Kutta Fourth Order Method. It provides the algorithms, MATLAB code, and an input/output table comparing the results of applying each method to equations for problems y'=-y+2cos(x) and y'=(x+y) with initial conditions.

    Copyright:

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

    NUMERICAL ANALYSIS

    MATLAB Practical (Autumn 2020)


    B.E. III Semester
    Thapar Institute of Engineering & Technology Patiala

    Name: Aarohan Verma


    Roll No.: 101903497
    Group: COE19
    Instructor: Mrs. Mamta Gulati
    (ASSIGNMENT 4- Modified Euler’s and
    Runge Kutta Fourth Order Method)

    Algorithms:
    Modified Euler’s Method:
    Given: a, b, f(x, y), either h or N, y(a) = yo
    Step 1: Assign x(1) = a, y(1) = yo and find h or N and f(x, y)
    Step 2: Assign 1 to i
    Start the loop and continue till i<=N
    Assign h*f(x(i), y(i)) to k1 and h*f(x(i)+h, y(i) + k1) to k2
    Assign x(i) + h to x(i+1) and y(i) + (k1+k2)/2 to y(i+1)
    Increment the value of i by 1
    End of loop
    Step 3: Print x and y array.

    Runge-Kutta Order Four Method:


    Step 1: Assign x(1) = a, y(1) = yo and find h or N and f(x, y)
    Step 2: Assign 1 to i
    Start the loop and continue till i<=N
    Assign h*f(x(i), y(i)) to k1 and h*f(x(i)+h/2, y(i) + k1/2) to k2
    Assign h*f(x(i) + h/2, y(i) + k2/2) to k3 and h*f(x(i)+h, y(i) + k3) to k4
    Assign x(i) + h to x(i+1) and y(i) + (k1+2*k2 + 2*k3 + k4)/6 to y(i+1)
    Increment the value of i by 1
    End of loop
    Step 3: Print x and y array.
    CODE:

    Modified Euler’s Method:

    clear all
    clc
    f = @(x,y) (-y+2*cos(x));
    a = 0;
    b = 1;
    h = 0.2;
    N = (b-a)/h;
    x(1) = a;
    y(1) = 1;
    for i=1:N
    k1 = h*f(x(i), y(i));
    k2 = h*f(x(i)+h, y(i)+k1);
    x(i+1) = x(i) + h;
    y(i+1) = y(i) + (k1+k2)/2;
    end
    fprintf("x y\n");
    for i=1:N
    fprintf("%f %f\n",x(i),y(i));
    end
    Runge-Kutta Order Four Method:

    clear all
    clc
    f = @(x,y) (-y+2*cos(x));
    a = 0;
    b = 1;
    h = 0.2;
    N = (b-a)/h;
    x(1) = a;
    y(1) = 1;
    for i=1:N
    k1 = h*f(x(i), y(i));
    k2 = h*f(x(i)+h/2, y(i)+k1/2);
    k3 = h*f(x(i)+h/2, y(i)+k2/2);
    k4 = h*f(x(i)+h, y(i)+k3);
    x(i+1) = x(i) + h;
    y(i+1) = y(i) + (k1 + 2*k2 + 2*k3 + k4)/6;
    end
    fprintf("x y\n");
    for i=1:N+1
    fprintf("%f %f\n",x(i),y(i));
    end
    Input & Output Table:

    Input Modified Euler’s Runge-Kutta Order Four


    Method Method
    3(a) x y x y
    y’ = -y +2*cos(x) 0.000000 1.000000 0.000000 1.000000
    y(0) = 1 0.200000 1.176013 0.200000 1.178732
    0.400000 1.305354 0.400000 1.310472
    0.600000 1.382827 0.600000 1.389968
    0.800000 1.405313 0.800000 1.414051
    1.000000 1.371890 1.000000 1.381759

    3(b) x y x y
    y’ = √𝑥 + 𝑦 0.000000 0.800000 0.000000 0.800000
    0.200000 0.998019 0.200000 0.999291
    y(0) = 0.8 0.400000 1.234632 0.400000 1.236964
    0.600000 1.507064 0.600000 1.510332
    0.800000 1.813386 0.800000 1.817506
    1.000000 2.152153 1.000000 2.157066

    You might also like