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 * wText5.Text = Text3.Text - Text4.TextIf 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 = "

17、unflipped"End IfEnd IfCase "文件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

18、 Then i0 = i0 + 2If 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 <= -2 Then k0 = k0 - 2End IfLoopClose #1 Open "E: 文件3" For Input As #1Do While Not EOF(1)Line Input #1, linei = Val(M

19、id(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 = 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) = Va

20、l(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(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 * e0E

21、nd IfLoopClose #1Next nText3.Text = -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

22、, 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 = 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

23、(2) / 2 Thenw = w + e * e0End IfLoopClose #1Next nText4.Text = -2 * wText5.Text = Text3.Text - Text4.TextIf 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 =

24、"unchanged"ElseLabel7.Caption = "unflipped"End IfEnd 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, 1

25、2, 4)k0 = Val(Mid(Text1, 17, 4)If i1 - i0 >= 2 Then i0 = i0 + 2If 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 <= -2 Then k0 = k0 - 2End IfLoopClose #1 Open "E: 文件4" For Input As

26、#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 = Sqr(3) / 4 Thennum = num + 1End IfLoopClose #1Text2.Text = numFor n = 1 To 381Open "E: 文件4" For Input As #1Do While Not EO

27、F(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(line, 17, 4)e = Val(Mid(line, 23, 4)End Ifdist = Sqr(i - i0) 2 + (j - j0) 2 + (k

28、- k0) 2)If dist = Sqr(3) / 4 Thenw = w + e * e0End IfLoopClose #1Next nText3.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,

29、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 = 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

30、 - i0) 2 + (j - j0) 2 + (k - k0) 2)If dist = Sqr(3) / 4 Thenw = w + e * e0End IfLoopClose #1Next nText4.Text = -2 * wText5.Text = Text3.Text - Text4.TextIf Val(Text5.Text) > 0 ThenLabel7.Caption = "flipped"Elseboltzmanfactor = Exp(Val(Text5.Text) / (300 * 1.38 * 10 -23)If Int(1000 * Rnd

31、) > Val(boltzmanfactor) ThenLabel7.Caption = "unchanged"ElseLabel7.Caption = "unflipped"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&qu

32、ot;GoTo 文件2Case "文件3"GoTo 文件3Case "文件4"GoTo文件4End Selectbasc:For k = 0 To 4For j = 0 To 4For i = 0 To 4l = l + 1m = Int(2 * Rnd)If m = 0 Thenm = -0.5Print #1, l; Tab(6); i; Tab(11); j; Tab(16); k; Tab(22); "=" Format(m, "0.0")Elsem = 0.5Print #1, l; Tab(6); i;

33、 Tab(11); j; Tab(16); k; Tab(22); "=" Format(m, "0.0")End IfNext iNext jNext kReturn文件1:Open "E:文件1" For Output As #1 '写入文本文件GoSub basc '调用函数，输入坐标点GoTo finish文件2:Open "E:文件2" For Output As #1GoSub basck = 0.5j = 0.5i = 0.5For k = k To 4For j = j To 4Fo

34、r i = i To 4l = l + 1m = Int(2 * Rnd)If m = 0 Thenm = -0.5Print #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")Elsem = 0.5Print #1, l; Tab(7); Format(i, "0.0"); Tab(12); F

35、ormat(j, "0.0"); Tab(17); Format(k, "0.0"); Tab(22); "=" Format(m, "0.0")End IfNext ii = 0.5Next jj = 0.5Next kGoTo finish文件3:Open "E:文件3" For Output As #1GoSub basck = 0j = 0.5i = 0.5For k = k To 4For i = i To 4For j = j To 4l = l + 1m = Int(2 * Rnd

36、)If m = 0 Thenm = -0.5Print #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")Elsem = 0.5Print #1, l; Tab(7); Format(i, "0.0"); Tab(12); Format(j, "0.0"); Tab(17);

37、Format(k, "0.0"); Tab(22); "=" Format(m, "0.0")End IfNext jj = 0.5Next ii = 0.5Next kk = 0.5j = 0.5i = 0For k = k To 4For j = j To 4For i = i To 4l = l + 1m = Int(2 * Rnd)If m = 0 Thenm = -0.5Print #1, l; Tab(7); Format(i, "0.0"); Tab(12); Format(j, "0.0&

38、quot;); Tab(17); Format(k, "0.0"); Tab(22); "=" Format(m, "0.0")Elsem = 0.5Print #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 ii = 0Next j

39、j = 0.5Next kk = 0.5j = 0i = 0.5For k = k To 4For j = j To 4For i = i To 4l = l + 1m = Int(2 * Rnd)If m = 0 Thenm = -0.5Print #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")Elsem =

40、0.5Print #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 ii = 0.5Next jj = 0Next kGoTo finish文件4:Open "E: 文件4" For Output As #1GoSub basck = 0.25j = 0.25i = 0.25

41、For k = k To 4For j = j To 4For i = i To 4l = l + 1m = Int(2 * Rnd)If m = 0 Thenm = -0.5Print #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")Elsem = 0.5Print #1, l; Tab(7); Forma

42、t(i, "0.00"); Tab(12); Format(j, "0.00"); Tab(17); Format(k, "0.00"); Tab(22); "=" Format(m, "0.0")End IfNext ii = 0.25Next jj = 0.25Next kk = 0.75j = 0.75i = 0.25For k = k To 4For j = j To 4For i = i To 4l = l + 1m = Int(2 * Rnd)If m = 0 Thenm = -0.5Print #1, l; Tab(7); Format(i, "0.00"); Tab(