void simp1(double a[], int mp, int np, int mm, int ll[], int nll,
		   int iabf, int& kp, double& bmax)
{
	double test;
    kp = ll[1];
    bmax = a[mm*np+kp + 1];
    for (int k = 2; k<=nll; k++)
	{
        if (iabf == 0)
		{
            test = a[mm*np+ll[k] + 1] - bmax;
		}
        else
		{
            test = fabs(a[mm*np+ll[k] + 1]) - fabs(bmax);
        }
        if (test > 0.0)
		{
            bmax = a[mm*np+ll[k] + 1];
            kp = ll[k];
        }
    }
}