mnk

import numpy as np
from pylab import *

x = array ([3.,6.,9.,12.,15.,18.])
y = array ([173.,245.,300.,346.,387.,424.])

lny=log(y)

koef = numpy.polyfit(x,lny,1)
print(koef)

RegPravac = numpy.poly1d(koef)
print(RegPravac)

a = exp(RegPravac[0])
b = RegPravac[1]

print('a =', a)
print('b =', b)

def f(x):
    return a*x**b

print('Jedn. krivulje y = %0.f*e^(%0.f x) ' %(a,b))

print('Vrijednost funkcije u x = 36', f(36))

xtocke =linspace(min(x), max(x),50)
yf = f(xtocke)
plot(x,y,'ro', xtocke,yf,'g-')
grid()
show()