麻烦各位高手帮忙根据下题编写MATLAB程序，谢谢！
 

以39个受辐射的老鼠在不同时间段的死亡个数为样本数据，将样本数据分成8个时间段，求每个时间段的死亡率值（此死亡率是一组归一化的概率，由对应时间段的样本死亡个数除以总样本个数得到），在计算中采用它们的前四阶矩为约束条件。在该例中，公式（3.2）中的n=8，其对偶规划的待求变量个数为4（而原问题变量个数为8）。通过求解最大熵修匀模型的对偶规划，即可得到对应各时间段的归一化的修匀死亡率 ，将其带入（3.3）式即可求出第 个时间段修匀后的死亡老鼠个数 ，然后根据（3.4）和（3.5），可求出死亡率的修匀值 。计算结果见表3.1。
由表3.1的计算结果，给出三种方法修匀结果的曲线图，如图3.1，其中横轴代表死亡的时间，纵轴代表死亡个数的修匀值：

因表中有三个方法得到的数据，只需帮忙算出用最大熵方法得到的数据即可！谢谢！
公式和图表都在附件中，请各位高手不要嫌麻烦，帮忙一下吧，实在紧急，在此谢过！