parte3.sce

// Respuesta del punto 3.
exec('data.sci', -1); // Carga de la data.
// Cálculo de la derivada hacia adelante
// Parámetros:
//      f           Función a evaluar
//      h           Diferencia entre los valores
//      xsDelante   Valor de x "adelantado"
//      xsCentrada  Valor de x de referencia
function [y]= DDelante(f, h, xsDelante, xsCentrada)
    Y = (f(xsDelante) -  f(xsCentrada))/h;
endfunction
// Cálculo de la derivada hacia atrás
// Parámetros:
//      f           Función a evaluar
//      h           Diferencia entre los valores
//      xsAtras     Valor de x "retrasado"
//      xsCentrada  Valor de x de referencia
function [y]= DAtras(f, h, xsAtras, xsCentrada)
    Y = (f(xsCentrada) -  f(xsAtras))/h;
endfunction
xsDelante= x + hs; // xs delanteras
xsDetras= x - hs;
DYAtras= zeros(hs);
DYDelante= zeros(hs);
for i = [1 : size(hs)(2)]
    DYDelante(i)= DDelante(f, hs(i), xsDelante(i), x);
    DYAtras(i)= DAtras(f, hs(i), xsDetras(i), x);
end