tablekkTxe
d
disp次牛顿插值多项式NxvpaNx4
4龙贝格积分myrombergmmyrombergmclcclearall积分区间a0b1aa1010a取极小值以使得f01，避免si
0o等faT_1bafafb2T_2T_12ba2faba2T_4T_22ba4faba4faba34T_8T_42ba8faba8faba38faba58aba78S_1T_2141T_2T_1S_2T_4141T_4T_2S_4T_8141T_8T_4C_1S_21421S_2S_1C_2S_41421S_4S_2R_1C_21431C_2C_1Tsum0R_kR_102k3开始迭代whileabsR_kR_11e6精度设为1e6
R_1R_kfori12k
TsumTsumfa2i1ba2k1e
dT_kT_82ba2k1TsumS_kT_k141T_kT_8C_kS_k1421S_kS_4R_kC_k1431C_kC_2
fT_8T_kS_4S_kC_2C_kTsum0kk1e
dkformatlo
gR_k
5四阶龙格库塔法myRKmmyRKmclcclearallh005
5hyzeros3
y1132K1zeros31K2zeros31K3zeros31K4zeros31fori1
xihforj13
K11hy2iK12hy3iK13hy3iy2iy1i2x3K21hy2iK122K22hy3iK132K23hy3iK132y2iK122y1iK1122xh23K31hy2iK222K32hy3iK232K33hy3iK232y2iK222y1iK2122xh23K41hy2iK32K42hy3iK33K43hy3iK33y2iK32y1iK312xh3e
dfork13yki1ykiK1k2K2k2K3kK4k6e
de
dxhhh
1hplotxhy1rholdo
xih001
1hyixiexpxi2xi1plotxiyikholdo
lege
d4级RKyxxex2x1
fr