在R中嵌套for循环 - 发布最终结果

I am in the midst of solving a problem in Reconstructing (or recovering) a probability distribution function when only the moments of the distribution are known. I have written codes in R for it and although the logic seems right to me, I am not getting the output that I want. I am in the midst of solving a problem in Recon