creating mach band

% creating mach band
clear all; close all; clc;
for i=1:256
for j=1:256
if j >224
f(i,j)= 224;
elseif j >192
f(i,j)=192;
elseif j > 160
f(i,j)=160;
elseif j >128
f(i,j)=128;
elseif j>96
f(i,j)=96;
elseif j>64
f(i,j)=64;
elseif j>32
f(i,j)=32;
else
f(i,j)=5;
end

end
end
f = uint8(f);
figure,imshow(f),title('Mach Band');
v = f(1,1:256);
figure,plot(v),ylabel('Intensity'),xlabel('Pixel Position');