嵌套不就解决了吗

@Abomt- @郑喵喵😺 @什么什么情况😂 @iew😘 @小可心丫 这几位是我哈佛的师姐，女博士

赋值HIJK，另外注意给定p初值p0。这里的求和项我估计是指定p(...k)=C，这样的形式，没有实际数据不清楚，如果是这样的形式，用4个行向量给出来:ph,pi,pj,pk。
for h=1:H
for i=1:I
for j=1:J
for k=1:K
p1=p0+pk(k)
end
p2=p1=pj(j)
end
以此类推，说实话，数学里的代码这种东西，一次编不太可能对得上，我估计十有八九是错的，慢慢调试才行