image=zeros(700,700,200);
p=0;
tic
for theta=0:2.81:357.19
   p=p+1;
   rad=theta*pi/180;
    prothe=proj(:,:,p);
    for z=1:1:200
    for i=1:1:500
    for j=1:1:500
    x=round((j+100)+((200-z)/200)*100*cos(rad));
        y=round((i+100)-((200-z)/200)*100*sin(rad));
        image(y,x,z)=image(y,x,z)+(1/200)*prothe(i,j);
    end
    end
    end
end
toc
di=max(max(max(image)));
image=image/di;
