odj

from numpy import *
from pylab import *
from scipy import integrate

def f(z,x):
    return array([z[1],2*z[1]-z[0]])
def funk(x):
    return (e**x)+x*e**x
y0=array([1.,2.])
xi=linspace(0,2,100)
yi=integrate.odeint(f,y0,xi)
print 'Aproksimativne vrijednosti yi:'
for i in range(0,100):
    print 'y%d=%.15f' % (i,yi[i,0])
print 'Vrijednost funkcije u tocki x=2:', funk(2)
print 'Greska za tocku x=2: ', abs(yi[99,0]-funk(2))

xt=linspace(0,2,10)
plot(xt,fun(xt),'co',xi,yi[:,0],'-b')
grid()
show()