Найти внутренние углы треугольника (Pascal)

[street_spirit]

Здравствуйте! В программировании я совсем новичок и столкнулся с проблемой.

Перед мной стоит следующая задача: нужно написать программу, которая высчитывает внутренние углы треугольника (пользователь задает координаты вершин треугольника).

Я придумал следующий алгоритм решения: берем координаты вершин и считаем координаты векторов (сторон треугольника) и потом по формуле угла между векторами считаем косинус угла. Узнав значение косинуса, пользуемся арккосинусом и получаем готовый результат (так проделываем 3 раза для разных углов).

Написав программу только до одного угла, я понял что она работает неправильно. В частности выдает значение косинуса больше единицы.

Код программы:

Program prog1;
Var x1, x2, x3, y1, y2, y3, xAB, yAB, xAC, yAC, cosA, arccosA: real;
Begin
Writeln('Vvedite koordinati tochki A (cherez probel)');
Readln(x1,y1);
Writeln('Vvedite koordinati tochki B (cherez probel)');
Readln(x2,y2);
Writeln('Vvedite koordinati tochki C (cherez probel)');
Readln(x3,y3);
Writeln('A(',x1,';',y1,')');
Writeln('B(',x2,';',y2,')');
Writeln('C(',x3,';',y3,')');
xAB:=x2-x1;
yAB:=y2-y1;
xAC:=x3-x1;
yAC:=y3-y1;
Writeln('AB{',xAB,';',yAB,'}');
Writeln('AC{',xAC,';',yAC,'}');
cosA:=(xAB*xAC+yAB*yAC)/((sqrt((xAB*xAB+yAB*yAB)))*(sqrt((xAC*xAC+yAC*yAC))));
Writeln('cosA=',cosA);
xAB:=x2-x1;
yAB:=y2-y1;
xAC:=x3-x1;
yAC:=y3-y1;
Writeln('AB{',xAB,';',yAB,'}');
Writeln('AC{',xAC,';',yAC,'}');
cosA:=(xAB*xAC+yAB*yAC)/((sqrt((xAB*xAB+yAB*yAB)))*(sqrt((xAC*xAC+yAC*yAC))));
Writeln('cosA=',cosA);
arccosA:= ARCTAN(SQRT(1-SQR(cosA))/cosA); 
Writeln('arccosA:=',arccosA);
end.

При тестировании с координатами А(3,2) В(7,7) С(9,3) выходит вот что:


Для наглядности (на всякий случай):


Помогите пожалуйста разобраться в чем моя ошибка. Изначально алгоритм решения задачи неправилен? Ошибка в коде? И еще мне кажется, что я что-то напутал с типом данных для переменных.

Заранее большое спасибо!

[eoln]

Что такое экспоненциальное представление
Writeln('cosA=',cosA:0:10);

[Rusat]

Для нахождения строн треугольника нужно пользоваться следующей формулой: s = ((x2 - x1)^2+(y2-y1)^2)^(1/2), а уже потом можно найти углы

[Namolem]

проверь скобки, похоже умножается на последний корень

[street_spirit]

Сообщение от eoln

Что такое экспоненциальное представление
Writeln('cosA=',cosA:0:10);

Спасибо! Очень помогло.