    function [p,pprim]=horn1(a,x)
    % Horners schema för beräkning av 
    % p=a(1) x^n + a(2) x^(n-1) + ... + a(n+1)
    n=length(a)-1;
    p=a(1);
    pprim=a(1);
    for i=1:n-1
      p=p*x+a(i+1);
      pprim=pprim*x+p;
    end
    p=p*x+a(n+1);

