Évaluation d'un polynôme

La function evalpoly permet l'évaluation d'un polynome en un point.

      
function evalpoly(P:polynome; x : real):real;
    var i : integer;
        somme : real;

    function puissance(x : real; n : integer):real;
      begin
        if (n=0) then puissance:=1
        else puissance:=x*puissance(x,n-1);
    end;

    begin
       somme:=P.coef[0];
       for i:=1 to P.degre do
           somme:=somme+P.coef[i]*puissance(x,i);
       evalpoly:=somme;
end;