Курсач. Площадь поверхности усеченного цилиндра.

uses crt;
function V(r1,r2,h:real):real;
begin
V:=pi*h*(sqr(r1)+r1*r2+sqr(r2))/3;
end;
function H(v,r1,r2:real):real;
begin
H:=3*v/(pi*(sqr(r1)+r1*r2+sqr(r2)));
end;
var h1,r1,r2,g,n,m,a,s,h2,y:real;
begin
clrscr;
repeat
writeln('Введите радиусы усечённого конуса, r2>r1:');
write('r1=');readln(r1);
write('r2=');readln(r2);
until (r1>0)and(r2>r1);
repeat
write('Введите длину стороны куба: ');
readln(a);
until a>0;
repeat
writeln('Введите интервал значений изменения высоты n>g:');
write('Начальное значение g<',H(a*a*a,r1,r2):0:2,' g=');readln(g);
write('Конечное значение n>',H(a*a*a,r1,r2):0:2,'=');readln(n);
until (g>=0)and(g<H(a*a*a,r1,r2))and(n>H(a*a*a,r1,r2));
repeat
write('Введите размер шага m<',H(V(r1,r2,n),r1,r2)-H(a*a*a,r1,r2):0:2,' m=');
readln(m);
until (m>0)and(m<H(V(r1,r2,n),r1,r2)-H(a*a*a,r1,r2));
h1:=g;
repeat
y:=V(r1,r2,h1);
h1:=h1+m;
until y>=a*a*a;
writeln('Высота усечённого конуса должна быть равна ',h1:0:2);
h2:=3*y/(pi*sqr(r2));
s:=sqrt(sqr(h2)+sqr(r2));
writeln('Высота конуса =',h2:0:2);
writeln('Его образующая=',s:0:2);
readln
End.