** Copyright (c) 1989, NVIDIA CORPORATION.  All rights reserved.
**
** Licensed under the Apache License, Version 2.0 (the "License");
** you may not use this file except in compliance with the License.
** You may obtain a copy of the License at
**
**     http://www.apache.org/licenses/LICENSE-2.0
**
** Unless required by applicable law or agreed to in writing, software
** distributed under the License is distributed on an "AS IS" BASIS,
** WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
** See the License for the specific language governing permissions and
** limitations under the License.

*--- Induction initial values, update reaching defs
*    

	program p
	parameter (N=16)
	integer i, j
	double precision expect(N)
	double precision result(100)
	common /c/ result

	call cnstnt(4,result)

	data expect /
     +    1.0, 1.0, 1.0, 1.0, 1.0, 2.0, 3.0, 0.0,
     +    1.0, 3.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0 /

	call checkd(result, expect, N)
	end
C***********************************************************************
      SUBROUTINE CNSTNT(N,C)
      IMPLICIT DOUBLEPRECISION(A-H,O-Z)
      DIMENSION C(N,N)
C
C.....MOLECULAR CONSTANTS FOR WATER IN ANGSTROM, RADIAN, AND A.M.U.
C
      C(1,1)=1.0D0
      DO 1000 N1=2,N
	  NN=N1-1
	  TN=NN
	  C(1,N1)=1.0D0
	  TK=1.0D0
	  DO 1100 K1=2,N1
	      C(K1,NN)=C(K1-1,NN+1)*TN/TK
	      NN=NN-1
	      TN=TN-1.0D0
	      TK=TK+1.0D0
1100      CONTINUE
1000  CONTINUE
	return
	end