function P1=NDDinterpol (x0,x1,f0,f1)
    x=poly(0,"x");
    f01=(f1-f0)/(x1-x0);
    P1=f0+(x-x0)*f01;
endfunction