cd /net/hi2/cpl/classes/cs7321_98_winter/images
rod = imread('rod.jpg');  
imshow(rod);
crod = improfile;
figure(2)
subplot(311)
plot(crod)
subplot(312)
plot(diff(crod))
subplot(313)
plot(diff(diff(crod)))