import numpy as np

def dy_dx(x): # definisanje diferencijalne jednačine  
    return 2*x  

x0 = 1 # početni uslovi  
y0 = 1 # početni uslovi  

xf = 1.5 # krajnja vrednost x, vrednost za koju zelimo da izracunamo vrednost jednačine  
n = 101 # broj tačaka na koje se deli interval od "x0" do "xf"  
m = n - 1 # broj podintervala (koraka) između "x0" i "xf"  

x = np.linspace(x0,xf,n) # interval od "x0" do "xf" deli na "n" tačaka  
y = np.zeros(n) # pravi skup od "n" clanova, pri cemu je svaki clan nula  
	  
h = (xf-x0)/m # ovako se definiše "h" na osnovu "np.linspace" komande  
              # a može i da se definiše "h" kao "h=x[1]-x[0]"   

# aproksimativno rešenje, Ojlerov metod  
for i in range (0,n-1):  
    y[0] = x0  
    x[0] = y0  
    y[i+1]= y[i] + h * dy_dx(x[i])

print("Vrednost koraka h je", h) 
print("Aproksimativna vrednost u tački", xf, "iznosi", y[i+1])