double *dvector(long nl, long nh);
void free_dvector(double *v, long nl, long nh);

z = dvector(0,2*n+10);