int n = Utilities.getInt(1, max);
kl = Math.min(n, Utilities.getInt(bmax));
ku = Math.min(n, Utilities.getInt(bmax));
L = new LowerSPDBandMatrix(n, kl);
Utilities.bandPopulate(L, kl, 0);
Utilities.addDiagonal(L, 1);
while (!Utilities.spd(L))
Utilities.addDiagonal(L, 1);