QR decomposition in R is incorrect?

Consider the following code, where lm is used for the purpose of quickly form a matrix only :

Error in ROI_solve(op) : no solver found for this signature when using R ROI package with matrix output

I’m trying to calculate the max expected value for the following scenario: