По заданным вещественным числам a(0),a(1) ,...,a(20), T вычислить значения многочлена a(20)*x^20+a(19)*x^19+...+a(1)*x+a(0) и значение его производной

{$N+}//для типа extended, если вводить большие Х
uses crt;
const n=20;//порядок полинома
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
a[i]:=5*random;//задаем коэффициенты, можно любые, можно ввод с экрана
write('x=');readln(x);
writeln('Summa polynoma:');
Vyvod(x,a,n);
Proiz(a,p);
writeln('Proizvodnaya:');
Vyvod(x,p,n-1);
readln
end.