自己顶 K~c^*;F
1 REM "G-2-2" Lf
>YdD
2 CLEAR : LPRINT "** "; $ye^uu;Z
3 LPRINT "弹性地基梁计算程序(G-2-2)"; /b.$jnqL
4 LPRINT " **" Nd`HB=ShJ
5 OPEN "I",#1,"G-2-2.DAT" *^'wFbaBO
40 INPUT #1, R,N,LM,L,E0,E,U,B ]Ap`
42 INPUT #1, NH: DIM X(NH),H(NH) RCzV5g
43 FOR I1=1 TO NH: INPUT #1, X(I1),H(I1): NEXT I1 ~}BJ0P(VMc
45 DIM Z(R,4),AK(N+2),Y(N),D(N+2),B1(N+2) wL;]1&Qq
47 DIM F((N+2)*(N+3)/2),X2(N*(N+1)/2) wGPotPdE2
50 FOR I=1 TO R: FOR J=1 TO 4: INPUT #1, Z(I,J): NEXT J: NEXT I (nLKQV 1
55 INPUT #1, RL: IF RL=0 THEN 65 a`#S|'oatC
60 DIM PL(RL),AL(RL): FOR I=1 TO RL: INPUT #1, PL(I),AL(I): NEXT I 8QV t,
'I
65 INPUT #1, RR: IF RR=0 THEN 75 Hk,lX r
70 DIM PR(RR),AR(RR): FOR I=1 TO RR: INPUT #1, PR(I),AR(I): NEXT I
n
]w7Zj
75 INPUT #1, U0 AzZJG v]H
76 INPUT "需要打印原始数据吗 (Y/N)";G$ zK=dzoy
80 IF G$="N" OR G$="n" THEN 152 MRQZIi
85 LPRINT "原始数据:" tWm> j
90 LPRINT "地基的压缩模量 E0=";E0;"(KN/m?)" d7* CwY9"
95 LPRINT "梁的弹性模量 E=";E;"(KN/m?)" XdpF&B&K7Q
100 LPRINT "梁的泊松比 U=";U; "梁宽 B=";B;"(m)" uJgI<l'|e3
105 IF NH=1 THEN 113 hp*<x4%*a"
107 LPRINT : LPRINT TAB(2);"X(m)";TAB(11);"梁高H(m)" *ksb?|<Ot
108 X(0)=0 : FOR I=1 TO NH: LPRINT TAB(1);X(I-1);"-";X(I);TAB(14);H(I) ( :{"C6x
110 NEXT I: LPRINT : GOTO 115 'Q5&5UrBr
113 LPRINT "梁高 H=";H;"(m)" 1!!\+
c2*
115 LPRINT TAB(1);"梁长L=";L;"(m)" RcQ>eZHl
116 LPRINT TAB(1);"分段数N=";N;TAB(13);"荷载个数R=";R E#8_hT]5
120 LPRINT TAB(1);"问题类型(空间/平面/文克勒尔(1/2/3))LM=";LM sxo;/~.p
122 LPRINT "梁上荷载要素:" 8x9Rm
125 LPRINT TAB(2);"Z(I,1)";TAB(11);"Z(I,2)";TAB(20);"Z(I,3)";TAB(29);"Z(I,4)" GR<c=
127 FOR I=1 TO R: LPRINT TAB(1);Z(I,1);TAB(10);Z(I,2);TAB(19);Z(I,3);TAB(28);Z(I,4) AM'gnP>
129 NEXT I: LPRINT *8PN!^
130 IF RL=0 THEN 135 ,+s e
131 LPRINT "边荷载要素:" <;%0T
xK|U
132 LPRINT TAB(2);"左PL(I)";TAB(11);"AL(I)" tw=K&/@^O
133 FOR I=1 TO RL: LPRINT TAB(1);PL(I);TAB(10);AL(I):NEXT I: LPRINT m4@MxQm
135 IF RR=0 THEN 140 m%pBXXfGYj
136 LPRINT TAB(2);"右PR(I)";TAB(11);"AR(I)" Zd[OWF
138 FOR I=1 TO RR: LPRINT TAB(1);PR(I);TAB(10);AR(I):NEXT I: LPRINT B| tzF0;c
140 LPRINT "U0=";U0;"(地基泊松比,文氏问题为梁下可压缩层厚度) ibXe"X/_
145 LPRINT "- - - - - - - - -": H'EY)s Hi
152 IF LM<>2 THEN 165 l7z6i*R
154 INPUT "按平面应力/平面应变计算 (L/B)";A$ ),=@q+{E{
156 IF A$="L" OR A$="l" THEN LPRINT "按平面应力计算":GOTO 165 Nsn~@.UuSW
158 IF A$="B" OR A$="b" THEN LPRINT "按平面应变计算":GOTO 160 l4KbTKm7
159 GOTO 154 ACb/ITu
160 E=E/(1-U^2):U=U/(1-U):E0=E0/(1-U0^2):U0=U0/(1-U0) >B_n/v3P(M
165 C=L/N: FOR K=1 TO N r0\?WoF2C
170 AK(K)=(K-.5)*C:S=0 bj}=8k0
175 FOR I=1 TO R:T=Z(I,4): ON T GOTO 180,185 ZHCr2^w6
180 P1=Z(I,1):B2=Z(I,2): GOSUB 1740:S=S-P1*S1: GOTO 350 %(`#A.yaE
185 P1=Z(I,1)*(Z(I,3)-Z(I,2)):L1=(Z(I,3)+Z(I,2))/2 [SkKz>rC
190 I1=1: IF AK(K)>Z(I,2) THEN 205 rnO0-h-;
192 M1=L1-X(I1-1):M3=AK(K)-X(I1-1):JK=E*B*H(I1)^3/12 DR@1z9 a
195 IF AK(K)>X(I1) THEN 200 kh3<V'k]
197 M2=L1-AK(K):S=S-P1*(2*M1*M3+M2*M3)*M3/6/JK: GOTO 350 4=xq:Tf
200 M2=L1-X(I1):M4=AK(K)-X(I1):S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK $o.Kn9\
202 I1=I1+1: GOTO 192 /THnfy\
205 IF AK(K)>Z(I,3) THEN 295 pj!:[d
210 IF X(I1)>Z(I,2) THEN 225 n_P3\Y|
212 M1=L1-X(I1-1):M2=L1-X(I1):M3=AK(K)-X(I1-1):M4=AK(K)-X(I1):JK=E*B*H(I1)^3/12 7d;pvhnH
215 S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK `4V"s-T'
220 I1=I1+1: GOTO 210 \0ov[T N.>
225 IF X(I1-1)>=Z(I,2) THEN 265 h*3{6X#(/
230 M1=L1-X(I1-1):M2=L1-Z(I,2):M3=AK(K)-X(I1-1):M4=AK(K)-Z(I,2) vZu~LW@1
231 JK=E*B*H(I1)^3/12 W`rMtzL5
235 S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK *"cD.)]#2
240 IF X(I1)<=AK(K) THEN 250 o>F*Itr{
245 A1=Z(I,2):A2=AK(K): GOSUB 1900:S=S-Z(I,1)*S2/JK: GOTO 350 gRAC d&)
250 A1=Z(I,2):A2=X(I1): GOSUB 1900:S=S-Z(I,1)*S2/JK FW#P*}#
255 IF X(I1)=AK(K) THEN 350 cwe1^SJ6y
260 I1=I1+1 F\:(*1C
265 JK=E*B*H(I1)^3/12 3yA2WW
270 IF X(I1)<=AK(K) THEN 280 <rIz Z'D
275 A1=X(I1-1):A2=AK(K): GOSUB 1900:S=S-Z(I,1)*S2/JK: GOTO 350 &dR=?bz-A
280 A1=X(I1-1):A2=X(I1): GOSUB 1900:S=S-Z(I,1)*S2/JK -#s [F S
285 IF X(I1)=AK(K) THEN 350 DmqSQA
290 I1=I1+1: GOTO 265 hs+kr?Pg`
295 IF X(I1)>Z(I,2) THEN 315 MYLsHIPC
300 M1=L1-X(I1-1):M2=L1-X(I1):M3=AK(K)-X(I1-1):M4=AK(K)-X(I1):JK=E*B*H(I1)^3/12 ?okx<'"[
305 S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK a9U_ug58
310 I1=I1+1: GOTO 295 p$E8Bn%[
315 IF X(I1-1)>=Z(I,2) THEN 337 =ll=)"O
320 M1=L1-X(I1-1):M2=L1-Z(I,2):M3=AK(K)-X(I1-1):M4=AK(K)-Z(I,2) +&|S'7&{
321 JK=E*B*H(I1)^3/12 Kw5Lhc1V
323 S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK 9GLb"6+PK
325 IF X(I1)<=Z(I,3) THEN 330 +iO/m
327 A1=Z(I,2):A2=Z(I,3): GOSUB 1900:S=S-Z(I,1)*S2/JK: GOTO 350 0?0Jz
330 A1=Z(I,2):A2=X(I1): GOSUB 1900:S=S-Z(I,1)*S2/JK E?)656F[
332 IF X(I1)=Z(I,3) THEN 350 j&R+2%
335 I1=I1+1 hk>;pU(
337 JK=E*B*H(I1)^3/12: IF X(I1)>Z(I,3) THEN 348 )|bC^{kH!l
340 A1=X(I1-1):A2=X(I1): GOSUB 1900:S=S-Z(I,1)*S2/JK a
W%5~3
342 IF X(I1)=Z(I,3) THEN 350 51BlM%
345 I1=I1+1: GOTO 337 \qh*E#j
347 S=S-AK(K)^2*Z(I,1)/(2*E*JK) MXD4|r(
348 A1=X(I1-1):A2=Z(I,3): GOSUB 1900:S=S-Z(I,1)*S2/JK %^]?5a!
350 NEXT I kAA>FI6
351 IF LM=3 THEN 355 H%F>@(U
352 IF RL<>0 THEN 540 ciQZHH2
354 IF RR<>0 THEN 580 8,"yNq
355 S=S+SL+SR +}xaQc:0|
360 D(K)=-S:Y(K)=S h"+ `13
370 NEXT K v(]]_h
380 S=0:Q=0 .TJEUK
390 FOR I=1 TO R zDTv\3rZ4X
400 T=Z(I,4) >rsqH+oL
410 ON T GOTO 420,430,450 !g!5_|
420 S=S+Z(I,1):Q=Q+Z(I,1)*Z(I,2): GOTO 460 g\SrO {*
430 S=S+Z(I,1)*(Z(I,3)-Z(I,2)) gB'fFkd
440 Q=Q+Z(I,1)*(Z(I,3)-Z(I,2))*(Z(I,2)+(Z(I,3)-Z(I,2))/2): GOTO 460 ^P'{U26
450 Q=Q+Z(I,1) WjSc/3Qy
460 NEXT I I@#;nyAj"
470 D(N+1)=-S:D(N+2)=-Q Dnf*7)X
480 FOR K=1 TO N LE6.nmvS
485 FOR I=1 TO K IiBD?}
490 IF I>K THEN 530 .Dg*\ h
495 B2=AK(I): GOSUB 1740 ~]C%/gEh
500 F(K*(K-1)/2+I)=S1 8flOq"uK^
510 X2(K*(K-1)/2+I)=F(K*(K-1)/2+I) Ev"|FTI/
520 NEXT I 8;14Q7,S
530 NEXT K: GOTO 645 ?:{sH#ua
540 T=LM: ON T GOTO 545,560,570 RDqFL.-S
545 SL=0: FOR I=1 TO RL:X=AL(I)+AK(K): GOSUB 1800 4NN$( S-W
550 SL=SL+PL(I)*S0*(1-U0^2)/ 3.141592654 /E0/C: NEXT I mEyK1h1G@
555 GOTO 570 ) 0 W`
560 SL=0:ST=2* LOG (100*L/C)+2+ LOG (2) vQ}6y
562 FOR I=1 TO RL:X=AL(I)+AK(K): GOSUB 1850 yifY%!@Xu
565 S0=ST+S0:SL=SL+PL(I)*S0/ 3.141592654 /E0/B: NEXT I jnB~sbyA
570 GOTO 354 pCh2SQ(Q>
580 T=LM: ON T GOTO 585,600,610 $Fkaa<9;P
585 SR=0: FOR I=1 TO RR:X=AR(I)+L-AK(K): GOSUB 1800 B~
S6R
590 SR=SR+PR(I)*S0*(1-U0^2)/ 3.141592654 /E0/C: NEXT I 'B5^P
595 GOTO 610 ?S$i?\Qh
600 SR=0:ST=2* LOG (100*L/C)+2+ LOG (2) 6z0@I*
602 FOR I=1 TO RR:X=AR(I)+L-AK(K): GOSUB 1850 WYB{% yf
605 S0=ST+S0:SR=SR+PR(I)*S0/ 3.141592654 /E0/B: NEXT I L"jY+{oLIJ
610 GOTO 355 9^Fz iM
645 T=LM uDw.|B2ui
650 ON T GOTO 660,730,820 g-c ;}qz
660 Q=(1-U0^2)/( 3.141592654 *E0*C) (Rt7%{*
665 FOR K=1 TO N FTc.]laO
670 IF K-1<1 THEN 700 zaTb~#c_
675 FOR I=1 TO K-1 @yd4$Mv8%
680 X=AK(K)-AK(I): GOSUB 1800 gKLyL]kAGz
685 S=S0 NA3\
690 F(K*(K-1)/2+I)=F(K*(K-1)/2+I)+S*Q ;V^pL((5J
695 NEXT I mE`kjmX{ E
700 S= SQR (C^2/B^2+1):S=B* LOG (C/B+S)/C+ LOG (1+S) c~QS9)=E
705 F(K*(K+1)/2)=F(K*(K+1)/2)+2*C*( LOG (B/C)+S)*Q/B "tjLc6Xl^
710 NEXT K =@,Q Dm]L
720 GOTO 860 qyC=(v
730 Q=1/( 3.141593 *E0*B) 'r1LSht'
740 FOR K=1 TO N 6 DG@?O
745 IF K-1<1 THEN 800 -=s(l.?Hm5
750 FOR I=1 TO K-1 ~x#w<0e>
760 X=AK(K)-AK(I): GOSUB 1850 y[B>~m8$
770 S=S0 nOx4<Wk&
780 F(K*(K-1)/2+I)=F(K*(K-1)/2+I)+S*Q G,Z^g|6
790 NEXT I 7dsefNPb
800 NEXT K Lip4)Y [
810 GOTO 860 ,p(<+6QZ
820 Q=U0/(E0*B*C) akw,P$i
830 FOR K=1 TO N bVP"(H]
840 F(K*(K+1)/2)=F(K*(K+1)/2)+Q
rc&%m
850 NEXT K _@S`5;4x
860 FOR I=1 TO N kKC9{^%)
870 F((N+1)*N/2+I)=-1 kmzH'wktt
880 F((N+2)*(N+1)/2+I)=-AK(I) 3(C\.oRc
890 NEXT I gs!(;N\j|
900 F((N+1)*(N+2)/2)=0 w 4[{2
910 F((N+1)*(N+4)/2)=0 oh#\]c\f
920 F((N+2)*(N+3)/2)=0 8-<:i
940 N1=N+2 ]smkTo/
980 FOR I=1 TO N+2 Z{&cuo.@<]
990 B1(I)=D(I) s0Z
uWVip
1000 NEXT I X7k.zlH7T
1010 GOSUB 1520 iq(
)8nxi
1020 DIM QK(N),MK(N) y.KO :P?5{
1030 FOR K=1 TO N )95f*wte
1040 QK(K)=0:MK(K)=0 p<=$&*
1050 FOR I=1 TO R {( r6e
1060 T=Z(I,4) L(&&26Y
1070 ON T GOTO 1080,1100,1180 quY:pqG38q
1080 IF AK(K)>=Z(I,2) THEN 1090 MSf;ZB
1085 GOTO 1200 df7wN#kO+
1090 MK(K)=MK(K)-(AK(K)-Z(I,2))*Z(I,1):QK(K)=QK(K)-Z(I,1): GOTO 1200 y;/VB,4V
1100 IF AK(K)>Z(I,2) THEN 1120 Zd"^</ S
1110 GOTO 1200 Ua!aaq&
1120 IF AK(K)<Z(I,3) THEN 1160 boDt`2=
1130 QK(K)=QK(K)-Z(I,1)*(Z(I,3)-Z(I,2)) .\>v0Du
1140 MK(K)=MK(K)-Z(I,1)*(Z(I,3)-Z(I,2))*(AK(K)-(Z(I,2)+Z(I,3))/2) (5]}5W*
1150 GOTO 1200 <b,~:9*?
1160 QK(K)=QK(K)-Z(I,1)*(AK(K)-Z(I,2)) oudxm[/U
1170 MK(K)=MK(K)-Z(I,1)*(AK(K)-Z(I,2))^2/2: GOTO 1200 Lemui)
1180 IF AK(K)>Z(I,2) THEN MK(K)=MK(K)+Z(I,1) <&C]sb
1200 NEXT I iY21Ql%
1205 NEXT K *z@>!8?
1210 FOR K=1 TO N j?'GZ d"B
1220 S=0:Q=0 8OgLn?"P
1230 FOR I=1 TO N 13taFVdU
1240 IF I>K THEN Y(K)=X2(I*(I-1)/2+K)*B1(I)+Y(K): GOTO 1260 $Xq!L
1250 Y(K)=X2(K*(K-1)/2+I)*B1(I)+Y(K) 1GzAG;UUo6
1260 NEXT I Ub6jxib
1265 Y(K)=B1(N+1)+AK(K)*B1(N+2)-Y(K) iLSr*`
o
1267 IF K-1<1 THEN 1300 Gz.|]:1
1270 FOR I=1 TO K-1 ;*MLRXq
1280 S=S+B1(I):Q=Q+B1(I)*(AK(K)-AK(I)) UX7t`l2R
1290 NEXT I <)1qt
9
1300 QK(K)=QK(K)+S+.5*B1(K) dAuJXGo
1310 MK(K)=MK(K)+Q+C*B1(K)/8 p5G?N(l
1315 NEXT K ` V##Y
1320 FOR K=1 TO N .V,@k7U,V
1330 D(K)=B1(K)/C vN;mPd~g
1340 NEXT K KCs[/]
1350 LPRINT TAB(2);"X(m)";TAB(9);"反力D(kN/m?)";TAB(23);"弯距M(kN-m)";TAB(36);"剪力Q(kN)";TAB(49);"变位Y(cm)" ]\|VpIg
1360 FOR I=1 TO N -B +4+&{T
1370 LPRINT TAB(1); INT (100*(I-.5)*C+.5)/100;TAB(8); INT (100*D(I)+.5)/100; zCk^B/j sM
1372 LPRINT TAB(23); INT (100*MK(I)+.5)/100;TAB(36); INT (100*QK(I)+.5)/100; F w?[lS
1373 LPRINT TAB(49); INT (100*(Y(I)*100)+.5)/100 `nu''B
H
1375 NEXT I FJMrs[
1378 GOTO 1710 \-g)T}g,I
1380 LPRINT " THE END": |ZmUNiAa
1420 END VVlr*`
1500 G1=2*X/C+1 Zu7)gf
1520 DIM C1(N1-1) AIZ]jq
1525 FOR I=1 TO N1 >7 ="8
1527 IF I-1<1 THEN 1567 CB^U6ZS
1530 I1=I*(I-1)/2 @{25xTt
1535 FOR J=1 TO I-1 wRVUu)
1540 S=0:J1=J*(J-1)/2:J2=J*(J+1)/2 'dn]rV0(C
1543 IF J-1<1 THEN 1560 DMO Mh#[
1545 FOR K=1 TO J-1 vX.VfY
1550 S=S+C1(K)*F(J1+K) %KLpig
1555 NEXT K Z.L c>7o
1560 C1(J)=F(I1+J)-S:F(I1+J)=C1(J)/F(J2) 'tH_p
1565 NEXT J [@.!~E)P
1567 S=0:I2=I*(I+1)/2 ')cMiX\v
1568 IF I-1<1 THEN 1585 > ;*b|Ik
1570 FOR K=1 TO I-1 uLV#SQ=bZN
1575 S=S+C1(K)*F(I1+K) `x*Pof!Io
1580 NEXT K YuO.yh_
1585 F(I2)=F(I2)-S p>huRp^w
1590 NEXT I h'{ C[d
1595 FOR I=1 TO N1 x<ZJb
1600 I1=I*(I-1)/2:S=0 ,>%}B3O:Y=
1603 IF I-1<1 THEN 1625 F'21jy&
1605 FOR K=1 TO I-1 BI%$c~wS
1610 S=S+F(I1+K)*B1(K) H:V2[y8\
1615 NEXT K *_d7E
1620 B1(I)=B1(I)-S #AJM6* G9
1625 NEXT I vQ
6^xvk]
1630 FOR I=1 TO N1 xA$XT[D
1635 I2=I*(I+1)/2:B1(I)=B1(I)/F(I2) cPlZXf
1640 NEXT I glDu2a,Q
1645 FOR I=N1 TO 1 STEP -1 ,
K~}\CR
1647 IF I+1>N1 THEN 1668 ZQV6xoN;r
1650 S=0 A
PEE~
1655 FOR K=I+1 TO N1 \XZ/v*d0
1660 S=S+F(K*(K-1)/2+I)*B1(K) Tw-;7Ae
1663 NEXT K |IzPgC
1667 B1(I)=B1(I)-S:D(I)=B1(I)/C 8<QdMkI
1668 NEXT I ;@oN s-
1700 RETURN ^'PWI{ O
1710 GOTO 1380 X4~y7
1740 S1=0:I1=1: IF AK(K)>=B2 THEN 1770 b0Ps5G\ u
1745 M1=B2-X(I1-1):M3=AK(K)-X(I1-1):JK=E*B*H(I1)^3/12 e w$B)W
1750 IF X(I1)>AK(K) THEN 1765 ,s"^kFl
1755 M2=B2-X(I1):M4=AK(K)-X(I1):S1=S1+(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK S4_YT@VD%
1758 IF X(I1)=AK(K) THEN 1790 a.k.n<
1760 I1=I1+1: GOTO 1745 :74y!
1765 M2=B2-AK(K):S1=S1+(2*M1*M3+M2*M3)*(M1-M2)/6/JK Po^?QVJ7
1768 GOTO 1790 zBzZxK>$
1770 M1=B2-X(I1-1):M3=AK(K)-X(I1-1):JK=E*B*H(I1)^3/12 yb\_zE\
1772 IF X(I1)>B2 THEN 1785 z7fp#>uw
1775 M2=B2-X(I1):M4=AK(K)-X(I1):S1=S1+(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK Jdj2~pTq
1780 IF X(I1)=B2 THEN 1790 I&x=;
1782 I1=I1+1: GOTO 1770 8}O lL,fP
1785 M4=AK(K)-B2:S1=S1+(2*M1*M3+M1*M4)*(M3-M4)/6/JK -lY6|79bF
1790 RETURN <Zmg#
1800 G1=2*X/C+1:G2=2*X/C-1 lR6@
xJd:@
1805 G3=2*X/B:G4=G3-C/B:G3=G3+C/B n{argI8wF
1810 S0=2* LOG (B/C)- LOG (G1*G2)-2*X/C* LOG (G1/G2) <h0?tv]
1815 G2= SQR (G3^2+1):G1= SQR (G4^2+1) Sw ig;`
1820 S0=S0+B* LOG ((G3+G2)/(G4+G1))/C B|C2lu
1825 S0=S0+2*X* LOG ((1+G2)/(1+G1))/C c(xrP/yOwi
1830 S0=(S0+ LOG ((1+G2)*(1+G1)))/B*C 286jI7 T
1835 RETURN vN;N/mL
1850 G1=2*X/C+1 2K/4Rf0;
1855 S0=2*X/C-1:S0=-2*X* LOG (G1/S0)/C- LOG (G1*S0) M P Y[X[
1860 RETURN n FHUy9q
1900 S2=AK(K)*Z(I,3)^2*(A2-A1)/2-Z(I,3)*(2*AK(K)+Z(I,3))*(A2^2-A1^2)/4 ^ B fC
1905 S2=S2+(AK(K)+2*Z(I,3))*(A2^3-A1^3)/6-(A2^4-A1^4)/8 )q8p k2
1910 RETURN 3YOq2pW72G