Hola estoy haciendo un programa que mediante el método de Ritz nos de las frecuencias principales de una placa. Para ello uno de los pasos es la multiplicación de dos funciones de x, y su posterior integración el problema viene en que la integracion no se realiza de forma correcta y no se por qué.
Aqui dejo parte del codigo.
import numpy as np
import scipy.integrate as integrate
import matplotlib.pyplot as plt
def X(x,i):
X= np.sqrt(2)*np.sin(i*np.pi*x/a)
return X
def I00(x,i,m):
I00=X(x,i)*X(x,m)
return I00
integrate.quad(I00,0,a,args=((mm+1),(i+1)))[0]
Tanto i como m son variables que se van variando dentro de un bucle.