function fac3d1(x,y,z,T,A,leg,flags,ebox)
	
	// Copyright INRIA
	 
	[lhs,rhs]=argn(0)
	if rhs <= 0 then   // demo
		
		title_demo = [
			'';
			'Demo of fac3d1()';
			'========================================';
			'';
			'WARNING : fac3d1 is obsolete see plot3d plot3d1 instead';
			''];
		
		s_mat=[
			'deff(''[x,y,z]=sph(alp,tet)'',[''x=r*cos(alp).*cos(tet)+orig(1)*ones(tet)'';';
			'''y=r*cos(alp).*sin(tet)+orig(2)*ones(tet)'';';
			'''z=r*sin(alp)+orig(3)*ones(tet)'']);';
			'r=1;orig=[0 0 0];';
			'[x1,y1,z1]=eval3dp(sph,linspace(-%pi/2,%pi/2,40),linspace(0,%pi*2,20));';
			'[n1,m1]=size(x1);';
			'r=1/2;orig=[-1 0 0];';
			'[x2,y2,z2]=eval3dp(sph,linspace(-%pi/2,%pi/2,40),linspace(0,%pi*2,20));';
			'[n2,m2]=size(x2);';
			'x=[x1 x2];y=[y1 y2];z=[z1 z2];';
			'plot3d1(x,y,z);';];
		
		write(%io(2),title_demo);
		write(%io(2),s_mat);
		write(%io(2),' ');
		execstr(s_mat);
		return
	end
	
	write(%io(2),'WARNING : fac3d1 is obsolete see plot3d plot3d1 instead');
	
	select rhs 
		case 3 then plot3d1(x,y,z) 
		case 4 then plot3d1(x,y,z,T) 
		case 5 then plot3d1(x,y,z,T,A) 
		case 6 then plot3d1(x,y,z,T,A,leg) 
		case 7 then plot3d1(x,y,z,T,A,leg,flags) 
		case 8 then plot3d1(x,y,z,T,A,leg,flags,ebox) 
	end
endfunction
