import numpy as np
import scipy as sci
from matplotlib import pyplot as pl
from scipy.special import factorial


def rozwiniecie(n, x):
  znak = np.ones(n+1)*(-1)
  znakn = np.linspace(1, n+1, n+1)
  znak = np.power(znak, znakn)
  print(znak)
  licznik = np.power(x, np.linspace(1, 2*n+1, n+1))
  mianownik = np.linspace(1, 2*n+1, n+1)
  print(mianownik)
  mianownik = factorial(mianownik) #silnia
  znak = np.multiply(znak,licznik)
  return np.sum(np.divide(licznik,mianownik))


test = rozwiniecie(10, 0.3)
bladw = 0.295520-test
bladbw = (0.295520-test)/0.295520

def tabela(analytic,approx,a,x):
    print("{}^{}".format(a,x))
    for item in ["n","analityczna","przyblizona","blad bezwzgledny","blad wzgledny (%)"]:
        print("{:<30}".format(item),end='')
    print()

def tab()
    for n in range(0,10):
    var = rozwiniecie(n, 0.5)
    bladw = 0.295520-var
    bladbw = (0.295520-var)/0.295520
    print("n", )



print("blad wzgl"+str(bladw))
print("blad bezwzgl"+str(bladbw*100))
# test = factorial(3)
print(test)