function img=sobel55(i);

i=imresize(i,[256 256]);
[r c]=size(i);
e=zeros(r+4,c+4);
for m=3:r+2
    for n=3:c+2
        e(m,n)=i(m-2,n-2);
        
    end
end
% figure, imshow(e,[])
gx=zeros(r,c);
gy=zeros(r,c);

for m=3:r+2
    for n=3:c+2
        gx(m-2,n-2)=(e(m-1,n-1)*1)+(e(m-1,n)*2)+(e(m-1,n+1)*1)+(e(m,n-1)*0)+(e(m,n)*0)+(e(m,n+1)*0)+(e(m+1,n-1)*-1)+(e(m+1,n)*-2)+(e(m+1,n+1)*-1)+(e(m-2,n-2)*2)+(e(m-2,n-1)*2)+(e(m-2,n)*4)+(e(m-2,n+1)*2)+(e(m-2,n+2)*2)+(e(m-1,n-2)*1)+(e(m-1,n+2)*1)+(e(m,n-2)*0)+(e(m,n+2)*0)+(e(m+1,n-2)*-1)+(e(m+1,n+2)*-1)+(e(m+2,n-2)*-2)+(e(m+2,n-1)*-2)+(e(m+2,n)*-4)+(e(m+2,n+1)*-2)+(e(m+2,n+2)*-2);
                    
        gy(m-2,n-2)=(e(m-1,n-1)*-1)+(e(m-1,n)*0)+(e(m-1,n+1)*1)+(e(m,n-1)*-2)+(e(m,n)*0)+(e(m,n+1)*2)+(e(m+1,n-1)*-1)+(e(m+1,n)*0)+(e(m+1,n+1)*1)+(e(m-2,n-2)*-2)+(e(m-2,n-1)*-1)+(e(m-2,n)*0)+(e(m-2,n+1)*1)+(e(m-2,n+2)*2)+(e(m-1,n-2)*-2)+(e(m-1,n+2)*2)+(e(m,n-2)*-4)+(e(m,n+2)*4)+(e(m+1,n-2)*-2)+(e(m+1,n+2)*2)+(e(m+2,n-2)*-2)+(e(m+2,n-1)*-1)+(e(m+2,n)*0)+(e(m+2,n+1)*1)+(e(m+2,n+2)*2);
        end
end
img=double((gx.^2 + gy.^2).^(0.5));