A0:=(2-120) 初始值 解方程组:Ac=y0 0.96146-3.8785i C: =Isolve(A, yo) C=0.96146+3.8785i 0.07708 下面生成解曲线的图形 E:=0.01t0:=0te:=6 t:= t0t0+8. te 2 1.401 1.042 下面使用 Mathcad的内部函数 rkfixed来求解: 定义初始条件 定义与未知解函数有关的向量Dty): 阶导数 二阶导数 y2ty1-yo 三阶导数 0 021 1:=0.n 需计算的点数 1188 31.64 F: =rkfixed(y, to, te, n, D) 解矩阵 51401 61281 F的第一列为自变量的值,第二列为对应于自变量值的解函数值 71.162 第三列和第四列分别为解函数的一阶导数和二阶导数值 90923 100805 0.686F的第一列为自变量的值, 第二列为对应于自变量值的解函数值, 第三列和第四列分别为解函数的一阶导数和二阶导数值. F i, 1 0 0 1 2 3 4 5 6 7 8 9 10 11 2 1.88 1.76 1.64 1.52 1.401 1.281 1.162 1.042 0.923 0.805 0.686 = F := rkfixed(y, t0, te, n,D) 解矩阵 n i := 0.. n 需计算的点数 te - t0 e := 三阶导数 D(t, y) 二阶导数 y 1 y 2 -2 1 4 - y 2 × y 1 + y 0 - æ ç è ö ÷ ø × é ê ê ê ê ë ù ú ú ú ú û := 一阶导数 定义与未知解函数有关的向量D(t,y): y 定义初始条件 2 -12 0 æ ç ç è ö ÷ ÷ ø y := y0 = 下面使用Mathcad的内部函数rkfixed 来求解: y(t,C) 2 1.88 1.76 1.64 1.52 1.401 1.281 1.162 1.042 0.923 = 0 1 2 3 4 5 6 10 y(t,C) t e := 0.01 t0 := 0 te := 6 t := t0, t0 + e.. te 下面生成解曲线的图形 C 0.96146 - 3.8785i 0.96146 + 3.8785i 0.07708 æ ç ç è ö ÷ ÷ ø C := lsolve(A, y0) = 解方程组: A*C=y0 y0 ( 2 -12 0 ) 初始值 T :=