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

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;

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