const double beta, <- 
\*= beta; <- = 0;
= beta .*; <- = 0;
