Вопрос о Python проекте о полиномах Лежандра и матрицах Грама

Question

Мой вопрос касается небольшого проекта Python (цифры, изучение математики).

Я относительно уверен, что матрицы Грамма и базы должны быть правильно алгоритмизованы. Он прекрасно работает с базой мономов, но выдает ошибку, которую я не понимаю, для других двух ...

Ошибка - объект модуля не вызывается: в функции prelambdalegendre для for для l oop

import math
from scipy import integrate as int
import numpy as np

def gaussklammer(n):
    if n%2==0:
        return n
    elif n%2==1:
        return n-1

monom= lambda x,n: x**n


def prelambdalegendre(x,k):
    pol=0
    for i in range (0,int(gaussklammer(k)/2)+1):
                    pol+= (-1)**i *math.factorial(2*k-2*i)/(math.factorial(k-i)*math.factorial(k-2*i)*math.factorial(i)*2**k)*(x**(k-2*i))

    return pol

legendre = lambda x,n:prelambdalegendre(x,n)


normlegendre =lambda x,k: math.factorial(2*k)/(2**k *math.factorial(k)**2) *legendre(x,k)

def grammatrix(baseofchoice,size):
    if baseofchoice=='monom':
        base =lambda x,k:monom(x,k)
    elif baseofchoice=='legendre':
        base =lambda x,k:legendre(x,k)
    elif baseofchoice=='normlegendre':
        base =lambda x,k:normlegendre(x,k)    
    #More elegant implementations didn't work , unfortunately.
    #To add another base, just add another elif statement

    A=np.zeros((size,size))


    for i in range (0,size):
        for j in range (0,size):
            f= lambda x : base(x,i)*base(x,j)
            A[i][j]=(int.quad(f,-1,1)[0])


    return A

print(grammatrix('monom',5))
print(grammatrix('legendre',5))

Answer 1

Нашел другой способ сделать это, спасибо за ваш вклад в любом случае;)

import math
from scipy import integrate
import numpy as np


monom= lambda x,n: x**n


def legendre(x,k):
    if k==0:
        return 1
    elif k==1: 
        return x
    else :

        pol=x*legendre(x,k-1)-((k-1)**2)/(4*(k-1)**2 -1)*legendre(x,k-2)

    return pol


normlegendre =lambda x,k: math.factorial(2*k)/(2**k *math.factorial(k)**2) *legendre(x,k)

def grammatrix(baseofchoice,size):
    if baseofchoice=='monom':
        base =lambda x,k:monom(x,k)
    elif baseofchoice=='legendre':
        base=lambda x,k:legendre(x,k)
    elif baseofchoice=='normlegendre':
        base=lambda x,k:normlegendre(x,k)    
    #Dieser Liste können nach Bedarf für neue Basen neue gleichförmige Clauses
    #hinzugefügt werden, elegantere Implementierungen scheiterten leider.

    A=np.zeros((size,size))


    for i in range (0,size):
        for j in range (0,size):
            f= lambda x : base(x,i)*base(x,j)
            A[i][j]=(integrate.quad(f,-1,1)[0])


    return A
A=grammatrix('monom',4)
print(A)
B=grammatrix('legendre',4)
print(B)
C=grammatrix('normlegendre',4)
print(C)
condA=np.linalg.cond(A)
condB=np.linalg.cond(B)
condC=np.linalg.cond(C)
print(condA)
print(condB)
print(condC)

Answer 2

Ваша проблема в том, что вы объявляете scipy.integrate как int,

и позже вызываете функцию int eger из python в вашей функции для l oop.

При импорте integrate попробуйте импортировать его под другим именем.