Pasted

# Bruk import math heller enn astropy og masse piss, it has all you need
import math
 
def f(x):
    return x**3 - 3.0
 
def df(x):
    return 3*x**2 
 
 
def newtonRaphson(f, df, x, nmax, tol):
    #x0 er startverdi, nmax er max antall iterasjoner
    #x = x0 - (f(x0)/df(x0))
    
    for i in range(nmax):
        
        x = x - f(x)/(df(x))
        
        if (f(x)==0):
            return x
        elif (df(x)==0):
            break
        
    return x
    
x_start = 2.5
max_iter = 3
eps = 1.E-1
eps2 = 9.E-6
 
root = newtonRaphson(f, df, x_start, max_iter, eps)
 
 
def g(x):
    # Igjen -> math.log(variable, base), se dokumentasjonen
    return 5*x + math.log(x, math.e) - 10000
 
def dg(x):
    return 5 + (1/x)
 
 
print(root)
print(newtonRaphson(g, dg, 1500 , 10, eps2))
 
# Wehuu it works