SUBROUTINE sort3(n,ra,rb,rc,wksp,iwksp) INTEGER n,iwksp(n) REAL ra(n),rb(n),rc(n),wksp(n) CU USES indexx INTEGER j call indexx(n,ra,iwksp) do 11 j=1,n wksp(j)=ra(j) 11 continue do 12 j=1,n ra(j)=wksp(iwksp(j)) 12 continue do 13 j=1,n wksp(j)=rb(j) 13 continue do 14 j=1,n rb(j)=wksp(iwksp(j)) 14 continue do 15 j=1,n wksp(j)=rc(j) 15 continue do 16 j=1,n rc(j)=wksp(iwksp(j)) 16 continue return END C (C) Copr. 1986-92 Numerical Recipes Software .