MACRO get_contents id=160011 * * get contents into vector contents * set fit 1; set stat 1111111 * CLOSE 0; HI/DEL 0; VE/DEL * * HI/FILE 1 [HBOOK].hbook npts1 = 121 npts = [npts1] - 1 xmin = 1.4 xmax = 2.0 sigma jj=array([npts],0#[npts]-1) vector/print jj vector/create contents([npts]) R [npts]*0. vector/create errors([npts]) R [npts]*0. vector/create x([npts]) R [npts]*0. sigma x = [xmin] + ([xmax]-[xmin])*jj/[npts] +0.0025 HI/GET_VECT/CONTENTS [id] contents hi/get_vect/errors [id] errors * vector/print contents * vector/print errors vector/print x vector/write x,contents,errors moskov_[id].dat (3f10.4) do j=1,[npts] y = $eval(contents([j])) yerr = $eval(errors([j])) * xx = [xmin] + ([j]-0.5)*[bin] * mess [xx] [y] [yerr] enddo RETURN