varbvszbin.m appears to be working correctly (I'm testing the program
in the script examplebin2.m.) So we are all set for the variational
inference method with covariates.