数值计算与程序设计.docVIP

1、用 Fortran或Basica在计算机上计算注意避免溢出算法。 Basic程序编写程序如下: !首先建立一个command,及一个text及两个说明的label Private Sub Form_Load() Command1.Caption = 计算 Label1.Caption = y=29^71/71! Label2.Caption = 避免溢出计算y值为: Text1.Text = End Sub Private Sub Command1_Click() Dim y y =1 For i = 1 To 71 y = y * 29 / i Next i Text1.Text = y End Sub 2、画出牛顿插值多项式程序结构流程图并用一种高级语言编制相应的程序计算函数f(x)的近似牛顿插值多项式N5(x), x=0.2+0.01k, k=0~20． 0．2 0．24 0．28 0．32 0．36 0．40 0.1987 0.2377 0.2764 0.3146 0.3523 0.3894 程序及框图如下: Private Sub Form_Load() N = 0 End Sub Dim x, y, y1 ,u(9), v(9), t(9) As Double Dim c ,e,N Private Sub Command1_Click() If N = 9 Then MsgBox (你已经输入了9组数据，牛顿插值的精度在九阶之内!) Exit Sub End If If Text1 = Then Text1 = 0 If Text2 = Then Text2 = N u(N) = Text1 v(N) = Text2 Text1 = Text2 = Label2.Caption = X( (N) )= (u(N)) Label3.Caption = Y( (N) )= (v(N)) N = N + 1 Text1.SetFocus End Sub Private Sub Command2_Click() c = 1 y1 = v(0) e = Text3 x = Text4 For k = 1 To N - 1 For i = k To N - 1 t(i) = (v(i) - v(i - 1)) / (u(i) - u(i - k)) Next c = c * (x - u(k - 1)) y1 = y1 + c * t(k) If Abs(y - y1) = e Then For i = k To N - 1 v(i) = t(i) Next y = y1 Else: Exit For End If Next Text5 = y1 End Sub 3 用选列主元的高斯消去法求解下列方程组： Private Sub Command1_Click() Dim a(3, 4), x(3), r, s, e As Double Dim i, j, k, p, flag As Single flag = 1 For i = 0 To 4 ‘输入增广矩阵A a(0, i) = Text1(i) a(1, i) = Text2(i) a(2, i) = Text3(i) a(3, i) = Text4(i) Next i 选主元过程 For k = 0 To 2 1 p = k: e = Abs(a(k, k)) For i = k + 1 To 3 2 If Abs(a(i, k)) e Then e = Abs(a(i, k)): p = i End If Next i For j = k To 4 3 ‘交换k行和p行 s = a(k, j): a(k, j) = a(p, j): a(p, j) = s Next j 3 If a(k, k) = 0 Then ‘消元 Print (GAUSS METHOD DOES NOT RUN!) flag = 0 Else For i = k + 1 To 3 4 r = a(i, k) / a(k, k) If a(i, k) 0 Then For j = 0 To 4 5 a(i, j) = a(i, j) - r * a(k, j) Next j End If Next i End If Next