Functions vsva.R and vsva_homo.R implement the VAR for the 
heteroscedastic and homoscedastic model respectively. Function 
vbhetlm1.R does the lower bound maximization and is a helper function 
for the others, examples.R gives some examples.
