Matlab - Derivada de un a imagen

function der = derivar_imagenes(imagen)


imagen=double(imagen);
[M N]=size(imagen);
for i=1:M;
der1(i,:)=diff(imagen(i,:));%funciion de derivar una señal, la aplicamos a cada una de las filas
end;
imagen_T=imagen'*(-1);
for i=1:N;
der2(i,:)=diff(imagen_T(i,:));% se transpone la matriz de la imagen y se deriva cada fila
end;

der3=der2';
der1=der1(1:M-1,:);
der3=der3(:,1:N-1);
der=sqrt(der1.^2+der3.^2);