вычисление объема параболойда по сегментам///need help

[Diller]

задача: необходимо разбить параболойд,заданной длины L, на части с определенным шагом (шаг задается пользователем) и посчитать обем каждого из элементов записав его в таблицу..
ниже приведена прога которая может служить примером.. но в ней есть ошибки.. я их найти не могу... (долго пытался)...

const pi=3.14;
var a,b,c,L,losk,x1,x2,D,xmax,ymin,ymax ,x,razr,rmax,r, xmin,
mindiam,maxdiam : real;
i,razmer:integer;
V: array of real;
begin
a:=strtofloat(Edit1.text);
b:=strtofloat(Edit2.text);
c:=strtofloat(Edit3.text);
losk:=strtofloat(Edit4.text);
L:=strtofloat(Edit5.text);

if a<0 then
begin
showmessage('Неправильно введено значение a. А -СТРОГО ПОЛОЖИТЕЛЬНОЕ ЧИСЛО');
end
else
begin
if abs(a)<0.000001 then
begin
if b=0 then
showmessage('При a=0 неправильно введено значение b. b - не равно "0"')
else
begin
razr:=round(L/losk);
edit6.Text:=floattostr(razr);
razmer:=strtoint(edit6.Text);
edit7.Text:=floattostr(razmer);
SetLength(V,razmer);
x:=10;
ymin:=b*(10)+c+losk;
edit8.Text:=floattostr(ymin);
ymax:=L+ymin;
xmax:=(ymax-c-losk)/b;
edit8.Text:=floattostr(xmax);
for i:=0 to razmer-1 do
begin
repeat
ymin:=b*x+c+losk;
ymax:=losk+ymin;
rmax:=(ymax-c)/b;
r:=(rmax+x)/2;
V[i]:=pi*r*r*losk;
ttable.cells[i+1,1]:=floattostr(V[i]);
x:=rmax;
until x>=xmax;
end;
end;
end
else
begin
razr:=round(L/losk);
edit6.Text:=floattostr(razr);
razmer:=strtoint(edit6.Text);
edit7.Text:=floattostr(razmer);
SetLength(V,razmer);
ymax:=L+a*(-10+(-b/(2*a)))*(-10+(-b/(2*a)))+b*(-10+(-b/(2*a)))+c+losk;
d:=b*b - 4*a*(c-ymax);
if d<0then
begin
showmessage('Неправильно введены a,b или c');
end
else
begin
x1:=(-b - sqrt(d))/(2*a);
{ x2:=(-b + sqrt(d))/(2*a); }
xmin:=x1;{до этого значения считаем в цикле}
mindiam:=20;{условие мин диаметра}
x:=-10+((-b)/(2*a));
for i:=0 to razmer-1 do
begin
repeat
ymin:=a*(x+(-b/(2*a)))*(x+(-b/(2*a)))+b*(x+(-b/(2*a)))+c+losk;
ymax:=ymin+losk;
d:=b*b - 4*a*(c-ymax);
if d<0then
begin
showmessage('Неправильно введены a,b или c');
end
else
begin
x1:=(-b - sqrt(d))/(2*a);
x2:=(-b + sqrt(d))/(2*a);
maxdiam:=x2-x1;
r:=((maxdiam/2)+(mindiam/2))/2 ;
V[i]:=pi*r*r*losk;
ttable.cells[i+1,2]:=floattostr(V[i]);
mindiam:=maxdiam;
x:=x1;
end;
until x>=xmin;
end;
end;
end;
end;
end;

[Diller]

losk- шаг
L - длина параболойда
a/b/c - параметры задающие форму парабулы вращением которой определяется параболойд

[Diller]

idies.doc вот идея того как я пытался решить эту проблему