num_angles = 400;
det_length = 1.5;
num_dets = 500;
res = 0.004;
filtertype = 1; % use '1' for convolution based filter, any other number
		% for FFT based filter


[angles,projs,l,det_space] = calcproj(num_angles,det_length,num_dets);
fprintf('Projections Generated\n');
[bpimage,fprojs] = backproj(projs,angles,det_space,l,res,filtertype);

rimage = flipud(bpimage');
load orig
error = orig - rimage;
sqerr = sum(sum(error.^2))


%save rimage rimage -ascii

%[orig] = projgraf(0,0,0,0.004);
%save orig orig -ascii
%quit
