改主程序所调用的函数中的表达式来实现对其它函数的任意初值条件求微分计算。
3程序中运用了大量的for循环语句，因为该公式中涉及大量的求和，且有不同的函数和对不同的数值求值，编程稍显繁琐。所以编写过程中一定要注意各循环的次数，以免出错。
14
f5
123841221152371061074111233601135840718719174238230678132031743
211523719141823312543210123452189736156384907841651112348
3123124

1061074312543215567914312384820314541836742105678103369931010103
111233610123453123848271084374101011374185621010230718280037585
A0113584218973620314544101011198979180431637311122321213141784317


071871915638491836742374185604316379789365010345811034560238417
17423820784165105678121010233111223010345814713846531237892213474
306781311123480336993071828212131411034563123789307193344446782
203174331231241010103003758517843170238417221347444467824000001
b21874369339923182517341708467169517843178661234311101230471934556784392T
用列主元消去法求解Axb。
51理论依据：
列主元素消元法是在应用Gauss消元法的基础上，凭借长期经验积累提出的，是线性方程组一般解法，目的是为避免在消元计算中使
误差的扩大，甚至严重损失了有效数字使数据失真，而在每次初等变
换前对矩阵作恰当的调整，以提高Gauss消元法的数字稳定性，进而提高计算所得数据的精确度。即在每主列中取绝对值最大的元素作主
元，再做对应的行交换然后消元求解的办法。具体做法如下：
将方阵A和向量b写成C（A，b）。将C的第1列中第1行的
元素与其下面的此列的元素逐一进行比较，找到最大的元素cj1，将第
j行的元素与第1行的元素进行交换，然后通过行变换，将第1列中
第2到第
个元素都消成0。将变换后的矩阵C1的第二列中第二行
的元素与其下面的此列的元素逐一进行比较，找到最大的元素
c1k2
，
将第k行的元素与第2行的元素进行交换，然后通过行变换，将第2
列中第3到第
个元素都消成0。以此方法将矩阵的左下部分全都消
15
f成0后再求解。最终形式如下：


a
110
Ab）～
0
52C语言程序代码

a
1

g1


0
a


g


（1）比较该列的元素的绝对值的大小，将绝对值最大的元素通过行
变换使其位于主对角线上；
（2）进行高斯消去法变换，把系数矩阵化成上三角形，然后回代求
i
clude