##USAGE	recf_zigf_X( { N=3 } )
##	select N=1 or N=2 or N=3
##	recf_zigf - demo
##AUTHOR (C) 2000 Rolf Fabian <fabian@tu-cottbus.de> 0519

function recf_zigf_X(N)

	if	nargin<1 N=3; end

	M=501; t=linspace(-4,4,M)'; 

	gset key;
	gset grid;
	xlabel('x');
	ylabel('y');

	if	(N==1)
		title('recf example')
		gset yrange [-2:2];
		plot(	t,recf(t)  ,';cos phase (M=0);',t,cos(t),';;',\
			t,recf(t,1),';sin phase (M=1);',t,sin(t),';;')
	elseif	(N==2)
		gset yrange [-2:2];
		title('zigf example')
		plot(	t,zigf(t)  ,';cos phase (M=0);',t,cos(t),';;',\
			t,zigf(t,1),';sin phase (M=1);',t,sin(t),';;')
	elseif	(N==3)
		gset yrange [-4:4];
		plot(	t,tan(t),';tan(t);',\
			t,zigf(t,1)./zigf(t),';zigf(t,1) / zigf(t);',\
			t,zigf(t,1)./recf(t),';zigf(t,1) / recf(t);',\
			t,recf(t,1)./zigf(t),';recf(t,1) / zigf(t);',\
			t,recf(t,1)./recf(t),';recf(t,1) / recf(t);');
	else
		disp('unknown example tag .. select N=1 -or- N=2 -or- N=3')
	end
