import matplotlib.pyplot as plt
import numpy as np


minMap = -np.pi*0
maxMap = np.pi*1
stepMap = 0.1
psi = 0.1
e = 10**(-3)

def f (x, y = np.pi):
	return np.cos(x) + np.sin(x) - y**3 - psi

def f1 (x, y):
	return (f(x, y + e) - f(x, y)) / e

def newtons_method(x, y0, f, f1):
    #f1 - производная
    y0 = float(y0)
    k = 0
    while True:
        y1 = y0 - (f(x, y0) / f1(x, y0))
        if abs(y1 - y0) < e:
        	print(y1)
        	return y1
        y0 = y1
        k += 1

def func(t):
	N = len(t)
	y = 0
	i = 0
	while i < N:
		y = newtons_method(t[i], 0.01, f, f1)
	return y

# Data for plotting
t = np.arange(minMap, maxMap, stepMap)
s = func(t)

fig, ax = plt.subplots()
ax.plot(t, s)

ax.set(xlabel='time (s)', 
	ylabel='voltage (mV)',
	title='About as simple as it gets, folks')
ax.grid()

plt.show()

# Data for plotting
s =[]
t = np.arange(minMap, maxMap, stepMap)
for val in t:
	s.append(func(val))