how to write a c program for numerical solution of swing equation of n dimensional system by RK4(Runge-Kutta 4th order method)

See more:

C++

C

The swing equation is of the form:
𝜃 ̇ = 𝜔
𝜔 ̇ =−𝛼𝜔+𝑃 −𝐾∗sin⁡(𝜃−𝜃𝑔𝑟𝑖𝑑) //Take 𝜃=x and 𝜔=y.


So take F1(x,y,t) = y
Second one is G1(x,y,t)= -alpha*y + p - k*sin(x-Xgrid) //alpha=0.1,p=1,k=8

So please take F1 & G1 for n dimensions(i.e. F1,G1,F2,G2,F3,G3.....Fn,Gn) and write the code for solving using C language.