Составить программу для решения этой формулы (Паскаль)

Программа (без проверки исходных данных):
var a,b,y,z:real;
begin
write(a, b: );
readln(a,b);
y:=cos(sqrt(a+b));
z:=ln(sin(y)/cos(y));
writeln(ln(tg(cos(sqrt(a+b)))) = ,z);
end.

Пример:
a, b: 0.5 0.8
ln(tg(cos(sqrt(a+b)))) = -0.813058267838025