load sigma6;
load sigma6blur;
pic=output;
%pic=blur3;

%img = double(imread('../peppers512x512.tif', 'tif'));
%img = img(4:4:512,4:4:512);

%load test2
%N = 1;
% subsample by using only every fourth element
%pic = test2(N:N:rows,N:N:cols);

% subsample by using only every fourth element
rows = size(pic,1);
cols = size(pic,2);

% Do the solver
U = pde_solve(pic);

disp('Displaying picture');
imagesc(uint8(full(U)));