Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- # -*- coding: utf-8 -*-
- """
- Spyder Editor
- This is a temporary script file.
- TODO: Make plots
- """
- import math
- def mclrn(x, prec, imax):
- est = 0
- errs = []
- ests = []
- for i in range(imax):
- est += (pow(x, i) / math.factorial(i))
- err = abs(math.exp(x) - est) / math.exp(x)
- errs.append(err)
- ests.append(est)
- if err < prec:
- break
- # This block stops calculations if error doesn't decrease
- if len(errs) > 2:
- if errs[len(errs) - 1] == errs[len(errs) - 2]:
- break
- return ests, errs, i
- x1, y1, term1 = mclrn(1, 1e-15, 30)
- x2, y2, term2 = mclrn(2, pow(10, -15), 300)
- x10, y10, term10 = mclrn(10, pow(10, -15), 300)
- for i in range(len(x1)):
- print("{}: {}".format(x1[i], y1[i]))
- print(term1)
- for i in range(len(x2)):
- print("{}: {}".format(x2[i], y2[i]))
- print(term2)
- for i in range(len(x10)):
- print("{}: {}".format(x10[i], y10[i]))
- print(term10)
- '''
- Report. Student: Alisher Shakhiyev
- In this assignment
- '''
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement