Implementation of a differential equation (Shooting method) in Python by varying the initial and eigenvalues

def sol(q_x, kappa, h, x0, epsilon):

    y = Symbol('y')

    defsol(x,y,epsilon):
        return f(x,y,q_x,kappa,h,epsilon)

    y = np.linspace (0.0, 1000, 10001)

    return order (fsol, x0, y, args=(epsilon,))

alpha=np.linspace(0,2*np.pi,101)
X = np.cos (alpha)
Y = np.sin(alpha)
ZEROS = [[0.0]*2 for i in range (len(alpha))]
X0 = np.c_ [X, Y, ZEROS ]

def disposition0(q_x, kappa, h, cutoff):

    defz(x0, epsilon):
        return sol(q_x, kappa, h, x0, epsilon)

    epsilon=np.linspace(0,q_x**2+kappa+h,101)

    Phi1 = np.empty ([alpha.shape[0], epsilon.shape[0]])
    Phi2 = np.empty ([alpha.shape[0], epsilon.shape[0]])
    for i in range(0,len(alpha)):
        for jin range(0,len(epsilon)):
            Phi1[i][j]=z(X0[i], epsilon[j])[1000][0]
            Phi2[i][j]=z(X0[i], epsilon[j])[1000][1]
    energy = 0
    forkin range(0,len(epsilon)):
        ifabs(Phi1[:,k])<cutoff&abs(Phi2[:,k])<cutoff:
            energy=epsilon(k)
            break            

    return energy