# Função para  a seleção  de  pares de elementos
# de valores comerciais que podem ser associados
# para obter-se  um valor próximo de um desejado
#
# por:	Nicolau Werneck<nwerneck@cefala.org
#	http://www.cefala.org/~nwerneck

function sai = comercial (quem)
	E12=[10 12 15 18 22 27 33 39 47 56 62 68 82];
	ex=floor(log10(quem));
	mant=quem/10**ex;
	escopo = [realmin 12/100 E12/10 E12 realmax].*  (10**ex); 
	arr = zeros((length(escopo))**2,6);
	int = 0;
	for A = (1:length(escopo))
		for B = (A:length(escopo))
			a = escopo(A); b=escopo(B);
			arr(++int,:)=[abs(log10((a+b)/quem))  ((a+b)/quem)    (a+b) a b 0];
			arr(++int,:)=[abs(log10((1/(1/a+1/b))/quem)) ((1/(1/a+1/b))/quem)   (1/(1/a+1/b))  a b 1];
		endfor;
	endfor;
	printf("\n");
	printf("==[Valor desejado: %9.3e ]======================\n", quem);
	[S, IN] =sort(arr(:,1));
	if S(1) > log10(1.05)
		printf("muito toscão! :(((\n");
	endif;

	printf("Desvio\tValor    \tErro  \t R_a     op  R_b    \n");
	int=0;
	odeiomatlab = [" +"; "//"];
	while (S(++int) <= log10(1.05))
		if arr(IN(int),6) 
			a = arr(IN(int),4);
			b = arr(IN(int),5);
			R = arr(IN(int),3);
ptot = (    (b/(a+b))**2*b*0.1 + (a/(a+b))**2*a*0.1     ) * ((1/a)+(1/b));
		else
			ptot = 0.8;
		endif;
		printf("%+.2f%%\t%.3e\t%05.2f%%\t%8.1e %s  %.1e\n",
			(1-arr(IN(int),2))*100,
			arr(IN(int),3),
			abs((1-ptot)*100),
			arr(IN(int),4),
			odeiomatlab(1+arr(IN(int),6),:),
			arr(IN(int),5));
	endwhile;
	
endfunction;