InputBox(输入晶体类型)

1、Private Sub command2_Click()Dim b, line, d, v As StringDim i, j, k, i0, j0, k0, i1, j1, k1, c, n, num, l, h As IntegerDim dist, e0, e, w, boltzmanfactor As Singleb = InputBox("输入晶体类型")Open "E:文件1" For Input As #1 '找出一个中心点Do While Not EOF(1)Line Input #1, lineIf Val(Left(line,

2、 4) = 63 Theni1 = Val(Mid(line, 7, 1)j1 = Val(Mid(line, 12, 1)k1 = Val(Mid(line, 17, 1)End IfLoopClose #1Select Case bCase "sc"Open "E:文件1" For Input As #1c = InputBox("输入序号")Do While Not EOF(1)Line Input #1, lineIf Val(Left(line, 4) = Val(c) Then '将输入序号对应的坐标点作为原点Te

3、xt1.Text = linei0 = Val(Mid(Text1, 7, 1)j0 = Val(Mid(Text1, 12, 1)k0 = Val(Mid(Text1, 17, 1)End IfIf i1 - i0 >= 2 Then i0 = i0 + 2 '边界判断If i1 - i0 <= -2 Then i0 = i0 - 2If j1 - j0 >= 2 Then j0 = j0 + 2If j1 - j0 <= -2 Then j0 = j0 - 2If k1 - k0 >= 2 Then k0 = k0 + 2If k1 - k0 <

4、= -2 Then k0 = k0 - 2LoopClose #1 Open "E:文件1" For Input As #1Do While Not EOF(1)Line Input #1, linei = Val(Mid(line, 7, 4)j = Val(Mid(line, 12, 4)k = Val(Mid(line, 17, 4)dist = Sqr(i - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = 1 Then '寻找近邻num = num + 1End IfLoopClose #1Text2.Text = nu

5、mFor n = 1 To 125 '求总能量Open "E:文件1" For Input As #1Do While Not EOF(1)Line Input #1, lineIf Val(Left(line, 4) = Val(n) Theni0 = Val(Mid(line, 7, 4)j0 = Val(Mid(line, 12, 4)k0 = Val(Mid(line, 17, 4)e0 = Val(Mid(line, 23, 4)Elsei = Val(Mid(line, 7, 4)j = Val(Mid(line, 12, 4)k = Val(Mid(l

6、ine, 17, 4)e = Val(Mid(line, 23, 4)End Ifdist = Sqr(i - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = 1 Thenw = w + e * e0End IfLoopClose #1Next nText3.Text = -2 * wh = Int(125 * Rnd) '给出一个随机数Label5.Caption = hFor n = 1 To 125Open "E:文件1" For Input As #1Do While Not EOF(1)Line Input #1, li

7、neVal (Mid(line, 21, 4)If Val(Left(line, 4) = Val(n) Theni0 = Val(Mid(line, 7, 4)j0 = Val(Mid(line, 12, 4)k0 = Val(Mid(line, 17, 4)e0 = Val(Mid(line, 23, 4)If Val(Left(line, 4) = Val(Label5) Then e0 = -Val(Mid(line, 23, 4) '改变随机数对应坐标点上磁子的方向Elsei = Val(Mid(line, 7, 4)j = Val(Mid(line, 12, 4)k = V

8、al(Mid(line, 17, 4)e = Val(Mid(line, 23, 4)If Val(Left(line, 4) = Val(Label5) Then e = -Val(Mid(line, 23, 4)End Ifdist = Sqr(i - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = 1 Thenw = w + e * e0 '求改变后的能量和End IfLoopClose #1Next nText4.Text = -2 * wText5.Text = Text3.Text - Text4.Text '输出能量差值If V

9、al(Text5.Text) > 0 Then '判断磁子能否改变方向Label7.Caption = "flipped"Elseboltzmanfactor = Exp(Val(Text5.Text) / (300 * 1.38 * 10 -23)If Int(1000 * Rnd) > Val(boltzmanfactor) ThenLabel7.Caption = "unchanged"ElseLabel7.Caption = "unflipped"End IfEnd IfCase "文件2&quo

10、t;Open "E:文件2" For Input As #1c = InputBox("输入序号")Do While Not EOF(1)Line Input #1, lineIf Val(Left(line, 4) = Val(c) ThenText1.Text = linei0 = Val(Mid(Text1, 7, 4)j0 = Val(Mid(Text1, 12, 4)k0 = Val(Mid(Text1, 17, 4)End IfIf i1 - i0 >= 2 Then i0 = i0 + 2If i1 - i0 <= -2 The

11、n i0 = i0 - 2If j1 - j0 >= 2 Then j0 = j0 + 2If j1 - j0 <= -2 Then j0 = j0 - 2If k1 - k0 >= 2 Then k0 = k0 + 2If k1 - k0 <= -2 Then k0 = k0 - 2LoopClose #1 Open "E:文件2" For Input As #1Do While Not EOF(1)Line Input #1, linei = Val(Mid(line, 7, 4)j = Val(Mid(line, 12, 4)k = Val(M

12、id(line, 17, 4)dist = Sqr(i - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = Sqr(3) / 2 Thennum = num + 1End IfLoopClose #1Text2.Text = numFor n = 1 To 189Open "E: 文件2" For Input As #1Do While Not EOF(1)Line Input #1, lineIf Val(Left(line, 4) = Val(n) Theni0 = Val(Mid(line, 7, 4)j0 = Val(Mid(li

13、ne, 12, 4)k0 = Val(Mid(line, 17, 4)e0 = Val(Mid(line, 23, 4)Elsei = Val(Mid(line, 7, 4)j = Val(Mid(line, 12, 4)k = Val(Mid(line, 17, 4)e = Val(Mid(line, 23, 4)End Ifdist = Sqr(i - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = Sqr(3) / 2 Thenw = w + e * e0End IfLoopClose #1Next nText3.Text = -2 * wh = In

14、t(189 * Rnd)Label5.Caption = hFor n = 1 To 189Open "E: 文件2" For Input As #1Do While Not EOF(1)Line Input #1, lineIf Val(Left(line, 4) = Val(n) Theni0 = Val(Mid(line, 7, 4)j0 = Val(Mid(line, 12, 4)k0 = Val(Mid(line, 17, 4)e0 = Val(Mid(line, 23, 4)If Val(Left(line, 4) = Val(Label5) Then e0 =

15、 -Val(Mid(line, 23, 4)Elsei = Val(Mid(line, 7, 4)j = Val(Mid(line, 12, 4)k = Val(Mid(line, 17, 4)e = Val(Mid(line, 23, 4)If Val(Left(line, 4) = Val(Label5) Then e = -Val(Mid(line, 23, 4)End Ifdist = Sqr(i - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = Sqr(3) / 2 Thenw = w + e * e0End IfLoopClose #1Next

16、 nText4.Text = -2 * wIf Val(Text5.Text) > 0 ThenLabel7.Caption = "flipped"Elseboltzmanfactor = Exp(Val(Text5.Text) / (300 * 1.38 * 10 -23)If Int(1000 * Rnd) > Val(boltzmanfactor) ThenLabel7.Caption = "unchanged"ElseLabel7.Caption = "unflipped"End IfEnd IfCase &quo

17、t;文件3"Open "E: 文件3" For Input As #1c = InputBox("输入序号")Do While Not EOF(1)Line Input #1, lineIf Val(Left(line, 4) = Val(c) ThenText1.Text = linei0 = Val(Mid(Text1, 7, 4)j0 = Val(Mid(Text1, 12, 4)k0 = Val(Mid(Text1, 17, 4)If i1 - i0 >= 2 Then i0 = i0 + 2If i1 - i0 <= -2

18、 Then i0 = i0 - 2If j1 - j0 >= 2 Then j0 = j0 + 2If j1 - j0 <= -2 Then j0 = j0 - 2If k1 - k0 >= 2 Then k0 = k0 + 2If k1 - k0 <= -2 Then k0 = k0 - 2End IfLoopClose #1 Open "E: 文件3" For Input As #1Do While Not EOF(1)Line Input #1, linei = Val(Mid(line, 7, 4)j = Val(Mid(line, 12,

19、4)k = Val(Mid(line, 17, 4)dist = Sqr(i - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = Sqr(2) / 2 Thennum = num + 1End IfLoopClose #1Text2.Text = num For n = 1 To 365Open "E: 文件3" For Input As #1Do While Not EOF(1)Line Input #1, lineIf Val(Left(line, 4) = Val(n) Theni0 = Val(Mid(line, 7, 4)j0

20、= Val(Mid(line, 12, 4)k0 = Val(Mid(line, 17, 4)e0 = Val(Mid(line, 23, 4)Elsei = Val(Mid(line, 7, 4)j = Val(Mid(line, 12, 4)k = Val(Mid(line, 17, 4)e = Val(Mid(line, 23, 4)End Ifdist = Sqr(i - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = Sqr(2) / 2 Thenw = w + e * e0End IfLoopClose #1Next nText3.Text =

21、-2 * wh = Int(365 * Rnd)Label5.Caption = hFor n = 1 To 365Open "E: 文件3" For Input As #1Do While Not EOF(1)Line Input #1, lineIf Val(Left(line, 4) = Val(n) Theni0 = Val(Mid(line, 7, 4)j0 = Val(Mid(line, 12, 4)k0 = Val(Mid(line, 17, 4)e0 = Val(Mid(line, 23, 4)If Val(Left(line, 4) = Val(Label

22、5) Then e0 = -Val(Mid(line, 23, 4)Elsei = Val(Mid(line, 7, 4)j = Val(Mid(line, 12, 4)k = Val(Mid(line, 17, 4)e = Val(Mid(line, 23, 4)If Val(Left(line, 4) = Val(Label5) Then e = -Val(Mid(line, 23, 4)End Ifdist = Sqr(i - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = Sqr(2) / 2 Thenw = w + e * e0End IfLoop

23、Close #1Next nText4.Text = -2 * wIf Val(Text5.Text) > 0 ThenLabel7.Caption = "flipped"Elseboltzmanfactor = Exp(Val(Text5.Text) / (300 * 1.38 * 10 -23)If Int(1000 * Rnd) > Val(boltzmanfactor) ThenLabel7.Caption = "unchanged"ElseLabel7.Caption = "unflipped"End IfEnd

24、 If Case "文件4"Open "E: 文件4" For Input As #1c = InputBox("输入序号")Do While Not EOF(1)Line Input #1, lineIf Val(Left(line, 4) = Val(c) ThenText1.Text = linei0 = Val(Mid(Text1, 7, 4)j0 = Val(Mid(Text1, 12, 4)k0 = Val(Mid(Text1, 17, 4)If i1 - i0 >= 2 Then i0 = i0 + 2If i1

25、- i0 <= -2 Then i0 = i0 - 2If j1 - j0 >= 2 Then j0 = j0 + 2If j1 - j0 <= -2 Then j0 = j0 - 2If k1 - k0 >= 2 Then k0 = k0 + 2If k1 - k0 <= -2 Then k0 = k0 - 2End IfLoopClose #1 Open "E: 文件4" For Input As #1Do While Not EOF(1)Line Input #1, linei = Val(Mid(line, 7, 4)j = Val(M

26、id(line, 12, 4)k = Val(Mid(line, 17, 4)dist = Sqr(i - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = Sqr(3) / 4 Thennum = num + 1End IfLoopClose #1Text2.Text = numFor n = 1 To 381Open "E: 文件4" For Input As #1Do While Not EOF(1)Line Input #1, lineIf Val(Left(line, 4) = Val(n) Theni0 = Val(Mid(li

27、ne, 7, 4)j0 = Val(Mid(line, 12, 4)k0 = Val(Mid(line, 17, 4)e0 = Val(Mid(line, 23, 4)Elsei = Val(Mid(line, 7, 4)j = Val(Mid(line, 12, 4)k = Val(Mid(line, 17, 4)e = Val(Mid(line, 23, 4)End Ifdist = Sqr(i - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = Sqr(3) / 4 Thenw = w + e * e0End IfLoopClose #1Next nT

28、ext3.Text = -2 * wh = Int(381 * Rnd)Label5.Caption = hFor n = 1 To 381Open "E: 文件4" For Input As #1Do While Not EOF(1)Line Input #1, lineIf Val(Left(line, 4) = Val(n) Theni0 = Val(Mid(line, 7, 4)j0 = Val(Mid(line, 12, 4)k0 = Val(Mid(line, 17, 4)e0 = Val(Mid(line, 23, 4)If Val(Left(line, 4)

29、 = Val(Label5) Then e0 = -Val(Mid(line, 23, 4)Elsei = Val(Mid(line, 7, 4)j = Val(Mid(line, 12, 4)k = Val(Mid(line, 17, 4)e = Val(Mid(line, 23, 4)If Val(Left(line, 4) = Val(Label5) Then e = -Val(Mid(line, 23, 4)End Ifdist = Sqr(i - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = Sqr(3) / 4 Thenw = w + e *

30、e0End IfLoopClose #1Next nText4.Text = -2 * wIf Val(Text5.Text) > 0 ThenLabel7.Caption = "flipped"Elseboltzmanfactor = Exp(Val(Text5.Text) / (300 * 1.38 * 10 -23)If Int(1000 * Rnd) > Val(boltzmanfactor) ThenLabel7.Caption = "unchanged"ElseLabel7.Caption = "unflipped&qu

31、ot;End IfEnd IfEnd SelectEnd SubPrivate Sub command1_Click()Dim i, j, k, l, m As SingleDim a As Stringa = InputBox("输入晶体类型")Select Case aCase "文件1"GoTo 文件1Case "文件2"GoTo 文件2Case "文件3"GoTo 文件3Case "文件4"GoTo文件4End Selectbasc:For k = 0 To 4For j = 0 To

32、4For i = 0 To 4l = l + 1m = Int(2 * Rnd)If m = 0 ThenPrint #1, l; Tab(6); i; Tab(11); j; Tab(16); k; Tab(22); "=" Format(m, "0.0")ElsePrint #1, l; Tab(6); i; Tab(11); j; Tab(16); k; Tab(22); "=" Format(m, "0.0")End IfNext iNext jNext kReturn文件1:Open "E:文件

33、1" For Output As #1 '写入文本文件GoSub basc '调用函数，输入坐标点GoTo finish文件2:Open "E:文件2" For Output As #1GoSub bascFor k = k To 4For j = j To 4For i = i To 4l = l + 1m = Int(2 * Rnd)If m = 0 ThenPrint #1, l; Tab(7); Format(i, "0.0"); Tab(12); Format(j, "0.0"); Tab(17);

34、 Format(k, "0.0"); Tab(22); "=" Format(m, "0.0")ElsePrint #1, l; Tab(7); Format(i, "0.0"); Tab(12); Format(j, "0.0"); Tab(17); Format(k, "0.0"); Tab(22); "=" Format(m, "0.0")End IfNext iNext jNext kGoTo finish文件3:Open &q

35、uot;E:文件3" For Output As #1GoSub basck = 0For k = k To 4For i = i To 4For j = j To 4l = l + 1m = Int(2 * Rnd)If m = 0 ThenPrint #1, l; Tab(7); Format(i, "0.0"); Tab(12); Format(j, "0.0"); Tab(17); Format(k, "0.0"); Tab(22); "=" Format(m, "0.0")E

36、lsePrint #1, l; Tab(7); Format(i, "0.0"); Tab(12); Format(j, "0.0"); Tab(17); Format(k, "0.0"); Tab(22); "=" Format(m, "0.0")End IfNext jNext iNext ki = 0For k = k To 4For j = j To 4For i = i To 4l = l + 1m = Int(2 * Rnd)If m = 0 ThenPrint #1, l; Tab

37、(7); Format(i, "0.0"); Tab(12); Format(j, "0.0"); Tab(17); Format(k, "0.0"); Tab(22); "=" Format(m, "0.0")ElsePrint #1, l; Tab(7); Format(i, "0.0"); Tab(12); Format(j, "0.0"); Tab(17); Format(k, "0.0"); Tab(22); "=&

38、quot; Format(m, "0.0")End IfNext ii = 0Next jNext kj = 0For k = k To 4For j = j To 4For i = i To 4l = l + 1m = Int(2 * Rnd)If m = 0 ThenPrint #1, l; Tab(7); Format(i, "0.0"); Tab(12); Format(j, "0.0"); Tab(17); Format(k, "0.0"); Tab(22); "=" Format(m

39、, "0.0")ElsePrint #1, l; Tab(7); Format(i, "0.0"); Tab(12); Format(j, "0.0"); Tab(17); Format(k, "0.0"); Tab(22); "=" Format(m, "0.0")End IfNext iNext jj = 0Next kGoTo finish文件4:Open "E: 文件4" For Output As #1GoSub bascFor k = k To

40、 4For j = j To 4For i = i To 4l = l + 1m = Int(2 * Rnd)If m = 0 ThenPrint #1, l; Tab(7); Format(i, "0.00"); Tab(12); Format(j, "0.00"); Tab(17); Format(k, "0.00"); Tab(22); "=" Format(m, "0.0")ElsePrint #1, l; Tab(7); Format(i, "0.00"); Tab(12); Format(j, "0.00"); Tab(17); Format(k, "0.00"); Tab(22); "=" Format(m, "0.0")End IfNext iNext jNext kFor k = k To 4For j = j To 4For i = i To 4l = l + 1m