subroutine zblas() d = dcabs1(z) call drotg(da,db,c,s) d = dzasum(n,zx,incx) d = dznrm2( n, zx, incx) call zdrot (n,zx,incx,zy,incy,c,s) i = izamax(n,zx,incx) call zaxpy(n,za,zx,incx,zy,incy) call zcopy(n,zx,incx,zy,incy) z = zdotc(n,zx,incx,zy,incy) z = zdotu(n,zx,incx,zy,incy) call zdscal(n,da,zx,incx) d = zmach(job) call zrotg(ca,cb,c,s) call zscal(n,za,zx,incx) call zswap (n,zx,incx,zy,incy) stop end