Malgo Joan

Thank you very much for explanation. I will try. Thanks again. On Friday 21 May 2021, 19:42:04 GMT+1, hilmeijer hilmeijer@users.sourceforge.net wrote: If you want to do bifurcation analysis on such a model, then formulate it as a return map, and use MatContM (the version for maps). Otherwise I'd keep it to plain Matlab for simulations. From: Malgo Joan joanna2007@users.sourceforge.net Sent: Friday, May 21, 2021 7:40 PM To: [matcont:discussion] Subject: [matcont:discussion] Re: Impulse function Hello,...

Hello, This is a dirac delta function, used to model the dose of the substance administered instantly (pharmacokinetics model). I have modelled a two dose regiment (one day apart). regards,

Hello, I have an impulse function in one of mine DE. Is it possible to use the MatCont for this ODE system and if yes, how should I write an impulse function? Are there any relevant examples? Thanks in advance for helping me.

Hello Hil, Thank you very much. Your guidance was the most useful. I know how to do it now. Thank you very much once again. Malgo On Tuesday 19 February 2019, 14:36:32 GMT, hilmeijer hilmeijer@users.sourceforge.net wrote: You add two additional equations u'=-wv+u(1-u^2-v^2) v'=wu+v(1-u^2-v^2) This is the Hopf-normal form with solution u=cos(w*t+phi), and phi corresponds to the initial condition. Now start with nonzero u,v, e.g. u=1,v=0, and then you can use u as periodic variable for your other 12...

Hello Hil, Thank you for your reply. I didn't check my private e-mail for a while and I've missed your reply.My ODE system is similar to the following ODE of 12 differential equations and contains 1 algebraic equation: L=2(18.311G^5-59.062G^4+68.983G^3-34.283G^2+5.9711G+0.6839); dxdt(1)=k_A(1/(1+(B/T_A_B)^n_A_B))-alpha_AA dxdt(2)=k_B((G/T_B_G)^n_B_G/(1+(G/T_B_G)^n_B_G))(1/(1+(D/T_B^D)^n_B_D))-alpha_B*B dxdt(3)=k_C((E/T_C_E)^n_C_E/(1+(E/T_C_E)^n_C_E))((D/T_C^D)^n_C_D/(1+(D/T_C^D)^n_C_D))+k_C_B(1/(1+(B/T_C^B)^n_C_B))-alpha_CC...

Hello hilmeijer, Thank you for your post. I've got through tutorial for bifurcation for periodic orbit (as you've kindly suggested): Codimension 2 bifurcations of periodic orbits in MatCont using my ODE system. However, I wasn't able to produce bifurcation diagram for chosen parameters. I wonder as in tutorial example was done for SIR model (2 equations) with 2 extra decoupled equations for u and v which I didn't add to my system. Are these equations required for matcont to produce bifurcation? If...

Hello, I'm new to Matcont. I have a problem with no convergence at x0 after clicking Type-Point and then using the last point in Type-Equilibrium. I try to find out how solution changes depending on the parameter value. I have ODE of 12 equations of periodic functions. I've tried to change initial condition values and Tolerance norms but with no luck. I came back to optimisation using lsqnonlin and I get local minimum possible info, however when playing with Algorithm type and Tolerance values, StepTolerance...