c	***********************************************************************
	subroutine sypvvv(a,b,c,d,n)
c	a(i)=b(i)+c(i)*d(i)
c       jul 2011 malv+fgg
c	***********************************************************************
	real*8 a(n),b(n),c(n),d(n)
	integer n,i
	do 1,i=2,n-1
	  a(i)= b(i) + c(i) * d(i)
 1	continue
	a(1) = 0.0d0
	a(n) = 0.0d0
	return
	end