数值计算编程20个.docVIP

1、P11，第一题 Private Sub Form_Click() Dim a, b, x, s As Double a = 1 b = 2 s = 0.001 Do x = (a + b) / 2 If (x ^ 3 - x - 1) * (a ^ 3 - a - 1) 0 Then b = x Else a = x End If Loop Until Abs(x ^ 3 - x - 1) s Print x Print x ^ 3 - x - 1 End Sub 结果：1.32470703125 2、P12，第五题 Private Sub Form_Click() Dim a, x, s As Double a = 8 s = 0.001 x = a Do p = x x = (p + a / p) / 2 Loop Until Abs(x - p) s Print x End Sub 结果：2.82842712504986 3、程序： Private Sub Form_Click() Dim i As Integer Dim t, a(0 To 3), p As Double t = 2 For i = 0 To 3 a(i) = InputBox(x( i )=) Next i p = a(3) For i = 2 To 0 Step -1 p = p * t + a(i) Next i Print p End Sub 结果：49 4、程序： Private Sub Form_Click() Dim a, b, s As Double a = 0.2 b = 1.34 s = (a + b) / 2 Do If a * (1 - a) * s * (1 - s) 0 Then b = s Else a = s End If s = (a + b) / 2 Loop Until b - a 0.01 Print s End Sub 结果：0.997109375 5、程序： Private Sub Form_click() Dim x, s, p As Double Dim i As Integer x = 101 s = 10 p = 1 / 2 / 10 i = 1 Do s = s + p p = p / 2 * (-1) / 100 * (2 * i - 1) i = i + 1 Loop Until p 0.01 Print s End Sub 结果：10.05 6、拉格朗日： Private Sub Form_Click() Dim i, j, n As Integer Dim x(1 To 100), y(1 To 100), p, s, t As Double t = 2.5 n = 3 For i = 1 To n x(i) = InputBox(x( i )=) y(i) = InputBox(y( i )=) Next i s = 0 For i = 1 To n p = 1 For j = 1 To n If i j Then p = p * (t - x(j)) / (x(i) - x(j)) End If Next j s = s + p * y(i) Next i Print s End Sub 结果：6 7、课堂埃尔米特插值程序： Private Sub Form_Activate() Dim x0, y0, x1, y1, g0, g1 As Double x0 = 0 y0 = 0 x1 = 1 y1 = 1 g0 = 0 g1 = 2 x = 0.5 a0 = -2 / (x0 - x1) ^ 3 b0 = 1 / (x0 - x1) ^ 2 + 2 * x0 / (x0 - x1) ^ 3 a1 = -2 / (x1 - x0) ^ 3 b1 = 2 * x1 / (x1 - x0) ^ 3 + 1 / (x1 - x0) ^ 2 k0 = 1 / (x0 - x1) ^ 2 k1 = 1 / (x1 - x0) ^ 2 l0 = (x - x1) ^ 2 * (a0 * x + b0) l1 = (x - x0) ^ 2 * (a1 * x + b1) gg0 = k0 * (x - x0) * (x - x1) ^ 2 gg1 = k1 * (x - x0) ^ 2 * (x - x1) y = l0 * y0 + l1 * y1 + gg0 * g0 + gg1 * g1 Print y End Sub 结果：0.25 8、第一个：样条函数： Private Sub Form_Activate() Dim x(3), y(3), m0, m1, m2