用矩阵的思路求解啊

把二次部分和一次部分转换为二次表达式，比如x^2+6x 转换为（x+3)^2然后把这个部分看出一个新的变量z，再用矩阵求系数的方法求解。用solve函数：
> rf<-c(1,2,3)
> lf<-matrix(c(1,23,43,54,1,23,4,5,6),3)
> solve(lf,rf)
[1] 0.037543801 0.001001168 0.227098281