I have a problem in runge kutta method for coupled differential equations .For example consider these equations:
The boundry conditions are:
Φ(x) and Ψ(x) are unknown functions and a,b,c,d,e,j,h are constants . I want to get answers by runge kutta numerical method .
i used this form of algorithm but i cant do it :
please help me to solve these types of equations.