Replace ex3.m

ex3.m

@@ -28,7 +28,7 @@ axis_length = 16;  % size of region of interest
 b = 3*epsilon;  % band-width of the ramp filter in filtered backprojection
 
 RadonDeltaS = pi/b;
-RadonDeltaOmega = sqrt(2)/axis_length * pi/b;
+RadonDeltaOmega = 360 / (ceil(360 / (sqrt(2)/axis_length*pi/b)));
 
 plot_flag = 0;
 
@@ -101,7 +101,7 @@ print plots/ex3_geom -depsc
 RadonRotAngles = (0:RadonDeltaOmega:360);
 RadonRotAngles = RadonRotAngles(1:end-1);
 
-SensorPos_max = 1.1*sqrt(2)*axis_length;
+SensorPos_max = ceil(1.1*sqrt(2)*axis_length/RadonDeltaS)*RadonDeltaS;
 RadonSensorPos = (-SensorPos_max:RadonDeltaS:SensorPos_max);
 
 RadonData = zeros(length(RadonSensorPos),length(RadonRotAngles));