#python #numpy #integration
#python #numpy #интеграция
Вопрос:
Я аппроксимирую интеграл, используя квадрат Гаусса-Лагерра, используя функцию numpy.polynomial.laguerre.laggauss в python. Мне нужно использовать большое количество точек, чтобы получить точные результаты, но, похоже, я не могу ввести значение deg> 185 точек, и у меня возникает следующая ошибка.
Файл «/usr/local/lib/python2.7/dist-packages/numpy/polynomial/laguerre.py «, строка 1726, в laggauss w = 1 /(fm * df) Предупреждение о времени выполнения: переполнение, обнаруженное в true_divide
Есть ли способ преодолеть это в python? Спасибо.
Ответ №1:
Кажется, теперь доступно:
import numpy as np
n = 186
(root, omega) = np.polynomial.laguerre.laggauss(n)