void sort(int n, double ra[])
{
    int i,j,l,ir;
	double rra;
	l = int(n / 2) + 1;
    ir = n;
    do
	{
        if (l > 1)
		{
            l = l - 1;
            rra = ra[l];
		}
        else
		{
            rra = ra[ir];
            ra[ir] = ra[1];
            ir = ir - 1;
            if (ir == 1)
			{
              ra[1] = rra;
              return;
            }
        }
        i = l;
        j = l + l;
        while (j <= ir)
		{
            if (j < ir)
			{
                if (ra[j] < ra[j + 1])
				{
					j+=1;
				}
            }
            if (rra < ra[j])
			{
                ra[i] = ra[j];
                i = j;
                j += j;
			}
            else
			{
                j = ir + 1;
            }
        }
        ra[i] = rra;
    }while(1);
}