数值计算编程20个.docx

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 A3-x - l)*(aA3-a- l)0Thcn b = x Else a = x End If Loop Until Abs(x A3-x- l)s Print x Print x人3?x?1 End Sub 结果：1. 2、 P12,第五题 Private Sub Fomi_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.986 3、 程序： Private Sub Form_Click() Dim i As Integer Dim L a(0 To 3). p As Double 1 = 2 For i = 0 lb 3 a(i) = InputBox(Hx(w i ”)=”) Next! P = a(3) For i = 2 To 0 Step -1 p = p ? t + a(i) Nexti 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-a0.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=l/2/10 i= I Do s = s + p p = p/2*(-l)/100*(2*i.l) i = i+l Loop Until p0.01 Print s End Sub 结果：10.05 6、 拉格朗日： Private Sub Form_Click() Din】i, j. n As Integer Dim x(l To 100). y(l To 100), p, s, I As Double t = 2.5 n = 3 For i = 1 lb n x(i) = InputBox(”x(” i ”)=”) y(i) = InputBox(”y(” i ”)=”) Nexti s = 0 For i = 1 lb n p= For j = 1 lb n IfiojThcn p = p*(t-x(j))/(x(i)-x。)) End If Nextj s = s + p*y(i) Next i Print s End Sub 结果：6 7、 课堂埃尔米特插值程序： Private Sub Form_Activate() Dim xO, yO, xl, yl. gO, gl As Double x0 = 0 y0 = 0 xl = 1 yl = l g0 = 0 gl =2 x = 0.5 a0 = ?2/(x0?xl)人3 bO = 1 / (xO - xl)A 2 + 2 ? xO / (xO - xl)A 3 al =-2/(x1 -xO) A3 bl =2*xl /(xl - xO) A3+ 1 / (xl - x0)A2 k0= l/(xO-xl)A2 kl = 1/(x1 -xO)A2 10=(x-xl)A2*(a0*x + b0) 11 =(x-xO)A2*(al *x + bl) ggO = kO*(x-xO)*(x-xl)A2 ggl =kl ?(x-xO)A2*(x-xl) y = 10*)0 +Il *yl + ggO * gO + ggl * gl Print y End Sub 结果：0.25 8、 第一个：样条函数： Private Sub Form_Activate() Dim x(3), y(3), mO. ml. m2 As Double Dim i As Integer For i = 0 lb 2 x(i) = i R)=i*i Nexti m0 = 0 m2 = 4 hO = x(l)-x(O) hl =x(2)?x(l) a = 2 / hO b = 4/h0 + 4/hl c = 2/hl d = (-6*y(l) + 6* y(2)