Суммирование

program pr;
uses crt;
const n=5;
type mas=array[1..n+1] of real;
 
function Pow(os:real;pk:byte):extended;
var i:byte;
    res:real;
begin
res:=1;
for i:=1 to pk do
res:=res*os;
Pow:=res;
end;
 
function Summ(x:real;m:mas;k:byte):extended;
var i:byte;
    sm:real;
begin
sm:=0;
for i:=n+1 downto 1 do
sm:=sm+Pow(x,k)*m[i];
Summ:=sm;
end;
 
procedure Proiz(m1:mas;var m2:mas);
var i:byte;
begin
for i:=1 to n do
m2[i]:=m1[i]*(n-i+1);
end;
 
procedure Vyvod(x:real;var m:mas;k:byte);
var i:byte;
begin
write(m[1]:0:2,'x^',k);
for i:=2 to k-1 do
if m[i]<0 then write(m[i]:0:2,'x^',k-i+1)
else write('+',m[i]:0:2,'x^',k-i+1);
if m[n]<0 then write(m[k]:0:2,'x')
else write('+',m[k]:0:2,'x');
if m[k+1]<0 then write(m[k+1]:0:2,'=')
else write('+',m[k+1]:0:2,'=');
write(Summ(x,m,i-1):0:8);
writeln;
end;
 
var a,p:mas;
    x:real;
    i:byte;
 
begin
clrscr;
randomize;
for i:=1 to n+1 do
begin
   writeln('Vvedite ',i,' koefficient polinoma');
   readln(a[i]);
end;
 
write('x=');
readln(x);
 
Proiz(a,p);
writeln;
writeln('Proizvodnaya:');
writeln;
Vyvod(x,p,n-1);
 
readln
end.