Найти площадь прямоугольника по стороне и диагонали. в паскале

Program e;
var 
s,a,b,d : real;
begin
a:= 1; {известная сторона}
d:= 3; {диагональ}
{по Т. Пифагора находим не известную сторону}
b: = sqrt(d*d - a*a);

{находим площадь }

s:= a*b;
writeln(s);
end;