Usage:

1)  CALL INIT(NO1,NV1)  : NO1 is the order of the Polynomial and NV1 the number of variables

2) CALL INIT(NO1=2,ND1=1,NP1=0,NDPT1 =0) 

Here ND1 is the number of degrees of freedom (1,2, or 3),  NP1  is the number of extra variables besides the phase space variables.

Total Number of Variables NV=2*ND1+NP1

NDPT1 is the position of the energy variable if and only if the variable is constant and part of the phase space.  It can take only two values in that case: NDPT1=2*ND1 or 2*ND1-1.

Otherwise NDPT1 is set to zero.