use subroutine in Netlib