n=100;
density=.1;
lb=-10;
ub=10;
A=sprandsym(n,density);
B=sprandsym(n,density);
[xv,lmb,iresult]=sptarn(A,B,lb,ub);
[xve,lmbe]=eigs(A,B,3);
