Scribd
Upload
55 views1 page

Che 560 HW5

This document provides instructions for homework assignment 5 regarding numerical methods for solving differential equations. Students are asked to modify an existing MATLAB code to solve a system of three coupled differential equations using three different methods: 1) an explicit Adams-Bashford 2nd-order method, 2) a 4th-order Runge-Kutta method, and 3) the built-in Ode45 solver. The document provides the equations, initial conditions, parameter values, and algorithm details for each method. Students are to submit three separate MATLAB code files, one for each method.

Uploaded by

Syed Samsamuddin
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
55 views1 page

Che 560 HW5

This document provides instructions for homework assignment 5 regarding numerical methods for solving differential equations. Students are asked to modify an existing MATLAB code to solve a system of three coupled differential equations using three different methods: 1) an explicit Adams-Bashford 2nd-order method, 2) a 4th-order Runge-Kutta method, and 3) the built-in Ode45 solver. The document provides the equations, initial conditions, parameter values, and algorithm details for each method. Students are to submit three separate MATLAB code files, one for each method.

Uploaded by

Syed Samsamuddin
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 1

King Fahd University of Petroleum & Minerals

Chemical Engineering Department


CHE 560 –Numerical Methods in Chemical Engineering
2018 - 2019 (182)

HW#5
Due: Tue. 19-Feb-2019

For the problem you have solved in the previous assignment:

𝑑𝑦1
= −𝑦2 𝑦10.5 − ∅ 𝑦3 𝑦10.5
𝑑𝑥

𝑑𝑦2
= −𝑦2 𝑦10.5
𝑑𝑥

𝑑𝑦3
= 𝑦2 𝑦10.5 − ∅ 𝑦3 𝑦10.5
𝑑𝑥

𝑥 = 0, 𝑦1 = 1, 𝑦2 = 0.5 𝑎𝑛𝑑 𝑦3 = 0

Using  = 0.55 and h = 0.1, modify the MATLAB code IVP.m to solve this problem for the
following three methods:

a) (20 points) 2nd-order explicit Adams-Bashford method, (send this code as


lastname_HW5_a.m):
𝑦𝑛+1 − 𝑦𝑛 3 1
= 𝑓𝑛 − 𝑓𝑛−1
ℎ 2 2

b) (20 points) 4th-oder Rung-Kutta method, (send this code as lastname_HW5_b.m):

yn 1  yn 
1
k1  k4   1 b k2  d k3 
6 3
k1  h f t n , yn 
 h 1 
k 2  h f  t n  , yn  k1 
 2 2 
 h 
k3  h f  t n  , yn  a k1  b k 2 
 2 
 h 
k 4  h f  t n  , yn  c k 2  d k3 
 2 
2 1 2 2  2 2
a , b , c , d  1
2 2 2 2

c) (10 points) Ode45 built-in solver, (send this code as lastname_HW5_c.m).

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy