1 subroutine xzdotu (n, zx, incx, zy, incy, retval) 2 double complex zdotu, zx(*), zy(*), retval 3 integer n, incx, incy 4 external zdotu 5 retval = zdotu (n, zx, incx, zy, incy) 6 return 7 end 8