' Generates Most Perfect Magic Squares of order 12, magic sum 870
' 1/2 Rows and 1/2 Columns sum to 435
' Tested with Office 365 under Windows 10
Sub MgcSqr12f()
Dim a1(144), a(144), b1(144), b(144), c(144)
y = MsgBox("Locked", vbCritical, "Routine MgcSqr12f")
End
n2 = 0: n3 = 0: n9 = 0: n10 = 0: k1 = 1: k2 = 1
Sheets("Klad1").Select
t1 = Timer
' Define variables
GoSub 2030 ' Define Natural Numbers
' Generate Squares
For j144 = m1 To m2 'a(144)
If b(a1(j144)) = 0 Then b(a1(j144)) = a1(j144): c(144) = a1(j144) Else GoTo 1440
a(144) = a1(j144)
a(66) = p12 - a(144): If b(a(66)) = 0 Then b(a(66)) = a(66): c(66) = a(66) Else GoTo 660
For j143 = m1 To m2 'a(143)
If b(a1(j43)) = 0 Then b(a1(j143)) = a1(j143): c(143) = a1(j143) Else GoTo 1430
a(143) = a1(j143)
a(65) = p12 - a(143): If b(a(65)) = 0 Then b(a(65)) = a(65): c(65) = a(65) Else GoTo 650
For j142 = m1 To m2 'a(142)
If b(a1(j142)) = 0 Then b(a1(j142)) = a1(j142): c(142) = a1(j142) Else GoTo 1420
a(142) = a1(j142)
a(64) = p12 - a(142): If b(a(64)) = 0 Then b(a(64)) = a(64): c(64) = a(64) Else GoTo 640
For j141 = m1 To m2 'a(141)
If b(a1(j141)) = 0 Then b(a1(j141)) = a1(j141): c(141) = a1(j141) Else GoTo 1410
a(141) = a1(j141)
a(63) = p12 - a(141): If b(a(63)) = 0 Then b(a(63)) = a(63): c(63) = a(63) Else GoTo 630
For j140 = m1 To m2 'a(140)
If b(a1(j140)) = 0 Then b(a1(j140)) = a1(j140): c(140) = a1(j140) Else GoTo 1400
a(140) = a1(j140)
a(62) = p12 - a(140): If b(a(62)) = 0 Then b(a(62)) = a(62): c(62) = a(62) Else GoTo 620
a(139) = s1 / 2 - a(140) - a(141) - a(142) - a(143) - a(144)
If a(139) < a1(m1) Or a(139) > a1(m2) Then GoTo 1390
If b1(a(139)) = 0 Then GoTo 1390
If b(a(139)) = 0 Then b(a(139)) = a(139): c(139) = a(139) Else GoTo 1390
a(61) = p12 - a(139): If b(a(61)) = 0 Then b(a(61)) = a(61): c(61) = a(61) Else GoTo 610
For j138 = m1 To m2 'a(138)
If b(a1(j138)) = 0 Then b(a1(j138)) = a1(j138): c(138) = a1(j138) Else GoTo 1380
a(138) = a1(j138)
a(72) = p12 - a(138): If b(a(72)) = 0 Then b(a(72)) = a(72): c(72) = a(72) Else GoTo 720
a(137) = s1 / 3 - a(138) - a(143) - a(144): If a(137) < a1(m1) Or a(137) > a1(m2) Then GoTo 1370
If b1(a(137)) = 0 Then GoTo 1370
If b(a(137)) = 0 Then b(a(137)) = a(137): c(137) = a(137) Else GoTo 1370
a(71) = p12 - a(137): If b(a(71)) = 0 Then b(a(71)) = a(71): c(71) = a(71) Else GoTo 710
a(136) = a(138) - a(142) + a(144): If a(136) < a1(m1) Or a(136) > a1(m2) Then GoTo 1360
If b1(a(136)) = 0 Then GoTo 1360
If b(a(136)) = 0 Then b(a(136)) = a(136): c(136) = a(136) Else GoTo 1360
a(70) = p12 - a(136): If b(a(70)) = 0 Then b(a(70)) = a(70): c(70) = a(70) Else GoTo 700
a(135) = s1 / 3 - a(138) - a(141) - a(144): If a(135) < a1(m1) Or a(135) > a1(m2) Then GoTo 1350
If b1(a(135)) = 0 Then GoTo 1350
If b(a(135)) = 0 Then b(a(135)) = a(135): c(135) = a(135) Else GoTo 1350
a(69) = p12 - a(135): If b(a(69)) = 0 Then b(a(69)) = a(69): c(69) = a(69) Else GoTo 690
a(134) = a(138) - a(140) + a(144): If a(134) < a1(m1) Or a(134) > a1(m2) Then GoTo 1340
If b1(a(134)) = 0 Then GoTo 1340
If b(a(134)) = 0 Then b(a(134)) = a(134): c(134) = a(134) Else GoTo 1340
a(68) = p12 - a(134): If b(a(68)) = 0 Then b(a(68)) = a(68): c(68) = a(68) Else GoTo 680
a(133) = -s1 / 6 - a(138) + a(140) + a(141) + a(142) + a(143)
If a(133) < a1(m1) Or a(133) > a1(m2) Then GoTo 1330
If b1(a(133)) = 0 Then GoTo 1330
If b(a(133)) = 0 Then b(a(133)) = a(133): c(133) = a(133) Else GoTo 1330
a(67) = p12 - a(133): If b(a(67)) = 0 Then b(a(67)) = a(67): c(67) = a(67) Else GoTo 670
For j132 = m1 To m2 'a(132)
If b(a1(j132)) = 0 Then b(a1(j132)) = a1(j132): c(132) = a1(j132) Else GoTo 1320
a(132) = a1(j132)
a(54) = p12 - a(132): If b(a(54)) = 0 Then b(a(54)) = a(54): c(54) = a(54) Else GoTo 540
a(131) = s1 / 3 - a(132) - a(143) - a(144): If a(131) < a1(m1) Or a(131) > a1(m2) Then GoTo 1310
If b1(a(131)) = 0 Then GoTo 1310
If b(a(131)) = 0 Then b(a(131)) = a(131): c(131) = a(131) Else GoTo 1310
a(53) = p12 - a(131): If b(a(53)) = 0 Then b(a(53)) = a(53): c(53) = a(53) Else GoTo 530
a(130) = a(132) - a(142) + a(144)
If a(130) < a1(m1) Or a(130) > a1(m2) Then GoTo 1300
If b1(a(130)) = 0 Then GoTo 1300
If b(a(130)) = 0 Then b(a(130)) = a(130): c(130) = a(130) Else GoTo 1300
a(52) = p12 - a(130): If b(a(52)) = 0 Then b(a(52)) = a(52): c(52) = a(52) Else GoTo 520
a(129) = s1 / 3 - a(132) - a(141) - a(144): If a(129) < a1(m1) Or a(129) > a1(m2) Then GoTo 1290
If b1(a(129)) = 0 Then GoTo 1290
If b(a(129)) = 0 Then b(a(129)) = a(129): c(129) = a(129) Else GoTo 1290
a(51) = p12 - a(129): If b(a(51)) = 0 Then b(a(51)) = a(51): c(51) = a(51) Else GoTo 510
a(128) = a(132) - a(140) + a(144)
If a(128) < a1(m1) Or a(128) > a1(m2) Then GoTo 1280
If b1(a(128)) = 0 Then GoTo 1280
If b(a(128)) = 0 Then b(a(128)) = a(128): c(128) = a(128) Else GoTo 1280
a(50) = p12 - a(128): If b(a(50)) = 0 Then b(a(50)) = a(50): c(50) = a(50) Else GoTo 500
a(127) = -s1 / 6 - a(132) + a(140) + a(141) + a(142) + a(143)
If a(127) < a1(m1) Or a(127) > a1(m2) Then GoTo 1270
If b1(a(127)) = 0 Then GoTo 1270
If b(a(127)) = 0 Then b(a(127)) = a(127): c(127) = a(127) Else GoTo 1270
a(49) = p12 - a(127): If b(a(49)) = 0 Then b(a(49)) = a(49): c(49) = a(49) Else GoTo 490
a(126) = a(132) - a(138) + a(144): If a(126) < a1(m1) Or a(126) > a1(m2) Then GoTo 1260
If b1(a(126)) = 0 Then GoTo 1260
If b(a(126)) = 0 Then b(a(126)) = a(126): c(126) = a(126) Else GoTo 1260
a(60) = p12 - a(126): If b(a(60)) = 0 Then b(a(60)) = a(60): c(60) = a(60) Else GoTo 600
a(125) = -a(132) + a(138) + a(143): If a(125) < a1(m1) Or a(125) > a1(m2) Then GoTo 1250
If b1(a(125)) = 0 Then GoTo 1250
If b(a(125)) = 0 Then b(a(125)) = a(125): c(125) = a(125) Else GoTo 1250
a(59) = p12 - a(125): If b(a(59)) = 0 Then b(a(59)) = a(59): c(59) = a(59) Else GoTo 590
a(124) = a(132) - a(138) + a(142): If a(124) < a1(m1) Or a(124) > a1(m2) Then GoTo 1240
If b1(a(124)) = 0 Then GoTo 1240
If b(a(124)) = 0 Then b(a(124)) = a(124): c(124) = a(124) Else GoTo 1240
a(58) = p12 - a(124): If b(a(58)) = 0 Then b(a(58)) = a(58): c(58) = a(58) Else GoTo 580
a(123) = -a(132) + a(138) + a(141): If a(123) < a1(m1) Or a(123) > a1(m2) Then GoTo 1230
If b1(a(123)) = 0 Then GoTo 1230
If b(a(123)) = 0 Then b(a(123)) = a(123): c(123) = a(123) Else GoTo 1230
a(57) = p12 - a(123): If b(a(57)) = 0 Then b(a(57)) = a(57): c(57) = a(57) Else GoTo 570
a(122) = a(132) - a(138) + a(140): If a(122) < a1(m1) Or a(122) > a1(m2) Then GoTo 1220
If b1(a(122)) = 0 Then GoTo 1220
If b(a(122)) = 0 Then b(a(122)) = a(122): c(122) = a(122) Else GoTo 1220
a(56) = p12 - a(122): If b(a(56)) = 0 Then b(a(56)) = a(56): c(56) = a(56) Else GoTo 560
a(121) = s1 / 2 - a(132) + a(138) - a(140) - a(141) - a(142) - a(143) - a(144)
If a(121) < a1(m1) Or a(121) > a1(m2) Then GoTo 1210
If b1(a(121)) = 0 Then GoTo 1210
If b(a(121)) = 0 Then b(a(121)) = a(121): c(121) = a(121) Else GoTo 1210
a(55) = p12 - a(121): If b(a(55)) = 0 Then b(a(55)) = a(55): c(55) = a(55) Else GoTo 550
For j120 = m1 To m2 'a(120)
If b(a1(j120)) = 0 Then b(a1(j120)) = a1(j120): c(120) = a1(j120) Else GoTo 1200
a(120) = a1(j120)
a(42) = p12 - a(120): If b(a(42)) = 0 Then b(a(42)) = a(42): c(42) = a(42) Else GoTo 420
a(119) = -a(120) + a(143) + a(144): If a(119) < a1(m1) Or a(119) > a1(m2) Then GoTo 1190
If b1(a(119)) = 0 Then GoTo 1190
If b(a(119)) = 0 Then b(a(119)) = a(119): c(119) = a(119) Else GoTo 1190
a(41) = p12 - a(119): If b(a(41)) = 0 Then b(a(41)) = a(41): c(41) = a(41) Else GoTo 410
a(118) = a(120) + a(142) - a(144): If a(118) < a1(m1) Or a(118) > a1(m2) Then GoTo 1180
If b1(a(118)) = 0 Then GoTo 1180
If b(a(118)) = 0 Then b(a(118)) = a(118): c(118) = a(118) Else GoTo 1180
a(40) = p12 - a(118): If b(a(40)) = 0 Then b(a(40)) = a(40): c(40) = a(40) Else GoTo 400
a(117) = -a(120) + a(141) + a(144): If a(117) < a1(m1) Or a(117) > a1(m2) Then GoTo 1170
If b1(a(117)) = 0 Then GoTo 1170
If b(a(117)) = 0 Then b(a(117)) = a(117): c(117) = a(117) Else GoTo 1170
a(39) = p12 - a(117): If b(a(39)) = 0 Then b(a(39)) = a(39): c(39) = a(39) Else GoTo 390
a(116) = a(120) + a(140) - a(144): If a(116) < a1(m1) Or a(116) > a1(m2) Then GoTo 1160
If b1(a(116)) = 0 Then GoTo 1160
If b(a(116)) = 0 Then b(a(116)) = a(116): c(116) = a(116) Else GoTo 1160
a(38) = p12 - a(116): If b(a(38)) = 0 Then b(a(38)) = a(38): c(38) = a(38) Else GoTo 380
a(115) = s1 / 2 - a(120) - a(140) - a(141) - a(142) - a(143)
If a(115) < a1(m1) Or a(115) > a1(m2) Then GoTo 1150
If b1(a(115)) = 0 Then GoTo 1150
If b(a(115)) = 0 Then b(a(115)) = a(115): c(115) = a(115) Else GoTo 1150
a(37) = p12 - a(115): If b(a(37)) = 0 Then b(a(37)) = a(37): c(37) = a(37) Else GoTo 370
a(114) = a(120) + a(138) - a(144): If a(114) < a1(m1) Or a(114) > a1(m2) Then GoTo 1140
If b1(a(114)) = 0 Then GoTo 1140
If b(a(114)) = 0 Then b(a(114)) = a(114): c(114) = a(114) Else GoTo 1140
a(48) = p12 - a(114): If b(a(48)) = 0 Then b(a(48)) = a(48): c(48) = a(48) Else GoTo 480
a(113) = s1 / 3 - a(120) - a(138) - a(143): If a(113) < a1(m1) Or a(113) > a1(m2) Then GoTo 1130
If b1(a(113)) = 0 Then GoTo 1130
If b(a(113)) = 0 Then b(a(113)) = a(113): c(113) = a(113) Else GoTo 1130
a(47) = p12 - a(113): If b(a(47)) = 0 Then b(a(47)) = a(47): c(47) = a(47) Else GoTo 470
a(112) = a(120) + a(138) - a(142): If a(112) < a1(m1) Or a(112) > a1(m2) Then GoTo 1120
If b1(a(112)) = 0 Then GoTo 1120
If b(a(112)) = 0 Then b(a(112)) = a(112): c(112) = a(112) Else GoTo 1120
a(46) = p12 - a(112): If b(a(46)) = 0 Then b(a(46)) = a(46): c(46) = a(46) Else GoTo 460
a(111) = s1 / 3 - a(120) - a(138) - a(141): If a(111) < a1(m1) Or a(111) > a1(m2) Then GoTo 1110
If b1(a(111)) = 0 Then GoTo 1110
If b(a(111)) = 0 Then b(a(111)) = a(111): c(111) = a(111) Else GoTo 1110
a(45) = p12 - a(111): If b(a(45)) = 0 Then b(a(45)) = a(45): c(45) = a(45) Else GoTo 450
a(110) = a(120) + a(138) - a(140): If a(110) < a1(m1) Or a(110) > a1(m2) Then GoTo 1100
If b1(a(110)) = 0 Then GoTo 1100
If b(a(110)) = 0 Then b(a(110)) = a(110): c(110) = a(110) Else GoTo 1100
a(44) = p12 - a(110): If b(a(44)) = 0 Then b(a(44)) = a(44): c(44) = a(44) Else GoTo 440
a(109) = -s1 / 6 - a(120) - a(138) + a(140) + a(141) + a(142) + a(143) + a(144)
If a(109) < a1(m1) Or a(109) > a1(m2) Then GoTo 1090
If b1(a(109)) = 0 Then GoTo 1090
If b(a(109)) = 0 Then b(a(109)) = a(109): c(109) = a(109) Else GoTo 1090
a(43) = p12 - a(109): If b(a(43)) = 0 Then b(a(43)) = a(43): c(43) = a(43) Else GoTo 430
For j108 = m1 To m2 'a(108)
If b(a1(j108)) = 0 Then b(a1(j108)) = a1(j108): c(108) = a1(j108) Else GoTo 1080
a(108) = a1(j108)
a(30) = p12 - a(108): If b(a(30)) = 0 Then b(a(30)) = a(30): c(30) = a(30) Else GoTo 300
a(107) = s1 / 3 - a(108) - a(143) - a(144): If a(107) < a1(m1) Or a(107) > a1(m2) Then GoTo 1070
If b1(a(107)) = 0 Then GoTo 1070
If b(a(107)) = 0 Then b(a(107)) = a(107): c(107) = a(107) Else GoTo 1070
a(29) = p12 - a(107): If b(a(29)) = 0 Then b(a(29)) = a(29): c(29) = a(29) Else GoTo 290
a(106) = a(108) - a(142) + a(144): If a(106) < a1(m1) Or a(106) > a1(m2) Then GoTo 1060
If b1(a(106)) = 0 Then GoTo 1060
If b(a(106)) = 0 Then b(a(106)) = a(106): c(106) = a(106) Else GoTo 1060
a(28) = p12 - a(106): If b(a(28)) = 0 Then b(a(28)) = a(28): c(28) = a(28) Else GoTo 280
a(105) = s1 / 3 - a(108) - a(141) - a(144): If a(105) < a1(m1) Or a(105) > a1(m2) Then GoTo 1050
If b1(a(105)) = 0 Then GoTo 1050
If b(a(105)) = 0 Then b(a(105)) = a(105): c(105) = a(105) Else GoTo 1050
a(27) = p12 - a(105): If b(a(27)) = 0 Then b(a(27)) = a(27): c(27) = a(27) Else GoTo 270
a(104) = a(108) - a(140) + a(144): If a(104) < a1(m1) Or a(104) > a1(m2) Then GoTo 1040
If b1(a(104)) = 0 Then GoTo 1040
If b(a(104)) = 0 Then b(a(104)) = a(104): c(104) = a(104) Else GoTo 1040
a(26) = p12 - a(104): If b(a(26)) = 0 Then b(a(26)) = a(26): c(26) = a(26) Else GoTo 260
a(103) = -s1 / 6 - a(108) + a(140) + a(141) + a(142) + a(143)
If a(103) < a1(m1) Or a(103) > a1(m2) Then GoTo 1030
If b1(a(103)) = 0 Then GoTo 1030
If b(a(103)) = 0 Then b(a(103)) = a(103): c(103) = a(103) Else GoTo 1030
a(25) = p12 - a(103): If b(a(25)) = 0 Then b(a(25)) = a(25): c(25) = a(25) Else GoTo 250
a(102) = a(108) - a(138) + a(144): If a(102) < a1(m1) Or a(102) > a1(m2) Then GoTo 1020
If b1(a(102)) = 0 Then GoTo 1020
If b(a(102)) = 0 Then b(a(102)) = a(102): c(102) = a(102) Else GoTo 1020
a(36) = p12 - a(102): If b(a(36)) = 0 Then b(a(36)) = a(36): c(36) = a(36) Else GoTo 360
a(101) = -a(108) + a(138) + a(143): If a(101) < a1(m1) Or a(101) > a1(m2) Then GoTo 1010
If b1(a(101)) = 0 Then GoTo 1010
If b(a(101)) = 0 Then b(a(101)) = a(101): c(101) = a(101) Else GoTo 1010
a(35) = p12 - a(101): If b(a(35)) = 0 Then b(a(35)) = a(35): c(35) = a(35) Else GoTo 350
a(100) = a(108) - a(138) + a(142): If a(100) < a1(m1) Or a(100) > a1(m2) Then GoTo 1000
If b1(a(100)) = 0 Then GoTo 1000
If b(a(100)) = 0 Then b(a(100)) = a(100): c(100) = a(100) Else GoTo 1000
a(34) = p12 - a(100): If b(a(34)) = 0 Then b(a(34)) = a(34): c(34) = a(34) Else GoTo 340
a(99) = -a(108) + a(138) + a(141): If a(99) < a1(m1) Or a(99) > a1(m2) Then GoTo 990
If b1(a(99)) = 0 Then GoTo 990
If b(a(99)) = 0 Then b(a(99)) = a(99): c(99) = a(99) Else GoTo 990
a(33) = p12 - a(99): If b(a(33)) = 0 Then b(a(33)) = a(33): c(33) = a(33) Else GoTo 330
a(98) = a(108) - a(138) + a(140): If a(98) < a1(m1) Or a(98) > a1(m2) Then GoTo 980
If b1(a(98)) = 0 Then GoTo 980
If b(a(98)) = 0 Then b(a(98)) = a(98): c(98) = a(98) Else GoTo 980
a(32) = p12 - a(98): If b(a(32)) = 0 Then b(a(32)) = a(32): c(32) = a(32) Else GoTo 320
a(97) = s1 / 2 - a(108) + a(138) - a(140) - a(141) - a(142) - a(143) - a(144)
If a(97) < a1(m1) Or a(97) > a1(m2) Then GoTo 970
If b1(a(97)) = 0 Then GoTo 970
If b(a(97)) = 0 Then b(a(97)) = a(97): c(97) = a(97) Else GoTo 970
a(31) = p12 - a(97): If b(a(31)) = 0 Then b(a(31)) = a(31): c(31) = a(31) Else GoTo 310
For j96 = m1 To m2 'a(96)
If b(a1(j96)) = 0 Then b(a1(j96)) = a1(j96): c(96) = a1(j96) Else GoTo 960
a(96) = a1(j96)
a(18) = p12 - a(96): If b(a(18)) = 0 Then b(a(18)) = a(18): c(18) = a(18) Else GoTo 180
a(95) = -a(96) + a(143) + a(144): If a(95) < a1(m1) Or a(95) > a1(m2) Then GoTo 950
If b1(a(95)) = 0 Then GoTo 950
If b(a(95)) = 0 Then b(a(95)) = a(95): c(95) = a(95) Else GoTo 950
a(17) = p12 - a(95): If b(a(17)) = 0 Then b(a(17)) = a(17): c(17) = a(17) Else GoTo 170
a(94) = a(96) + a(142) - a(144): If a(94) < a1(m1) Or a(94) > a1(m2) Then GoTo 940
If b1(a(94)) = 0 Then GoTo 940
If b(a(94)) = 0 Then b(a(94)) = a(94): c(94) = a(94) Else GoTo 940
a(16) = p12 - a(94): If b(a(16)) = 0 Then b(a(16)) = a(16): c(16) = a(16) Else GoTo 160
a(93) = -a(96) + a(141) + a(144): If a(93) < a1(m1) Or a(93) > a1(m2) Then GoTo 930
If b1(a(93)) = 0 Then GoTo 930
If b(a(93)) = 0 Then b(a(93)) = a(93): c(93) = a(93) Else GoTo 930
a(15) = p12 - a(93): If b(a(15)) = 0 Then b(a(15)) = a(15): c(15) = a(15) Else GoTo 150
a(92) = a(96) + a(140) - a(144): If a(92) < a1(m1) Or a(92) > a1(m2) Then GoTo 920
If b1(a(92)) = 0 Then GoTo 920
If b(a(92)) = 0 Then b(a(92)) = a(92): c(92) = a(92) Else GoTo 920
a(14) = p12 - a(92): If b(a(14)) = 0 Then b(a(14)) = a(14): c(14) = a(14) Else GoTo 140
a(91) = s1 / 2 - a(96) - a(140) - a(141) - a(142) - a(143)
If a(91) < a1(m1) Or a(91) > a1(m2) Then GoTo 910
If b1(a(91)) = 0 Then GoTo 910
If b(a(91)) = 0 Then b(a(91)) = a(91): c(91) = a(91) Else GoTo 910
a(13) = p12 - a(91): If b(a(13)) = 0 Then b(a(13)) = a(13): c(13) = a(13) Else GoTo 130
a(90) = a(96) + a(138) - a(144): If a(90) < a1(m1) Or a(90) > a1(m2) Then GoTo 900
If b1(a(90)) = 0 Then GoTo 900
If b(a(90)) = 0 Then b(a(90)) = a(90): c(90) = a(90) Else GoTo 900
a(24) = p12 - a(90): If b(a(24)) = 0 Then b(a(24)) = a(24): c(24) = a(24) Else GoTo 240
a(89) = s1 / 3 - a(96) - a(138) - a(143): If a(89) < a1(m1) Or a(89) > a1(m2) Then GoTo 890
If b1(a(89)) = 0 Then GoTo 890
If b(a(89)) = 0 Then b(a(89)) = a(89): c(89) = a(89) Else GoTo 890
a(23) = p12 - a(89): If b(a(23)) = 0 Then b(a(23)) = a(23): c(23) = a(23) Else GoTo 230
a(88) = a(96) + a(138) - a(142): If a(88) < a1(m1) Or a(88) > a1(m2) Then GoTo 880
If b1(a(88)) = 0 Then GoTo 880
If b(a(88)) = 0 Then b(a(88)) = a(88): c(88) = a(88) Else GoTo 880
a(22) = p12 - a(88): If b(a(22)) = 0 Then b(a(22)) = a(22): c(22) = a(22) Else GoTo 220
a(87) = s1 / 3 - a(96) - a(138) - a(141): If a(87) < a1(m1) Or a(87) > a1(m2) Then GoTo 870
If b1(a(87)) = 0 Then GoTo 870
If b(a(87)) = 0 Then b(a(87)) = a(87): c(87) = a(87) Else GoTo 870
a(21) = p12 - a(87): If b(a(21)) = 0 Then b(a(21)) = a(21): c(21) = a(21) Else GoTo 210
a(86) = a(96) + a(138) - a(140): If a(86) < a1(m1) Or a(86) > a1(m2) Then GoTo 860
If b1(a(86)) = 0 Then GoTo 860
If b(a(86)) = 0 Then b(a(86)) = a(86): c(86) = a(86) Else GoTo 860
a(20) = p12 - a(86): If b(a(20)) = 0 Then b(a(20)) = a(20): c(20) = a(20) Else GoTo 200
a(85) = -s1 / 6 - a(96) - a(138) + a(140) + a(141) + a(142) + a(143) + a(144)
If a(85) < a1(m1) Or a(85) > a1(m2) Then GoTo 850
If b1(a(85)) = 0 Then GoTo 850:
If b(a(85)) = 0 Then b(a(85)) = a(85): c(85) = a(85) Else GoTo 850
a(19) = p12 - a(85): If b(a(19)) = 0 Then b(a(19)) = a(19): c(19) = a(19) Else GoTo 190
a(84) = s1 / 2 - a(96) - a(108) - a(120) - a(132) - a(144)
If a(84) < a1(m1) Or a(84) > a1(m2) Then GoTo 840
If b1(a(84)) = 0 Then GoTo 840
If b(a(84)) = 0 Then b(a(84)) = a(84): c(84) = a(84) Else GoTo 840
a(6) = p12 - a(84): If b(a(6)) = 0 Then b(a(6)) = a(6): c(6) = a(6) Else GoTo 60
a(83) = -s1 / 6 + a(96) + a(108) + a(120) + a(132) - a(143)
If a(83) < a1(m1) Or a(83) > a1(m2) Then GoTo 830
If b1(a(83)) = 0 Then GoTo 830
If b(a(83)) = 0 Then b(a(83)) = a(83): c(83) = a(83) Else GoTo 830
a(5) = p12 - a(83): If b(a(5)) = 0 Then b(a(5)) = a(5): c(5) = a(5) Else GoTo 50
a(82) = s1 / 2 - a(96) - a(108) - a(120) - a(132) - a(142)
If a(82) < a1(m1) Or a(82) > a1(m2) Then GoTo 820
If b1(a(82)) = 0 Then GoTo 820
If b(a(82)) = 0 Then b(a(82)) = a(82): c(82) = a(82) Else GoTo 820
a(4) = p12 - a(82): If b(a(4)) = 0 Then b(a(4)) = a(4): c(4) = a(4) Else GoTo 40
a(81) = -s1 / 6 + a(96) + a(108) + a(120) + a(132) - a(141)
If a(81) < a1(m1) Or a(81) > a1(m2) Then GoTo 810
If b1(a(81)) = 0 Then GoTo 810
If b(a(81)) = 0 Then b(a(81)) = a(81): c(81) = a(81) Else GoTo 810
a(3) = p12 - a(81): If b(a(3)) = 0 Then b(a(3)) = a(3): c(3) = a(3) Else GoTo 30
a(80) = s1 / 2 - a(96) - a(108) - a(120) - a(132) - a(140)
If a(80) < a1(m1) Or a(80) > a1(m2) Then GoTo 800
If b1(a(80)) = 0 Then GoTo 800
If b(a(80)) = 0 Then b(a(80)) = a(80): c(80) = a(80) Else GoTo 800
a(2) = p12 - a(80): If b(a(2)) = 0 Then b(a(2)) = a(2): c(2) = a(2) Else GoTo 20
a(79) = -2 * s1 / 3 + a(96) + a(108) + a(120) + a(132) + a(140) + a(141) + a(142) + a(143) + a(144)
If a(79) < a1(m1) Or a(79) > a1(m2) Then GoTo 790
If b1(a(79)) = 0 Then GoTo 790
If b(a(79)) = 0 Then b(a(79)) = a(79): c(79) = a(79) Else GoTo 790
a(1) = p12 - a(79): If b(a(1)) = 0 Then b(a(1)) = a(1): c(1) = a(1) Else GoTo 10
a(78) = s1 / 2 - a(96) - a(108) - a(120) - a(132) - a(138)
If a(78) < a1(m1) Or a(78) > a1(m2) Then GoTo 780
If b1(a(78)) = 0 Then GoTo 780
If b(a(78)) = 0 Then b(a(78)) = a(78): c(78) = a(78) Else GoTo 780
a(12) = p12 - a(78): If b(a(12)) = 0 Then b(a(12)) = a(12): c(12) = a(12) Else GoTo 120
a(77) = -s1 / 2 + a(96) + a(108) + a(120) + a(132) + a(138) + a(143) + a(144)
If a(77) < a1(m1) Or a(77) > a1(m2) Then GoTo 770
If b1(a(77)) = 0 Then GoTo 770
If b(a(77)) = 0 Then b(a(77)) = a(77): c(77) = a(77) Else GoTo 770
a(11) = p12 - a(77): If b(a(11)) = 0 Then b(a(11)) = a(11): c(11) = a(11) Else GoTo 110
a(76) = s1 / 2 - a(96) - a(108) - a(120) - a(132) - a(138) + a(142) - a(144)
If a(76) < a1(m1) Or a(76) > a1(m2) Then GoTo 760
If b1(a(76)) = 0 Then GoTo 760
If b(a(76)) = 0 Then b(a(76)) = a(76): c(76) = a(76) Else GoTo 760
a(10) = p12 - a(76): If b(a(10)) = 0 Then b(a(10)) = a(10): c(10) = a(10) Else GoTo 100
a(75) = -s1 / 2 + a(96) + a(108) + a(120) + a(132) + a(138) + a(141) + a(144)
If a(75) < a1(m1) Or a(75) > a1(m2) Then GoTo 750
If b1(a(75)) = 0 Then GoTo 750
If b(a(75)) = 0 Then b(a(75)) = a(75): c(75) = a(75) Else GoTo 750
a(9) = p12 - a(75): If b(a(9)) = 0 Then b(a(9)) = a(9): c(9) = a(9) Else GoTo 90
a(74) = s1 / 2 - a(96) - a(108) - a(120) - a(132) - a(138) + a(140) - a(144)
If a(74) < a1(m1) Or a(74) > a1(m2) Then GoTo 740
If b1(a(74)) = 0 Then GoTo 740
If b(a(74)) = 0 Then b(a(74)) = a(74): c(74) = a(74) Else GoTo 740
a(8) = p12 - a(74): If b(a(8)) = 0 Then b(a(8)) = a(8): c(8) = a(8) Else GoTo 80
a(73) = a(96) + a(108) + a(120) + a(132) + a(138) - a(140) - a(141) - a(142) - a(143)
If a(73) < a1(m1) Or a(73) > a1(m2) Then GoTo 730
If b1(a(73)) = 0 Then GoTo 730
If b(a(73)) = 0 Then b(a(73)) = a(73): c(73) = a(73) Else GoTo 730
a(7) = p12 - a(73): If b(a(7)) = 0 Then b(a(7)) = a(7): c(7) = a(7) Else GoTo 70
n9 = n9 + 1 'Total
' GoSub 2645 'Print results (Lines)
GoSub 2650 'Print results (Squares)
' Cells(2, 1).Value = n9 'Counting
b(c(7)) = 0: c(7) = 0
70 b(c(73)) = 0: c(73) = 0
730 b(c(8)) = 0: c(8) = 0
80 b(c(74)) = 0: c(74) = 0
740 b(c(9)) = 0: c(9) = 0
90 b(c(75)) = 0: c(75) = 0
750 b(c(10)) = 0: c(10) = 0
100 b(c(76)) = 0: c(76) = 0
760 b(c(11)) = 0: c(11) = 0
110 b(c(77)) = 0: c(77) = 0
770 b(c(12)) = 0: c(12) = 0
120 b(c(78)) = 0: c(78) = 0
780 b(c(1)) = 0: c(1) = 0
10 b(c(79)) = 0: c(79) = 0
790 b(c(2)) = 0: c(2) = 0
20 b(c(80)) = 0: c(80) = 0
800 b(c(3)) = 0: c(3) = 0
30 b(c(81)) = 0: c(81) = 0
810 b(c(4)) = 0: c(4) = 0
40 b(c(82)) = 0: c(82) = 0
820 b(c(5)) = 0: c(5) = 0
50 b(c(83)) = 0: c(83) = 0
830 b(c(6)) = 0: c(6) = 0
60 b(c(84)) = 0: c(84) = 0
840 b(c(19)) = 0: c(19) = 0
190 b(c(85)) = 0: c(85) = 0
850 b(c(20)) = 0: c(20) = 0
200 b(c(86)) = 0: c(86) = 0
860 b(c(21)) = 0: c(21) = 0
210 b(c(87)) = 0: c(87) = 0
870 b(c(22)) = 0: c(22) = 0
220 b(c(88)) = 0: c(88) = 0
880 b(c(23)) = 0: c(23) = 0
230 b(c(89)) = 0: c(89) = 0
890 b(c(24)) = 0: c(24) = 0
240 b(c(90)) = 0: c(90) = 0
900 b(c(13)) = 0: c(13) = 0
130 b(c(91)) = 0: c(91) = 0
910 b(c(14)) = 0: c(14) = 0
140 b(c(92)) = 0: c(92) = 0
920 b(c(15)) = 0: c(15) = 0
150 b(c(93)) = 0: c(93) = 0
930 b(c(16)) = 0: c(16) = 0
160 b(c(94)) = 0: c(94) = 0
940 b(c(17)) = 0: c(17) = 0
170 b(c(95)) = 0: c(95) = 0
950 b(c(18)) = 0: c(18) = 0
180 b(c(96)) = 0: c(96) = 0
960 Next j96
b(c(31)) = 0: c(31) = 0
310 b(c(97)) = 0: c(97) = 0
970 b(c(32)) = 0: c(32) = 0
320 b(c(98)) = 0: c(98) = 0
980 b(c(33)) = 0: c(33) = 0
330 b(c(99)) = 0: c(99) = 0
990 b(c(34)) = 0: c(34) = 0
340 b(c(100)) = 0: c(100) = 0
1000 b(c(35)) = 0: c(35) = 0
350 b(c(101)) = 0: c(101) = 0
1010 b(c(36)) = 0: c(36) = 0
360 b(c(102)) = 0: c(102) = 0
1020 b(c(25)) = 0: c(25) = 0
250 b(c(103)) = 0: c(103) = 0
1030 b(c(26)) = 0: c(26) = 0
260 b(c(104)) = 0: c(104) = 0
1040 b(c(27)) = 0: c(27) = 0
270 b(c(105)) = 0: c(105) = 0
1050 b(c(28)) = 0: c(28) = 0
280 b(c(106)) = 0: c(106) = 0
1060 b(c(29)) = 0: c(29) = 0
290 b(c(107)) = 0: c(107) = 0
1070 b(c(30)) = 0: c(30) = 0
300 b(c(108)) = 0: c(108) = 0
1080 Next j108
b(c(43)) = 0: c(43) = 0
430 b(c(109)) = 0: c(109) = 0
1090 b(c(44)) = 0: c(44) = 0
440 b(c(110)) = 0: c(110) = 0
1100 b(c(45)) = 0: c(45) = 0
450 b(c(111)) = 0: c(111) = 0
1110 b(c(46)) = 0: c(46) = 0
460 b(c(112)) = 0: c(112) = 0
1120 b(c(47)) = 0: c(47) = 0
470 b(c(113)) = 0: c(113) = 0
1130 b(c(48)) = 0: c(48) = 0
480 b(c(114)) = 0: c(114) = 0
1140 b(c(37)) = 0: c(37) = 0
370 b(c(115)) = 0: c(115) = 0
1150 b(c(38)) = 0: c(38) = 0
380 b(c(116)) = 0: c(116) = 0
1160 b(c(39)) = 0: c(39) = 0
390 b(c(117)) = 0: c(117) = 0
1170 b(c(40)) = 0: c(40) = 0
400 b(c(118)) = 0: c(118) = 0
1180 b(c(41)) = 0: c(41) = 0
410 b(c(119)) = 0: c(119) = 0
1190 b(c(42)) = 0: c(42) = 0
420 b(c(120)) = 0: c(120) = 0
1200 Next j120
b(c(55)) = 0: c(55) = 0
550 b(c(121)) = 0: c(121) = 0
1210 b(c(56)) = 0: c(56) = 0
560 b(c(122)) = 0: c(122) = 0
1220 b(c(57)) = 0: c(57) = 0
570 b(c(123)) = 0: c(123) = 0
1230 b(c(58)) = 0: c(58) = 0
580 b(c(124)) = 0: c(124) = 0
1240 b(c(59)) = 0: c(59) = 0
590 b(c(125)) = 0: c(125) = 0
1250 b(c(60)) = 0: c(60) = 0
600 b(c(126)) = 0: c(126) = 0
1260 b(c(49)) = 0: c(49) = 0
490 b(c(127)) = 0: c(127) = 0
1270 b(c(50)) = 0: c(50) = 0
500 b(c(128)) = 0: c(128) = 0
1280 b(c(51)) = 0: c(51) = 0
510 b(c(129)) = 0: c(129) = 0
1290 b(c(52)) = 0: c(52) = 0
520 b(c(130)) = 0: c(130) = 0
1300 b(c(53)) = 0: c(53) = 0
530 b(c(131)) = 0: c(131) = 0
1310 b(c(54)) = 0: c(54) = 0
540 b(c(132)) = 0: c(132) = 0
1320 Next j132
b(c(67)) = 0: c(67) = 0
670 b(c(133)) = 0: c(133) = 0
1330 b(c(68)) = 0: c(68) = 0
680 b(c(134)) = 0: c(134) = 0
1340 b(c(69)) = 0: c(69) = 0
690 b(c(135)) = 0: c(135) = 0
1350 b(c(70)) = 0: c(70) = 0
700 b(c(136)) = 0: c(136) = 0
1360 b(c(71)) = 0: c(71) = 0
710 b(c(137)) = 0: c(137) = 0
1370 b(c(72)) = 0: c(72) = 0
720 b(c(138)) = 0: c(138) = 0
1380 Next j138
b(c(61)) = 0: c(61) = 0
610 b(c(139)) = 0: c(139) = 0
1390 b(c(62)) = 0: c(62) = 0
620 b(c(140)) = 0: c(140) = 0
1400 Next j140
b(c(63)) = 0: c(63) = 0
630 b(c(141)) = 0: c(141) = 0
1410 Next j141
b(c(64)) = 0: c(64) = 0
640 b(c(142)) = 0: c(142) = 0
1420 Next j142
b(c(65)) = 0: c(65) = 0
650 b(c(143)) = 0: c(143) = 0
1430 Next j143
b(c(66)) = 0: c(66) = 0
660 b(c(144)) = 0: c(144) = 0
1440 Next j144
t2 = Timer
t10 = Str(t2 - t1) + " sec., " + Str(n9) + " Solutions for sum" + Str(s1)
y = MsgBox(t10, 0, "Routine MgcSqr12f")
End
' Define Natural Numbers
2030 Erase b1
For i1 = 1 To 144
a1(i1) = i1: b1(i1) = i1
Next i1
m1 = 1: m2 = 144
s1 = 870: p12 = s1 / 6
Return
' Print results (selected numbers)
2645 For i1 = 1 To 144
Cells(n9, i1).Value = a(i1)
Next i1
Cells(n9, 145).Value = n9
Cells(1, 146).Value = n9
Return
' Print results (squares)
2650 n2 = n2 + 1
If n2 = 3 Then
n2 = 1: k1 = k1 + 13: k2 = 1
Else
If n9 > 1 Then k2 = k2 + 13
End If
Cells(k1, k2 + 1).Select
Cells(k1, k2 + 1).Font.Color = -4165632
Cells(k1, k2 + 1).Value = s1
i3 = 0
For i1 = 1 To 12
For i2 = 1 To 12
i3 = i3 + 1
Cells(k1 + i1, k2 + i2).Value = a(i3)
Next i2
Next i1
Return
End Sub
