ValueError: kesalahan domain matematika

Saya baru saja menguji contoh dari Metode Numerik dalam Teknik dengan 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)

Ketika saya menjalankannya, itu menunjukkan kesalahan berikut:

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

Saya telah mempersempitnya ke log seperti ketika saya menghapus log dan menambahkan fungsi yang berbeda, itu berhasil. Saya berasumsi itu karena semacam gangguan dengan pangkalan, saya tidak tahu caranya. Adakah yang bisa menyarankan solusi?

Kode Anda melakukan logangka yang kurang dari atau sama dengan nol. Itu tidak terdefinisi secara matematis, jadi logfungsi Python memunculkan pengecualian. Berikut contohnya:

>>> from math import log
>>> log(-1)
Traceback (most recent call last):
  File "<pyshell#59>", line 1, in <module>
    log(-1)
ValueError: math domain error

Tanpa mengetahui apa newtonRaphson2fungsi Anda, saya tidak yakin saya bisa menebak dari mana nilai yang tidak valid x[2]itu berasal, tapi mudah-mudahan ini akan membawa Anda ke jalur yang benar.

Anda mencoba melakukan logaritma dari sesuatu yang tidak positif.

Logaritma mencari tahu basis setelah diberi nomor dan kekuatannya dinaikkan. log(0)berarti bahwa sesuatu yang dipangkatkan 2adalah 0. Eksponen tidak pernah bisa menghasilkan 0*, yang berarti log(0)tidak ada jawaban, sehingga membuangmath domain error

* Catatan: 0^0bisa menghasilkan 0, tapi bisa juga menghasilkan 1pada saat bersamaan. Masalah ini sangat diperdebatkan.

Anda juga dapat menggunakan math.log1p.

Menurut dokumentasi resmi :

math.log1p (x)

Kembalikan logaritma natural dari 1 + x (basis e). Hasilnya dihitung dengan cara yang akurat untuk x mendekati nol.

Anda dapat mengonversi kembali ke nilai asli menggunakan math.expm1yang mengembalikan epangkat x, minus 1.

Anda mendapatkan kesalahan domain matematika karena salah satu alasannya: Anda mencoba menggunakan angka negatif di dalam fungsi log atau nilai nol.