MATLAB can numerically solve the equilibrium of your model using the 'fsolve' function. You need to write a function that returns the steady state equations for your model and then use 'fsolve' to find the roots. Here are the steps to solve for equilibrium points for your model:
- Set up your system and parameters
- Create a function that returns the vector of equations at steady state
- Choose a reasonable initial guess and supply the function handle and initial guesses to the 'fsolve' function
x0 = [1, 1, 1, 1, 1]; % Initial guess
F = @(x) glioma_equilibrium(<Add parameter values>);
[x_equil, fval, exitflag] = fsolve(F, x0);
- Interpret the equilibria stored in the x_equil variable
- Repeat this process using different initial guesses to arrive at a meaningful solution
Since your model is non-linear, it may have multiple equilibria and hence the initial guess you provide to the 'fsolve' function will determine the chances of convergence. To ensure that you receive meaningful equilibria, make your parameters and initial guesses physically meaningful.
For more information regarding the 'fsolve' function, refer to the following documentation: https://www.mathworks.com/help/optim/ug/fsolve.html
