假设一组常见的性常微分方程组,如下所述.docVIP

假设一组常见的线性常微分方程组，如下所述 (5.1) 有如下边值条件 (5.2) 这里和是上连续解析函数，是特征值，是实常数。等式(5.1)和(5.2)构成了常微分方程特征值问题。表示。 SUBROUTINE FMN(……) SUBROUTINE IMMEI (……) 输入参数: A ,B: 区间的两端点值，见(5.1)式 N: 区间的等分数 IR: 方程的个数，即(5.1)式中的 IT： 边界条件的数目，即(5.2)式中的 NUM: 特征值按实部从小到达排列，要求取前面NUM个实部较小的特征值 Mpart:要求打印在节点上的值 Model:积分矩阵的选择 1．model=9 拉格朗日多项式插值，此时n=3或4 2．model=1 线性插值 3．model=2 抛物线插值 4．model=23 改进的抛物线插值 5．model=32 三次样条插值 IRE(J): 记录变量在方程中不含特征值项的最高阶次，即(5.1)式中的 ILQ(J): 记录变量在方程中含特征值项的最高阶次，即(5.1)式中的 IKG(I,J,K): 对于第i中方程变量第k阶倒不含特征值项系数为零时，即(5.1)式中时，IKG(I,J,K)=0；否则等于1 IKQ(I,J,K): 对于第i中方程变量第k阶倒含特征值项系数为零时，即(5.1)式中时，IKQ(I,J,K)=0；否则等于1 GMN(I,J,K,L): 记录第i个方程变量第k阶倒在第L节点上不含特征值项的系数, 对(5.1)式,GMN(I,J,K,N)＝， QLN(I,J,K,L): 记录第i个方程变量第k阶倒在第L节点上含特征值项的系数, 对(5.1)式,QLN(I,J,K,N)＝， EM(L,J):记录第L个边界条件变量的阶次信息，EM(L,J)中的元素为0.XXXXXXXX，小数点后第1~4位记录该边界条件中不含特征值的每项的阶次k，项数不足4位时，以数字8充之(截止)；小数点后第5~8位记录该边界条件中含特征值的每项的阶次k，项数不足8位时，以数字8充之(截止)。 MX(L,J): 记录第L个边界条件变量的节点号c，如果，MX(L,J)＝0，如果，MX(L,J)＝-1或n。 EMX(L,J,K): 记录第L个边界条件中不含特征值的每项前的系数 BMX(L,J,K): 记录第L个边界条件中含特征值的每项前的系数 EM(L,J) ,MX(L,J) ,EMX(L,J,K) 和BMX(L,J,K)中的L和J的顺序相同。EMX(L,J,K) 和BMX(L,J,K)中的K顺序相同。 输出参数: EVR(I):存放按特征值实部从小到达排列的第I个特征值的实部 EVI(I): 存放按特征值实部从小到达排列的第I个特征值的虚部 YKNB(J,K,L):每个特征值中变量在节点号L处的特征向量的实部 YKMN(J,K,L): 个特征值中变量在节点号L处的特征向量的虚部