using 4th order runge kutta method to solved coupled euqations

any one can assist me ,that how we can apply 4th order Runge kutta method for 4 coupled equation?the system is given below;

dx/dt=−ax − eω + yz

dy/dt= by + xz

dz/dt= cz + fω − xy

dω/dt = dω – gz

a = 50, b =−16, c = 10, d = 0.2, e = 10, f = 16, g = 0.5

Step size 0.001 .

regards

faiz