BCESolve
fparandomrp.cpp