4th order Runge-Kutta Method

I want to solve this equation for project

𝔉(𝑦1,𝑦2,…,𝑦𝑁+2,𝑡)=𝑄(𝑦1,𝑦2,…,𝑦𝑁+2,𝑡)−𝐶(𝑦1,𝑦3,…,𝑦𝑁+2)𝑦22𝐽(𝑦1,𝑦3,…,𝑦𝑁+2) using

Excess Differential Equations Method

Secondary coordinates are used in the equation of motion and N+2 first order differential equations are solved at each time step where N is the number of the secondary coordinates.

𝑦1=𝑞 ⇒ 𝑦̇1=𝑞̇=𝑦2

𝑦2=𝑞̇ ⇒ 𝑦̇2=𝑞̈=𝔉

𝑦3=𝜙1 ⇒ 𝑦̇3=𝜙̇1=𝑔1𝑞̇=𝑔1𝑦2

𝑦4=𝜙2 ⇒ 𝑦̇4=𝜙̇2=𝑔2𝑞̇=𝑔2𝑦2

: : :

𝑦𝑁+2=𝜙𝑁 ⇒ 𝑦̇𝑁+2=𝜙̇𝑁=𝑔𝑁𝑞̇=𝑔𝑁𝑦2

I hope for your help

[Arpit Bhatia]

Hi Adil,

Please include some information on what you have tried as of now and what portions do you specifically need MATLAB related help in. Refer to the following post on how to seek help on homework questions: https://www.mathworks.com/matlabcentral/answers/8626-how-do-i-get-help-on-homework-questions-on-matlab-answers