static int irbit2(int& iseed)
{
	int ib1,ib3,ib5,ib14,temp,mask=21;
    ib1 = 1;
    ib3 = 4; 
	ib5 = 16;
	ib14 = 8192;
    if ((iseed & ib14) != 0)
	{
        iseed = iseed ^ mask;
        if (iseed > pow(2,14 ))
			 iseed = iseed - int(pow(2, 14));
        iseed = 2 * iseed | ib1;
        temp = 1;
	}
    else
	{
        if (iseed >pow(2,14 ) )
			iseed = iseed -int( pow(2,14 ));
        iseed = 2 * iseed & (~ ib1);
        temp = 0;
    }
	return temp;
}