figure(1)
clf
a = imread('Mammo CC P1.tif',  'tif');
b = a(:,:,2);
imagesc(a(:,:,2))
colormap(gray)
print -djpeg Lab1_2_1

figure(2)
clf
theta = 0:180;
R = radon(b,theta);
imagesc(R)
colormap(gray)
print -djpeg Lab1_2_2

figure(3)
clf
I2 = iradon(R,theta);
imagesc(I2)
colormap(gray)
print -djpeg Lab1_2_3

I = zeros(100,100);
I(25:75,25:75) = 1;
figure(4)
clf
imagesc(I)
colormap(gray)
print -djpeg Lab1_2_4

figure(5)
theta = 0:180;
R = radon(I,theta);
imagesc(R)
colormap(gray)
print -djpeg Lab1_2_5

figure(6)
I2 = iradon(R,theta);
imagesc(I2)
colormap(gray)
print -djpeg Lab1_2_6

