def sink(i): if i%2==0: return 0 return rat(1,factorial(i)*(-1)**((i-1)/2)) a = [ sink(i) for i in range(n) ] out = open("sin-funcroot.py","w") out.write("from rational import rat\n") out.write("a=[]\n") try: idx = 1 for r in funcrootgen(a, 1): out.write("a.append(%s)\n"%repr(r)) out.write("\n") out.flush() print idx, float(r) idx+=1 except Exception: out.close()
def sink(i): if i%2==0: return 0 return rat(1,factorial(i)*(-1)**((i-1)/2))