/* ---------------------------------- */
/* save as x_init.h                  */
/* ---------------------------------- */
double *I_p(
int col)
{
double *P;
int i;

    ++col;

    P = (double*) malloc( col * sizeof(*P));

    if(!P) {
            printf(" I was unable to allocate "
                   "the memory you requested.\n\n"
                   " double *I_p();\n\n"
                   " P = malloc(col * sizeof(*P));\n\n");
            fflush(stdout);
            getchar();
            exit(EXIT_FAILURE);
            }

    P[0] = col;
    for(i=1; i<col; i++) P[i] = 0;

 return(P);
}
/* ----------------------------------------------------- */
void c_a_P(
double  *a,
double  *Px
)
{
int i=1;

 	for (; i<Px[0]; i++)

 	         Px[i]=a[i-1];
}