Implement solver for stiff system of differential equations

Reference no: EM13908797

In this problem set, you will implement a solver for the stiff system of differential equations

y' = -Ay,

y(₀) = y₀,

where A is a symmetric positive definite n x n matrix, y is a function from [0, ∞) to Rⁿ, and y₀ is a vector in Rⁿ. You can choose one of the following two algorithms to implement.

Method 1. Solve the differential equation using the backward Euler method. Then, apply one step of Richardson extrapolation.

Method 2. Diagonalize the matrix A. Then, find an explicit formula for the solution. In C, your calling sequence should be void stiff_solve(double* a, int n, double* y0, double t, int m, double* y) Input arguments:

a - The n x n symmetric positive definite matrix of doubles A

n - The dimension n of the matrix A

y₀ - An array of doubles of length n containing the initial value yo

t - The final time t

m - The number of steps used in backward Euler (the step size should be t/m). This argument should be ignored if you are using method 2.

Output arguments:

y - An array of doubles of length n containing the solution y(t) at time t

Attachment:- Assignment.rar

Reference no: EM13908797

Questions Cloud

Bonds would eliminate interest-rate risk of equity holders : As a newly hired CEO of the People Trust Co., the first you do is to study the firm’s balance sheet. You find that your firm has $100 million in three-year loans (total assets), $70 million in one-year deposits (total liabilities) and $30 million in ..

After evaluating null company manufacturing process : After evaluating Null Company's manufacturing process, management decides to establish standards of 3 hours of direct labor per unit of product and $ 15 per hour for the labor rate.

Issuing new common stock to raise investment funds : Issuing new common stock to raise investment funds is more expensive than using retained earnings as the funding source

Dividends are expected to grow at a constant rate : The Jackson–Timberlake Wardrobe Co. just paid a dividend of $1.20 per share on its stock. The dividends are expected to grow at a constant rate of 4 percent per year indefinitely. Investors require a return of 10 percent on the company's stock. What ..

Implement solver for stiff system of differential equations : Implement a solver for the stiff system of differential equations - Array of doubles of length n containing the initial value.

State of nature probability return : State of Nature Probability Return on A Return on B I 0.4 6% 11% II 0.2 9% 6% III 0.4 12% 15% The correlation between A and B is 0.35. The portfolio weight of A is 25%. The portfolio weight of B is 75%. 22.value: The expected return in percent for in..

Draw rate of product formation graph using acid phosphatase : In an acid phosphatase you have determined Ki to be 10.5 mu M for NaF, a non-competitive inhibitor. How much sodium fluoride is required to bring 65% inhibition of the Vmax?

What will be the level of accounts receivable : Harrelson Inc currently has $750,000 in accounts receivable and it's days sales outstanding (DSO) is 55 days. It wants to reduce its DSO to 35 days by pressuring more of its customers to pay their bills on time. If this policy is adopted the company'..

Sedona company set the following standard costs : Sedona Company set the following standard costs for one unit of its product for 2013.

User Account

All Pages