O pacote `lpSolveAPI` aceita restrições de três tipos, `"<="`, `">="` e `"="`. Creio que na pergunta a demanda deve ser entendida com *demanda máxima*, não como demanda fixa. Então as restrições a serem usadas serão de menoridade, `"<="`

    #x1 = 50, pelo x2 = 50, pelo x3 = 80 e pelo x4 = 120
    add.constraint(lprec, 1, "<=", 50, indices = 1) #restricao 4a
    add.constraint(lprec, 1, "<=", 50, indices = 2) #restricao 4b
    add.constraint(lprec, 1, "<=", 80, indices = 3) #restricao 4c
    add.constraint(lprec, 1, "<=", 120, indices = 4) #restricao 4d

Agora resolve-se o problema tal como está na pergunta. As soluções serão dadas por `get.variables(lprec)`.

    solve(lprec) #resolver ppl
    #[1] 0
     
    get.objective(lprec) #funcao objetivo resolvida
    #[1] 1948.62
    
    get.variables(lprec) #qto de cada variavel otimiza o sistema
    #[1] 2.6654343 0.5221156 0.0000000 0.2073184

Verificar estas soluções multiplicando os coeficientes da função objetivo pelo vetor solução.

    c(600, 550, 400, 300) %*% get.variables(lprec)
    #        [,1]
    #[1,] 1948.62