Написать программу (паскаль) для нахождения корней квадратного уравнения и ответ

x²-x-6=0

1. Вот так это пишется с функцией

// PascalABC.NET 3.1, сборка 1200 от 13.03.2016
function IsPrime(n:integer):boolean;
begin
  if n<4 then Result:=True
  else begin
    var found:= (n mod 2 = 0);
    var p:=3;
    while (not found) and (sqr(p)<=n) do
      begin
      found:=(n mod p = 0);
      p+=2
      end;
    Result:=not found
    end
end;

begin
  Writeln('k=',ArrRandom(ReadInteger('n='),1,999).Println.
    Where(x->IsPrime(x)).Count)
end.

Тестовое решение:
n= 10
401 828 780 444 694 965 23 341 673 875
k=3

2. А вот так это пишется с процедурой

// PascalABC.NET 3.1, сборка 1200 от 13.03.2016
procedure IsPrime(n:integer; var res:boolean);
begin
  if n<4 then res:=True
  else begin
    var found:= (n mod 2 = 0);
    var p:=3;
    while (not found) and (sqr(p)<=n) do
      begin
      found:=(n mod p = 0);
      p+=2
      end;
    res:=not found
    end
end;

begin
  var a:=ArrRandom(ReadInteger('n='),1,999); a.Println;
  var k:=0;
  var prime:boolean;
  foreach var e in a do begin
    IsPrime(e,prime);
    if Prime then Inc(k)
    end;
  Writeln('k=',k)
end.

Тестовое решение:
n= 12
199 43 71 365 417 904 170 212 694 103 161 689
k=4

# Задача 1

def calc(num: int):

   res = 1

   for i in range(1, num + 1):

       res += res * i

   return res

user_num = int(input())

res = 0

for i in range(user_num):

   res += calc(i)

print(f'ответ: {res}')

# Задача 2

num = int(input())

res = []

while num > 0:

   res.append(str(num % 2))

   num //= 2

print(''.join(res[::-1]))

# bin(num) - встроенная функция в Python для перевода в двоичную систему

# Задача 3

command = input('Введите номер фигуры:\n1 - Прямоугольник\n2 - Круг\n3- Треугольник')

if command == '1':

   side_a = int(input('Введите длину первой стороны: '))

   side_b = int(input('Введите длину второй стороны: '))

   print(f'S = {side_a} * {side_b} = {side_a * side_b}')

   # Ещё есть формула через диагонали и синус угла, можно это реализовать, не сильно меняет

elif command == '2':

   rad = int(input('Введите радиус окружности: '))

   print(f'S = 3.14 * {rad}^2 = {3.14 * (rad ** 2)}')

   # Есть ещё формулы через диаметр и длину окружности

elif command == '3':

   h = int(input('Введите длину высоты: '))

   side = int(input('Введите длину стороны, к кторой проведена высота: '))

   print(f'S = {h} * {side} = {h * side}')

   # А у треугольника этих формул хренова туча