Не работает медианный фильтр

%Возьмем некоторое исходное изображение 
L=imread('lili.bmp');
figure, imshow(L);
% Наложим на него импульсный шум с некоторыми характеристиками:
L=imnoise(L,'gaussian');
figure, imshow(L);
% Далее для каждой цветовой составляющей применим операцию медианной фильтрации.
for k=1; % обработка отдельно по каждой составляющей
   L=Lin(:,:,k);
     for i=1+n1:N+n1;     
         disp(i)  
         for j=1+m1:M+m1;
                if j==1+m1;
                           D=0;
                            for a=-n1:n1;   
                                for b=-m1:m1;
                                 D(n1+1+a,m1+1+b)=L(i+a,j+b);
                            end;                   end;
               end;
                     if j>1+m1;
               for a=-n1:n1;     
                   D(n1+1+a,m+1)=L(i+a,j+m1);    
               end;
                            D=D(1:n,2:m+1);
                     end;
              Lres(i,j)=median(D(:));
      end;                                  
     end;
end;

Error in Untitled (line 5)
L = imnoise(I,'salt &
pepper',0.02);
 
Error in run (line 63)
evalin('caller', [script ';']);