y_eq_yx.F
c**********************************************************************
#include "author.inc"
c* $Id: y_eq_yx.F,v 1.5 1995/07/30 17:49:43 turner Exp $
c*
c* Computes y = y*x
c*
c* <PARAMETER LIST>
c*
c* Input:
c* n - number of elements
c* x - x-vector
c*
c* In/Out:
c* y - y-vector
c*
c* Output:
c* status - return status
c*
#include "copyright.inc"
c**********************************************************************
subroutine JT_y_eq_yx (n, x, y, status)
implicit none
c
c ... Input:
integer n
real x(n)
c
c ... In/Out:
real y(n)
c
c ... Output:
integer status
c
c ... Local:
integer i
c
do i=1,n
y(i) = y(i)*x(i)
enddo
c
status = 0
return
end