Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import math
- def func(x,y):
- a = math.cos(x)
- b = math.sin(y - 1)
- znam = a * b + 4
- Matr = [[b / znam, 2 / znam], [2 / znam, a / (-znam)]]
- F = [float(math.sin(x) + 2 * y - 2), 2 * x + math.cos(y) - 0.7]
- X1 = [0,0]
- Xk = [x, y]
- for i in range(0,2):
- helper = Matr[i][0] * F[0] + Matr[i][1] * F[1]
- X1[i] = Xk[i] - helper
- return X1
- x = float (input())
- y = float (input())
- while(True):
- Xk1 = func(x,y)
- print(Xk1)
- if math.sqrt(abs(Xk1[0]**2 - x**2) + abs(Xk1[1]**2 - y**2)) < 0.001:
- break
- else:
- x = Xk1[0]
- y = Xk1[1]
- print(x,y)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement