function Vo=voxel(tau1,tau2,zeta)
%
% IN  : a pixel (tau1,tau2,zeta)
% OUT : the voxels in the sample that are connected to the pixel P
%       V(i,n) is the i-th component of the n-th voxel
% All the coordinates are given in the frame (C,nr,np,nc)
%
% (c) Stephane Balac - MAPLY - INSA Lyon - mars 2001


global r1
global r2
global g
global adim

Vo = [];
R=[];
% computation of the roots of F
R = fndrc(tau1/adim,tau2/adim,zeta/adim,r1/adim,r2/adim,adim);
R = adim*R;
% computation of the coordinates of the voxels Vo linked to P
for i=1:length(R)
  xs = R(i);	
  xr = tau1 + (xs-zeta)*g(3)/g(1);
  xp = tau2; 
  Vo(1,i) = xr;
  Vo(2,i) = xp;
  Vo(3,i) = xs;
end



