Matlab Code For Power Law Transformation

%Power law

a=imread('cameraman.tif');
subplot(1,2,1);
imshow(a);
[m,n]=size(a);
for i=1:m-1
for j=1:n-1
s(i,j)=a(i,j).^4;
end;
end;
subplot(1,2,2);
imshow(a);