Переписать алгоритм в код

# ЧМК
from sympy import *
 
def CheckProstoeChislo(n): #проверка на простое число
    lst = []
    # в k будем хранить количество делителей
    k = 0
    # пробегаем все числа от 2 до N
    for i in range(2, n + 1):
        # пробегаем все числа от 2 до текущего
        for j in range(2, i):
            # ищем количество делителей
            if i % j == 0:
                k = k + 1
        # если делителей нет, добавляем число в список
        if k == 0:
            if n == i:
                return
            lst.append(i)
        else:
            k = 0
    # выводим на экран список
    print('=============================')
    print(f'Вы ввели не простое число.\n*[ПОДСКАЗКА]* Список простых числе до {n} {lst}\nВыбирите число.')
    print('=============================')
    n = int(input('\tВведите простое число p='))
    CheckProstoeChislo(n) #
 
def CheckPNaturalChislo(a):
    if a <= 1:
        print('=============================')
        print(f'Вы ввели число не больше 1.\n*[ПОДСКАЗКА]* Введите простое число больше 1')
        print('=============================')
        a = int(input('\tВведите простое число a='))
        CheckPNaturalChislo(a)  #
 
def korni():
    e = [0, 0]
    for i in range(0, 100):
        if (f[3] * i ** 3 + f[2] * i ** 2 + f[1] * i + f[0]) % p == 0:
            if e[0] == 0:
                e[0] = i
            elif e[1] == 0:
                e[1] = i
    return e
 
def proizvodnaia():
 
    x = Symbol('x')
    y = int(f[3]) * x ** 3 + int(f[2]) * x ** 2 + int(f[1]) * x + int(f[0])
 
    proizvodnaia = y.diff(x)
 
    return proizvodnaia
 
def proizvodanaiaE(proizvodnaia,e,p):
    x = Symbol('x')
    A1 = (proizvodnaia.subs(x, e[0]))
    A11 = A1 % p
    A2 = (proizvodnaia.subs(x, e[1]))
    A22 = A2 % p
 
    print(f'f({e[0]}) = {A1} mod {p} = {A11}\t\tf({e[1]}) = {A2} mod {p} = {A22}')
 
if __name__ == '__main__':
    p = int(input('\tВведите простое число p='))
    CheckProstoeChislo(p)
 
    a = int(input('\tВведите натуральное число a='))
    CheckPNaturalChislo(a)
 
    f = [0, 0, 0, 0, 0]
    f[3] = int(input('\tВведите число f='))
    f[2] = int(input(f'\tВведите число f={f[3]}x^3 + '))
    f[1] = int(input(f'\tВведите число f={f[3]}x^3 + {f[2]}x^2 + '))
    f[0] = int(input(f'\tВведите число f={f[3]}x^3 + {f[2]}x^2 + {f[1]}x + '))
    print(f'Результат f(x)={f[3]}x^3 + {f[2]}x^2 + {f[1]}x + {f[0]}')
    print('=============================')
 
    e = korni()
    print(f'Корни {e}')
 
    proizvodnaia = proizvodnaia()
    print(f"Производная (** это ^ степень) f'={proizvodnaia}")
 
    proizvodanaiaE(proizvodnaia,e,p)
 
# вычислить многочлен