Я просто тестировал пример из Numericical Methods in Engineering с Python.
from numpy import zeros, array
from math import sin, log
from newtonRaphson2 import *
def f(x):
f = zeros(len(x))
f[0] = sin(x[0]) + x[1]**2 + log(x[2]) - 7.0
f[1] = 3.0*x[0] + 2.0**x[1] - x[2]**3 + 1.0
f[2] = x[0] + x[1] + x[2] -5.0
return f
x = array([1.0, 1.0, 1.0])
print newtonRaphson2(f,x)
Когда я запустил его, появится следующая ошибка:
File "example NR2method.py", line 8, in f
f[0] = sin(x[0]) + x[1]**2 + log(x[2]) - 7.0
ValueError: math domain error
Я сузил его до журнала, как при удалении журнала и добавлении другой функции, он работает. Я предполагаю, что это связано с каким-то вмешательством в базу, я не могу понять, как это сделать. Может ли кто-нибудь предложить решение?