Алгоритм не работает с n> 5, в чем проблема?

Question

Я пытаюсь решить математическую задачу.Мой алгоритм работает с n <5 и застрял с n = 6, 7 ... Возможно, я выбрал неправильный решатель, или, может быть, слишком много циклов.Где ошибка?Помогите, пожалуйста. </p>

import numpy as np
import sympy
n=int(input("Number of it.: "))
nu=sympy.Symbol('nu')
y=sympy.Symbol('y')
dy=np.linspace(0.0,1.0,n)
F=sympy.symbols('b0:%d' %n)
Eq15=1
psi=1

for i in range (n):
    Eq15=(Eq15+F[i])
    psi=(psi+F[i]*y**(i+2))
equations=[sympy.Eq(Eq15,0)]
Eq12=psi.diff(y,y)+nu*psi*(1-y**2) 

for i in range (n-1):
    Eq12_=Eq12.subs(y,dy[i])
    equations.append(sympy.Eq(Eq12_,0))
sol=sympy.solve(equations,sympy.symbols('b0:%d'%n))

psi_final=1
for key, value in sol.items(): 
   psi_final=psi_final+sol[key]*(y**((list(sol).index(key))+2))
print(psi_final)

V120=psi_final.diff(y,y)+nu*psi_final*(1-y**2)
Z=(sympy.integrate(V120,(y,0,1)))
Z1=sympy.solve(Z,nu)

Я пытаюсь решить эту проблему, вот пример для 4 значений y: imgur.com/FU9i7wA