vector linear( double xmin, xmax, int n)
/*
 * create a vector of elements in a linear scale of n points from xmin to xmax
 */
{
    return xmin + [0:n-1]*((xmax - xmin)/(n - 1));
}
