% Adds vertical blur to image I
% using the D from makeD.m.
% I must be defined in the workspace.

J = D*I;

figure('Units', 'Normalized', 'Position', [.05, .2, .4, .6])

subplot(2,1,1)
imagesc(I)
colormap(gray)
axis image
colorbar
title('Original Image I')

subplot(2,1,2)
imagesc(J)
colormap(gray)
axis image
colorbar
title('D-Blurred Image J')