自己顶 y*!8[wASHq
1 REM "G-2-2" (Rsf;VPO
2 CLEAR : LPRINT "** "; K0|:+s@ u
3 LPRINT "弹性地基梁计算程序(G-2-2)"; S5\KI+;PW
4 LPRINT " **" rP IAu[],g
5 OPEN "I",#1,"G-2-2.DAT" -`D<OSt7
40 INPUT #1, R,N,LM,L,E0,E,U,B K!(hj '0.
42 INPUT #1, NH: DIM X(NH),H(NH) q;.LK8M
43 FOR I1=1 TO NH: INPUT #1, X(I1),H(I1): NEXT I1 d%:J-UtG"
45 DIM Z(R,4),AK(N+2),Y(N),D(N+2),B1(N+2) vi]cl=S
47 DIM F((N+2)*(N+3)/2),X2(N*(N+1)/2) hcoZ5!LvT
50 FOR I=1 TO R: FOR J=1 TO 4: INPUT #1, Z(I,J): NEXT J: NEXT I Fg0!2MKq*
55 INPUT #1, RL: IF RL=0 THEN 65 Q _`EKz;N{
60 DIM PL(RL),AL(RL): FOR I=1 TO RL: INPUT #1, PL(I),AL(I): NEXT I hjz`0AS
65 INPUT #1, RR: IF RR=0 THEN 75 a H\A
70 DIM PR(RR),AR(RR): FOR I=1 TO RR: INPUT #1, PR(I),AR(I): NEXT I _k#!^AJ}x
75 INPUT #1, U0 1[!7xA0 j
76 INPUT "需要打印原始数据吗 (Y/N)";G$ &aPR" X
80 IF G$="N" OR G$="n" THEN 152 "gt1pf~y
85 LPRINT "原始数据:" Y+eDE:4
90 LPRINT "地基的压缩模量 E0=";E0;"(KN/m?)" Ro=dgQ0:t
95 LPRINT "梁的弹性模量 E=";E;"(KN/m?)" v oO7W "
100 LPRINT "梁的泊松比 U=";U; "梁宽 B=";B;"(m)" k\Y*tY#2
105 IF NH=1 THEN 113 y?Pw6;e.
107 LPRINT : LPRINT TAB(2);"X(m)";TAB(11);"梁高H(m)" [WI'o y
108 X(0)=0 : FOR I=1 TO NH: LPRINT TAB(1);X(I-1);"-";X(I);TAB(14);H(I) 6Zx 5^f(qd
110 NEXT I: LPRINT : GOTO 115 Q]<6voyy
113 LPRINT "梁高 H=";H;"(m)" .:O($9^Ho
115 LPRINT TAB(1);"梁长L=";L;"(m)" Bm 4$
116 LPRINT TAB(1);"分段数N=";N;TAB(13);"荷载个数R=";R e?p Q uF~
120 LPRINT TAB(1);"问题类型(空间/平面/文克勒尔(1/2/3))LM=";LM I~4!8W-Y
122 LPRINT "梁上荷载要素:" <{i1/"k?X
125 LPRINT TAB(2);"Z(I,1)";TAB(11);"Z(I,2)";TAB(20);"Z(I,3)";TAB(29);"Z(I,4)" zm3$)*p1
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) S
R s
129 NEXT I: LPRINT pTT00 `R
130 IF RL=0 THEN 135 sl:1P^b
131 LPRINT "边荷载要素:" T.W^L'L `
132 LPRINT TAB(2);"左PL(I)";TAB(11);"AL(I)" v6ei47-
133 FOR I=1 TO RL: LPRINT TAB(1);PL(I);TAB(10);AL(I):NEXT I: LPRINT n. vrq-
135 IF RR=0 THEN 140 rONz*ly|i
136 LPRINT TAB(2);"右PR(I)";TAB(11);"AR(I)" HX ,\a`
138 FOR I=1 TO RR: LPRINT TAB(1);PR(I);TAB(10);AR(I):NEXT I: LPRINT z:= E- +
140 LPRINT "U0=";U0;"(地基泊松比,文氏问题为梁下可压缩层厚度) c]"w0a-`^@
145 LPRINT "- - - - - - - - -": .)<l69ZD Z
152 IF LM<>2 THEN 165 _ jH./ @G
154 INPUT "按平面应力/平面应变计算 (L/B)";A$ Jp=
)L
156 IF A$="L" OR A$="l" THEN LPRINT "按平面应力计算":GOTO 165
Tj}%G
158 IF A$="B" OR A$="b" THEN LPRINT "按平面应变计算":GOTO 160 [=cYsW%WG
159 GOTO 154 & Z js
160 E=E/(1-U^2):U=U/(1-U):E0=E0/(1-U0^2):U0=U0/(1-U0) #Z,@yJ2wl
165 C=L/N: FOR K=1 TO N WWC& -Ni
170 AK(K)=(K-.5)*C:S=0 ([#'G+MC&
175 FOR I=1 TO R:T=Z(I,4): ON T GOTO 180,185 p) VMYu
180 P1=Z(I,1):B2=Z(I,2): GOSUB 1740:S=S-P1*S1: GOTO 350 O:ACp< @
185 P1=Z(I,1)*(Z(I,3)-Z(I,2)):L1=(Z(I,3)+Z(I,2))/2 =Y*@8=V
190 I1=1: IF AK(K)>Z(I,2) THEN 205 IC1nR
u2I
192 M1=L1-X(I1-1):M3=AK(K)-X(I1-1):JK=E*B*H(I1)^3/12 (}s& 8 4!
195 IF AK(K)>X(I1) THEN 200 y9k'jEZ"oh
197 M2=L1-AK(K):S=S-P1*(2*M1*M3+M2*M3)*M3/6/JK: GOTO 350 ^^< C9
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 eg H,7f(yP
202 I1=I1+1: GOTO 192 6n~) R
205 IF AK(K)>Z(I,3) THEN 295 V+r&Z<&
210 IF X(I1)>Z(I,2) THEN 225 dnV& U%fO
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 ].sD#~L_
215 S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK {R"mvB`
220 I1=I1+1: GOTO 210 Yj'"Wg
225 IF X(I1-1)>=Z(I,2) THEN 265 H?,Dv> .#*
230 M1=L1-X(I1-1):M2=L1-Z(I,2):M3=AK(K)-X(I1-1):M4=AK(K)-Z(I,2) :4)(Qa(
231 JK=E*B*H(I1)^3/12 *j=58d`n
235 S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK *9| p}q9n
240 IF X(I1)<=AK(K) THEN 250 ,5j3 (Lk
245 A1=Z(I,2):A2=AK(K): GOSUB 1900:S=S-Z(I,1)*S2/JK: GOTO 350 e,d}4 jy
250 A1=Z(I,2):A2=X(I1): GOSUB 1900:S=S-Z(I,1)*S2/JK XA;f.u
255 IF X(I1)=AK(K) THEN 350 2vhP'?;K
260 I1=I1+1 <[w5M?n8
265 JK=E*B*H(I1)^3/12 X_!km-{
270 IF X(I1)<=AK(K) THEN 280 Wf-XH |j[
275 A1=X(I1-1):A2=AK(K): GOSUB 1900:S=S-Z(I,1)*S2/JK: GOTO 350 / %gMzF
280 A1=X(I1-1):A2=X(I1): GOSUB 1900:S=S-Z(I,1)*S2/JK 1 ^30]2'_
285 IF X(I1)=AK(K) THEN 350 )2# q i/
290 I1=I1+1: GOTO 265 RJKi98xwJ
295 IF X(I1)>Z(I,2) THEN 315 f>e0 l'\
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 X&i;WI
305 S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK PHJHW#sv
310 I1=I1+1: GOTO 295 *!E~4z=
315 IF X(I1-1)>=Z(I,2) THEN 337 fs-LaV
0
320 M1=L1-X(I1-1):M2=L1-Z(I,2):M3=AK(K)-X(I1-1):M4=AK(K)-Z(I,2) IusZY B
321 JK=E*B*H(I1)^3/12
lnK
323 S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK k7j[t B#
325 IF X(I1)<=Z(I,3) THEN 330 f= .!/e7 0
327 A1=Z(I,2):A2=Z(I,3): GOSUB 1900:S=S-Z(I,1)*S2/JK: GOTO 350 ]{|lGtK %
330 A1=Z(I,2):A2=X(I1): GOSUB 1900:S=S-Z(I,1)*S2/JK [Jj@A(Cz
332 IF X(I1)=Z(I,3) THEN 350 h<bhH=6~
335 I1=I1+1 w'X N<RWA
337 JK=E*B*H(I1)^3/12: IF X(I1)>Z(I,3) THEN 348 @/: 7G.
340 A1=X(I1-1):A2=X(I1): GOSUB 1900:S=S-Z(I,1)*S2/JK EhW"s%Q
342 IF X(I1)=Z(I,3) THEN 350 q8n@ fi6
345 I1=I1+1: GOTO 337 W*Ow%$%2
347 S=S-AK(K)^2*Z(I,1)/(2*E*JK) F1BXu@~e(
348 A1=X(I1-1):A2=Z(I,3): GOSUB 1900:S=S-Z(I,1)*S2/JK $G{j[iLY
350 NEXT I VMe
351 IF LM=3 THEN 355 ?\(qA+iP0
352 IF RL<>0 THEN 540 ^Wn+G8n
354 IF RR<>0 THEN 580 o^Lq8u;i*
355 S=S+SL+SR 2-jXj9kp`
360 D(K)=-S:Y(K)=S kR:kn:
370 NEXT K !" &-k:|g
380 S=0:Q=0 6R guUDRQ
390 FOR I=1 TO R s= 5 k7
400 T=Z(I,4) (h= ]Ox
410 ON T GOTO 420,430,450 a)c;z@r
420 S=S+Z(I,1):Q=Q+Z(I,1)*Z(I,2): GOTO 460 !YI< A\P
430 S=S+Z(I,1)*(Z(I,3)-Z(I,2)) sI h5 cT
440 Q=Q+Z(I,1)*(Z(I,3)-Z(I,2))*(Z(I,2)+(Z(I,3)-Z(I,2))/2): GOTO 460 JYmYX-
450 Q=Q+Z(I,1) H=w):kL|
460 NEXT I lw99{y3<<
470 D(N+1)=-S:D(N+2)=-Q A{ M 7
480 FOR K=1 TO N -~4r 6ZcA
485 FOR I=1 TO K &[
],rT
490 IF I>K THEN 530 b)M- q{
495 B2=AK(I): GOSUB 1740 w w[||
=
500 F(K*(K-1)/2+I)=S1 <QC7HR
510 X2(K*(K-1)/2+I)=F(K*(K-1)/2+I) 7a^D[f0V
520 NEXT I U OGjil{.
530 NEXT K: GOTO 645 ,>za|y<n
540 T=LM: ON T GOTO 545,560,570 ^fH)E"qq5
545 SL=0: FOR I=1 TO RL:X=AL(I)+AK(K): GOSUB 1800 :I7mM y*
550 SL=SL+PL(I)*S0*(1-U0^2)/ 3.141592654 /E0/C: NEXT I {Qa O\{J=
555 GOTO 570 pLCS\AUTsv
560 SL=0:ST=2* LOG (100*L/C)+2+ LOG (2) i$`o,m#
562 FOR I=1 TO RL:X=AL(I)+AK(K): GOSUB 1850 jSH.e?
565 S0=ST+S0:SL=SL+PL(I)*S0/ 3.141592654 /E0/B: NEXT I #:P$a%V
570 GOTO 354 R4P&r=?
580 T=LM: ON T GOTO 585,600,610 B3D }'<
585 SR=0: FOR I=1 TO RR:X=AR(I)+L-AK(K): GOSUB 1800 BFhEDkk
590 SR=SR+PR(I)*S0*(1-U0^2)/ 3.141592654 /E0/C: NEXT I
kPiY |EH
595 GOTO 610 f F)M'C
600 SR=0:ST=2* LOG (100*L/C)+2+ LOG (2) s6Dkh}:d
602 FOR I=1 TO RR:X=AR(I)+L-AK(K): GOSUB 1850 GB<.kOGQ[
605 S0=ST+S0:SR=SR+PR(I)*S0/ 3.141592654 /E0/B: NEXT I /U0Hk>$~(
610 GOTO 355 z^j mf_
645 T=LM 023uAaI^3r
650 ON T GOTO 660,730,820 = :zPT;K
660 Q=(1-U0^2)/( 3.141592654 *E0*C) Bg[yn<)
]
665 FOR K=1 TO N -Eig#]Se3
670 IF K-1<1 THEN 700 t'0&n3
675 FOR I=1 TO K-1 \8m9^Z7IfK
680 X=AK(K)-AK(I): GOSUB 1800 7U1 M;@y
685 S=S0 "Nb2[R
690 F(K*(K-1)/2+I)=F(K*(K-1)/2+I)+S*Q 3$q#^UvD
695 NEXT I NZ&ZK@h}.
700 S= SQR (C^2/B^2+1):S=B* LOG (C/B+S)/C+ LOG (1+S) b 9"t%R9/Q
705 F(K*(K+1)/2)=F(K*(K+1)/2)+2*C*( LOG (B/C)+S)*Q/B );_ /0:
710 NEXT K >~]|o
720 GOTO 860 kB=B?V~ #
730 Q=1/( 3.141593 *E0*B) @$ Nti>
740 FOR K=1 TO N jkta]#O
745 IF K-1<1 THEN 800 F$jy~W_
750 FOR I=1 TO K-1 bvrXz-j
760 X=AK(K)-AK(I): GOSUB 1850 n4M
Xa()P1
770 S=S0 I )vR
780 F(K*(K-1)/2+I)=F(K*(K-1)/2+I)+S*Q o| vL:| 8Q
790 NEXT I 0 #Pa; (
800 NEXT K l4.ql1BX@y
810 GOTO 860 MRK=\qjD
820 Q=U0/(E0*B*C) (*=>YE'V{
830 FOR K=1 TO N _-TW-{7bh
840 F(K*(K+1)/2)=F(K*(K+1)/2)+Q >o v# \
850 NEXT K r QGInzYp
860 FOR I=1 TO N !lL
`L \
870 F((N+1)*N/2+I)=-1
T3<1{"&
880 F((N+2)*(N+1)/2+I)=-AK(I) 4lpcJ+:o
890 NEXT I H!hd0.
900 F((N+1)*(N+2)/2)=0 q~5zv4NX
910 F((N+1)*(N+4)/2)=0 9-q> W
920 F((N+2)*(N+3)/2)=0 P0(LdZH6u
940 N1=N+2 = 'w 2"4
980 FOR I=1 TO N+2 OFmHj] I7=
990 B1(I)=D(I) jO5R ~O`
1000 NEXT I s8 MQ:eAP
1010 GOSUB 1520 !G\gqkSL
1020 DIM QK(N),MK(N) V|B4lGS&
1030 FOR K=1 TO N o _l_Yi
1040 QK(K)=0:MK(K)=0 8 >LDo"<
1050 FOR I=1 TO R vu.?@k@
1060 T=Z(I,4) @ #hvQ6u
1070 ON T GOTO 1080,1100,1180 e f&8L
1080 IF AK(K)>=Z(I,2) THEN 1090 ! .q99DB
1085 GOTO 1200 Wa.xm_4s2
1090 MK(K)=MK(K)-(AK(K)-Z(I,2))*Z(I,1):QK(K)=QK(K)-Z(I,1): GOTO 1200 "4Q_F3?_`
1100 IF AK(K)>Z(I,2) THEN 1120 3.hFYA w
1110 GOTO 1200 oqysfLJ
1120 IF AK(K)<Z(I,3) THEN 1160 &4} =@'G@
1130 QK(K)=QK(K)-Z(I,1)*(Z(I,3)-Z(I,2)) hY%} x5ntU
1140 MK(K)=MK(K)-Z(I,1)*(Z(I,3)-Z(I,2))*(AK(K)-(Z(I,2)+Z(I,3))/2) 3 {t[>O;
1150 GOTO 1200 Vp~ cN
1160 QK(K)=QK(K)-Z(I,1)*(AK(K)-Z(I,2)) ,%N[FZ`|
1170 MK(K)=MK(K)-Z(I,1)*(AK(K)-Z(I,2))^2/2: GOTO 1200 C96*,.j~'
1180 IF AK(K)>Z(I,2) THEN MK(K)=MK(K)+Z(I,1) 4e
eh+T
1200 NEXT I 2gkN\w6zQ
1205 NEXT K %+,*$wk#*
1210 FOR K=1 TO N YEaT_zWG0
1220 S=0:Q=0 xyM|q 9Gf@
1230 FOR I=1 TO N !X$19"
1240 IF I>K THEN Y(K)=X2(I*(I-1)/2+K)*B1(I)+Y(K): GOTO 1260 %,z;W-#gnY
1250 Y(K)=X2(K*(K-1)/2+I)*B1(I)+Y(K) E@xrn+L>-
1260 NEXT I .I_<\h7
1265 Y(K)=B1(N+1)+AK(K)*B1(N+2)-Y(K) :aaX Y:<
1267 IF K-1<1 THEN 1300 Oso**WUOZ&
1270 FOR I=1 TO K-1
kSc~gJrne
1280 S=S+B1(I):Q=Q+B1(I)*(AK(K)-AK(I)) BY [7`@
1290 NEXT I ^VG].6
1300 QK(K)=QK(K)+S+.5*B1(K) z,E` +a;
1310 MK(K)=MK(K)+Q+C*B1(K)/8 EirZ}fDJzB
1315 NEXT K 73 .+0x
1320 FOR K=1 TO N 4lc|~Fj++
1330 D(K)=B1(K)/C ]1>R8
1340 NEXT K uKXD(lzX
1350 LPRINT TAB(2);"X(m)";TAB(9);"反力D(kN/m?)";TAB(23);"弯距M(kN-m)";TAB(36);"剪力Q(kN)";TAB(49);"变位Y(cm)" :: #[lw
1360 FOR I=1 TO N 7%? bl
1370 LPRINT TAB(1); INT (100*(I-.5)*C+.5)/100;TAB(8); INT (100*D(I)+.5)/100; OaY8 9 ko
1372 LPRINT TAB(23); INT (100*MK(I)+.5)/100;TAB(36); INT (100*QK(I)+.5)/100; }0]uA|lH*
1373 LPRINT TAB(49); INT (100*(Y(I)*100)+.5)/100 s !2Iui
@
1375 NEXT I p U !:
1378 GOTO 1710 k[ D,du')
1380 LPRINT " THE END": ru 9@|FgAE
1420 END .81Y/Gad_
1500 G1=2*X/C+1 ~&RrlF h
1520 DIM C1(N1-1) x*'H@!! G
1525 FOR I=1 TO N1 dL 'oKh,
1527 IF I-1<1 THEN 1567 : 7"Q
1530 I1=I*(I-1)/2 Y G8C<g6E7
1535 FOR J=1 TO I-1 )D/ ,QWk
1540 S=0:J1=J*(J-1)/2:J2=J*(J+1)/2 qWWt5rJ
1543 IF J-1<1 THEN 1560 u*I'c2m
1545 FOR K=1 TO J-1 & C) 1(
1550 S=S+C1(K)*F(J1+K) }i^|.VZZ
1555 NEXT K b-#{O=B
1560 C1(J)=F(I1+J)-S:F(I1+J)=C1(J)/F(J2) [ei~Xkzkj
1565 NEXT J I}|a7,8
1567 S=0:I2=I*(I+1)/2 3T0-RP*
1568 IF I-1<1 THEN 1585 o/V T"cT
1570 FOR K=1 TO I-1 =fnBE`Uc
1575 S=S+C1(K)*F(I1+K) /J_ ],KdU
1580 NEXT K ~@l4T_,k
1585 F(I2)=F(I2)-S bfoTGi
1590 NEXT I =j- ,yxBvJ
1595 FOR I=1 TO N1 #>)z}a]
1600 I1=I*(I-1)/2:S=0
]ilLed
1603 IF I-1<1 THEN 1625 j[A(@ w"
1605 FOR K=1 TO I-1 W]7<PL*u
1610 S=S+F(I1+K)*B1(K) f"j9C% '*
1615 NEXT K zz( |V
1620 B1(I)=B1(I)-S 9oGcbD4*
1625 NEXT I ' fb\t ,
1630 FOR I=1 TO N1 R?:Q=7K
1635 I2=I*(I+1)/2:B1(I)=B1(I)/F(I2) d^6-P
R_
1640 NEXT I i6n,N)%H
1645 FOR I=N1 TO 1 STEP -1 "p;tj7 4O9
1647 IF I+1>N1 THEN 1668 ,'`yh|}G\
1650 S=0 EZI#CLT[
1655 FOR K=I+1 TO N1 Y %bb-|\W
1660 S=S+F(K*(K-1)/2+I)*B1(K) Ot\[Ya ''
1663 NEXT K nM .g8d K
1667 B1(I)=B1(I)-S:D(I)=B1(I)/C 4D<C;>*/b
1668 NEXT I pQ0*)}l,
1700 RETURN l P= I0A-
1710 GOTO 1380 YU[#4f~
1740 S1=0:I1=1: IF AK(K)>=B2 THEN 1770 #'> )?]tn
1745 M1=B2-X(I1-1):M3=AK(K)-X(I1-1):JK=E*B*H(I1)^3/12 d;gs1]E50
1750 IF X(I1)>AK(K) THEN 1765 E8#aE\'t
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 t'At9<ib
1758 IF X(I1)=AK(K) THEN 1790 Ym\<@[3+!
1760 I1=I1+1: GOTO 1745 =0 EKrG
1765 M2=B2-AK(K):S1=S1+(2*M1*M3+M2*M3)*(M1-M2)/6/JK ([rn.b]
1768 GOTO 1790 > ]ux3F3\
1770 M1=B2-X(I1-1):M3=AK(K)-X(I1-1):JK=E*B*H(I1)^3/12 .4"BN<9
1772 IF X(I1)>B2 THEN 1785 nyR<pnuC'
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 RDHK'PGA
1780 IF X(I1)=B2 THEN 1790 8M,o)oH
1782 I1=I1+1: GOTO 1770 1irSI,j%z
1785 M4=AK(K)-B2:S1=S1+(2*M1*M3+M1*M4)*(M3-M4)/6/JK #\Y`?
1790 RETURN {X?1}5ry
1800 G1=2*X/C+1:G2=2*X/C-1 $},XRo&R
1805 G3=2*X/B:G4=G3-C/B:G3=G3+C/B W
= Bw*o-
1810 S0=2* LOG (B/C)- LOG (G1*G2)-2*X/C* LOG (G1/G2) l#^?sb G
1815 G2= SQR (G3^2+1):G1= SQR (G4^2+1) {p/Yz#
1820 S0=S0+B* LOG ((G3+G2)/(G4+G1))/C j[dZ*Jr_
1825 S0=S0+2*X* LOG ((1+G2)/(1+G1))/C zo "L9&Hzo
1830 S0=(S0+ LOG ((1+G2)*(1+G1)))/B*C U n)Xe
1835 RETURN ^ /G ;
1850 G1=2*X/C+1 3$YbEl@#
1855 S0=2*X/C-1:S0=-2*X* LOG (G1/S0)/C- LOG (G1*S0) {vW0O &[
1860 RETURN qQDe'f~
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 U7n#TPet
1905 S2=S2+(AK(K)+2*Z(I,3))*(A2^3-A1^3)/6-(A2^4-A1^4)/8 e
}C,)
1910 RETURN [DeD U:
