function div_dif(x0, x1)

div_dif =(f(x1) -f(x0))/(x1-x0)
end


function div_dif_2(x0, x1, x2)

div_dif_2 =( div_dif(x2,x1) -div_dif(x1,x0) )/(x2-x0)
end
