Here we aim to find the steady-state solution without solving an ODE to 'exhaustion'. ss 0=-BSssISS...

 

Transcribed Image Text:

Here we aim to find the steady-state solution without solving an ODE to 'exhaustion'. ss 0=-BSssISS +VRS 0= BSssIss - YssI 0=YISS-VRSS Now we need to solve these 3 equations for the 3 steady state unknowns, Sss, Iss, Rss. We will need a multivariable Newton solver for this such as broyden. or fsolve. (a) Solve Eqn. 6 for the steady-state. (b) It is a little disingenuous to use as starting values for the algebraic equation solver the values you found by 'integrating to exhaustion', since if you are going to bother with the former, then you may as well not bother with the latter. Repeat the algebraic solution procedure from part a, but this time use a quite dif- ferent starting condition, say ass = [2,0,55]T. Explain what happens, and if there are any issues how you are going to fix them. Here we aim to find the steady-state solution without solving an ODE to 'exhaustion'. ss 0=-BSssISS +VRS 0= BSssIss - YssI 0=YISS-VRSS Now we need to solve these 3 equations for the 3 steady state unknowns, Sss, Iss, Rss. We will need a multivariable Newton solver for this such as broyden. or fsolve. (a) Solve Eqn. 6 for the steady-state. (b) It is a little disingenuous to use as starting values for the algebraic equation solver the values you found by 'integrating to exhaustion', since if you are going to bother with the former, then you may as well not bother with the latter. Repeat the algebraic solution procedure from part a, but this time use a quite dif- ferent starting condition, say ass = [2,0,55]T. Explain what happens, and if there are any issues how you are going to fix them.

Answer rating: 100% (QA)

clear close all clc x0 6414054050 x fsolve fx0 7 Fu... View the full answer