subroutine dblas() d = dasum(n,dx,incx) call daxpy(n,da,dx,incx,dy,incy) call dcopy(n,dx,incx,dy,incy) d = ddot(n,dx,incx,dy,incy) d = dmach(job) d = dnrm2 ( n, dx, incx) call drot (n,dx,incx,dy,incy,c,s) call drotg(da,db,c,s) call dscal(n,da,dx,incx) call dswap (n,dx,incx,dy,incy) i = idamax(n,dx,incx) stop end