Нахождение дирекционного угла

from turtle import *
from math import sin,cos,tan,atan,degrees,pi
from random import randint
 
 
R = 100 # радиус окружности построения кругов
r = 30 # радиус кругов в окружности
n = 4 # количество строк матрицы
Fi = 2*pi/n # угол смещения в радианах [1.5707963267948966]
p = [] # координаты кругов по оси x & y
matrix = [[randint(0,9) for _ in range(n)] for _ in range(n)] # генерация случайной двумерной матрицы
 
speed(10)
 
 
for i in range(n):
    """Вычисление координат x & y"""
    p += [[(R * cos(Fi * i)), (R * sin(Fi * i))]] 
 
 
# print(matrix)
# print(p)
 
for i in range(n):
    """Отрисовка кругов и их нумерация"""
    penup()
    goto(p[i % n][0], p[i % n][1] - r)
    pendown()
    circle(r)
    penup()
    goto(p[i % n][0], p[i % n][1] + r)
    pendown()
 
    write(i+1)
 
speed(1)
 
for i in range(n):
    for j in range(n):
        if matrix[i][0] != 0 and matrix[i][j] != 0 and i != j:
            """Нахождение зависимостей элементов матрицы A & B"""
            penup()
            setpos(p[i][0], p[i][1]) # круг A
            pendown()
            goto(p[j][0], p[j][1]) # путь к B
 
            angle = tan((p[i][1] + 1 - p[j][1]) / (p[i][0] + 1 - p[j][0])) # вычисление угла наклона
            angle = atan(angle)
            angle = degrees(angle) # приведение из меридианы в градусы
 
            print(f"A: {i+1}, B: {j+1}, Градус черепашки: {heading()}, Дирекциональнный угол: {angle}")
            
            setheading(angle) # установка черепашки в заданом градусе
            stamp() # отрисовка стрелки
done()