Aproximação das derivadas de 1ª ordem

> restart;

> #====================================================

> DifFinPro:=proc(f::procedure,x0::numeric,h::positive)
local res;
res:=evalf((f(x0+h)-f(x0))/h);
RETURN(res);
end:

> #====================================================

> DifFinReg:=proc(f::procedure,x0::numeric,h::positive)
local res;
res:=evalf((f(x0)-f(x0-h))/h);
RETURN(res);
end:

> #====================================================

> DifFinCen:=proc(f::procedure,x0::numeric,h::positive)
local res;
res:=evalf((f(x0+h/2)-f(x0-h/2))/h);
RETURN(res);
end:

> #====================================================