Question:differential quadrature method stuck at rk4 method coding after successfully derived system of odes using dqm

Question:differential quadrature method stuck at rk4 method coding after successfully derived system of odes using dqm

Maple 2022

im trying to write differential quadrature method i have successfully  generated system of equations using code but now im stuck with rk4 method to solve these equation in a loop kindly give me suggestions or modify the code

 >
 >
 >
 (1)
 >
 >
 (2)
 >
 (3)
 >
 >
 (4)
 >
 >
 u_1 = -u[1]*(-10.1340811913091*u[1]+15.4039041703445*u[2]-8.08708750877537*u[3]+3.92079823166652*u[4]-1.10353370192667*u[5])+.046910077030668 u_2 = -u[2]*(-1.92051204726391*u[1]-1.51670643433575*u[2]+4.80550130432862*u[3]-1.85711605328765*u[4]+.488833230558754*u[5])+.230765344947158 u_3 = -u[3]*(.602336319455661*u[1]-2.87077648466948*u[2]-1.11022302462516e-015*u[3]+2.87077648466942*u[4]-.602336319455684*u[5])+.5 u_4 = -u[4]*(-.488833230558737*u[1]+1.85711605328767*u[2]-4.8055013043286*u[3]+1.51670643433578*u[4]+1.92051204726404*u[5])+.769234655052844 u_5 = -u[5]*(1.1035337019266*u[1]-3.92079823166643*u[2]+8.08708750877511*u[3]-15.4039041703442*u[4]+10.1340811913086*u[5])+.95308992296933
 >
 >
 Initial condition for u_1: u_1(0) = 0.010000 Initial condition for u_2: u_2(0) = 0.020000 Initial condition for u_3: u_3(0) = 0.030000 Initial condition for u_4: u_4(0) = 0.040000 Initial condition for u_5: u_5(0) = 0.050000
 >
 >
 >