Función y derivada numérica (método Euler adelante)

% AUTOR:
% FECHA:
ndata=4096;
x=linspace(-10,10,ndata);
fun=x.^3+x+2;
fun1=fun+rand(1,ndata);
hold on;
plot(x,fun,'k-');
plot(x,fun1,'r+');
%graf2=figure(2);
%graf2=plot(x,fun-fun1,'r+'); % gráfica de la derivada de la función
%afectada por el random
%DERIVADA NUMERICA
deltax= x(2)-x(1);
for k=1:(ndata-1);
deri(k)= ((fun(k+1)-fun(k))/deltax);
end;
graf3=figure(3);
graf3=plot(x(1:ndata-1),deri,'b-','LineWidth',1.5);
pos=[30 100 550 350]; %determino la posición de la gráfica
set(figure(3),'Position',pos); % constituye la posición elegida