自己顶 ~.3v\Q
1 REM "G-2-2"
!D['}%
2 CLEAR : LPRINT "** "; E`_T_O=P
3 LPRINT "弹性地基梁计算程序(G-2-2)"; s.7=!JQ#]p
4 LPRINT " **" _ab8z]H
5 OPEN "I",#1,"G-2-2.DAT" %`k [xz
40 INPUT #1, R,N,LM,L,E0,E,U,B iw M xTty
42 INPUT #1, NH: DIM X(NH),H(NH) AR( gI]1
43 FOR I1=1 TO NH: INPUT #1, X(I1),H(I1): NEXT I1 +0U=UV)U
45 DIM Z(R,4),AK(N+2),Y(N),D(N+2),B1(N+2) bu6Sp3g
47 DIM F((N+2)*(N+3)/2),X2(N*(N+1)/2) nxhlTf>3
50 FOR I=1 TO R: FOR J=1 TO 4: INPUT #1, Z(I,J): NEXT J: NEXT I A{;"e^a-^l
55 INPUT #1, RL: IF RL=0 THEN 65 :y7K3:d3
60 DIM PL(RL),AL(RL): FOR I=1 TO RL: INPUT #1, PL(I),AL(I): NEXT I M`al~9
65 INPUT #1, RR: IF RR=0 THEN 75 P9
HKev?y
70 DIM PR(RR),AR(RR): FOR I=1 TO RR: INPUT #1, PR(I),AR(I): NEXT I !y XGA g,
75 INPUT #1, U0 &1<[@:;
76 INPUT "需要打印原始数据吗 (Y/N)";G$ 8GW+:
80 IF G$="N" OR G$="n" THEN 152 s? ;8h &]=
85 LPRINT "原始数据:" (rhlK}
C
90 LPRINT "地基的压缩模量 E0=";E0;"(KN/m?)" 5FJLDT2Lg
95 LPRINT "梁的弹性模量 E=";E;"(KN/m?)" + F o^NT
100 LPRINT "梁的泊松比 U=";U; "梁宽 B=";B;"(m)" " |I.j)
105 IF NH=1 THEN 113 gk| %
4.
107 LPRINT : LPRINT TAB(2);"X(m)";TAB(11);"梁高H(m)" $=d iG
108 X(0)=0 : FOR I=1 TO NH: LPRINT TAB(1);X(I-1);"-";X(I);TAB(14);H(I) o3+s.7 "
110 NEXT I: LPRINT : GOTO 115 hO[_ _j8
113 LPRINT "梁高 H=";H;"(m)" rP]| `*B
115 LPRINT TAB(1);"梁长L=";L;"(m)" J|{50?S{^
116 LPRINT TAB(1);"分段数N=";N;TAB(13);"荷载个数R=";R
^cw9Yjh6
120 LPRINT TAB(1);"问题类型(空间/平面/文克勒尔(1/2/3))LM=";LM ?WXftzdf6u
122 LPRINT "梁上荷载要素:" v|~=rvXFC
125 LPRINT TAB(2);"Z(I,1)";TAB(11);"Z(I,2)";TAB(20);"Z(I,3)";TAB(29);"Z(I,4)" S|| W
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) AJ6l#j-
129 NEXT I: LPRINT eEBNO*2
130 IF RL=0 THEN 135 Kw" e4 a
131 LPRINT "边荷载要素:" O F`J{`{r
132 LPRINT TAB(2);"左PL(I)";TAB(11);"AL(I)" <vWP_yy
133 FOR I=1 TO RL: LPRINT TAB(1);PL(I);TAB(10);AL(I):NEXT I: LPRINT nZnqXclzxn
135 IF RR=0 THEN 140 v3cMPN
136 LPRINT TAB(2);"右PR(I)";TAB(11);"AR(I)" TO89;O
138 FOR I=1 TO RR: LPRINT TAB(1);PR(I);TAB(10);AR(I):NEXT I: LPRINT ^z!=,M<+{
140 LPRINT "U0=";U0;"(地基泊松比,文氏问题为梁下可压缩层厚度) O EaL2T
145 LPRINT "- - - - - - - - -": BA1H)%
152 IF LM<>2 THEN 165 6oLOA}q
154 INPUT "按平面应力/平面应变计算 (L/B)";A$ n*$g1 HG6
156 IF A$="L" OR A$="l" THEN LPRINT "按平面应力计算":GOTO 165 eb`3'&zV&)
158 IF A$="B" OR A$="b" THEN LPRINT "按平面应变计算":GOTO 160 nuWQ3w
p[e
159 GOTO 154 )B
$Q
160 E=E/(1-U^2):U=U/(1-U):E0=E0/(1-U0^2):U0=U0/(1-U0) I1m[M?
165 C=L/N: FOR K=1 TO N QWa@?BO2p
170 AK(K)=(K-.5)*C:S=0 @P~%4:!Hr
175 FOR I=1 TO R:T=Z(I,4): ON T GOTO 180,185 .F$}a%
180 P1=Z(I,1):B2=Z(I,2): GOSUB 1740:S=S-P1*S1: GOTO 350 Pl/Xh03E
185 P1=Z(I,1)*(Z(I,3)-Z(I,2)):L1=(Z(I,3)+Z(I,2))/2 U9T}iI
190 I1=1: IF AK(K)>Z(I,2) THEN 205 /7"V~c6
192 M1=L1-X(I1-1):M3=AK(K)-X(I1-1):JK=E*B*H(I1)^3/12 $\*Z
195 IF AK(K)>X(I1) THEN 200 VsSAb%
197 M2=L1-AK(K):S=S-P1*(2*M1*M3+M2*M3)*M3/6/JK: GOTO 350 glCpA$;VPu
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 h[qZM
202 I1=I1+1: GOTO 192 az![u)
205 IF AK(K)>Z(I,3) THEN 295 ?7wcv$K5
210 IF X(I1)>Z(I,2) THEN 225 4G I3|{
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 k^ |z.$+
215 S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK F% a&|X
220 I1=I1+1: GOTO 210 }*rS g .
225 IF X(I1-1)>=Z(I,2) THEN 265 >vo 6X]p~
230 M1=L1-X(I1-1):M2=L1-Z(I,2):M3=AK(K)-X(I1-1):M4=AK(K)-Z(I,2) !;8Y?c-D
231 JK=E*B*H(I1)^3/12 -){6ynqv
235 S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK 'cc8 xC
240 IF X(I1)<=AK(K) THEN 250 {.k IC@^O
245 A1=Z(I,2):A2=AK(K): GOSUB 1900:S=S-Z(I,1)*S2/JK: GOTO 350 $"NH{%95}
250 A1=Z(I,2):A2=X(I1): GOSUB 1900:S=S-Z(I,1)*S2/JK }Fu1Y@M%
255 IF X(I1)=AK(K) THEN 350 er24}G8
260 I1=I1+1 X<IW5*
265 JK=E*B*H(I1)^3/12 ?"
4X&6xl
270 IF X(I1)<=AK(K) THEN 280 oS$7k3s
fj
275 A1=X(I1-1):A2=AK(K): GOSUB 1900:S=S-Z(I,1)*S2/JK: GOTO 350 8y6dT
280 A1=X(I1-1):A2=X(I1): GOSUB 1900:S=S-Z(I,1)*S2/JK 40MKf/9
285 IF X(I1)=AK(K) THEN 350 @ "NP`#
290 I1=I1+1: GOTO 265 PcA^ jBgGl
295 IF X(I1)>Z(I,2) THEN 315 xltN-<n7
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 EpG9t9S9
305 S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK %XN;S29d5W
310 I1=I1+1: GOTO 295 bL* b>R[x
315 IF X(I1-1)>=Z(I,2) THEN 337 -h7ssf'u[
320 M1=L1-X(I1-1):M2=L1-Z(I,2):M3=AK(K)-X(I1-1):M4=AK(K)-Z(I,2) Gr\jjf`
321 JK=E*B*H(I1)^3/12 ;b 65s9n^b
323 S=S-P1*(2*M1*M3+2*M2*M4+M1*M4+M2*M3)*(M1-M2)/6/JK [;IE Z/ZX
325 IF X(I1)<=Z(I,3) THEN 330 *w0|`[P+h
327 A1=Z(I,2):A2=Z(I,3): GOSUB 1900:S=S-Z(I,1)*S2/JK: GOTO 350 *cM=>3ws/
330 A1=Z(I,2):A2=X(I1): GOSUB 1900:S=S-Z(I,1)*S2/JK *(5;5r
332 IF X(I1)=Z(I,3) THEN 350 u QH]
335 I1=I1+1 ;UgwV/d
337 JK=E*B*H(I1)^3/12: IF X(I1)>Z(I,3) THEN 348 {+#{Cha
340 A1=X(I1-1):A2=X(I1): GOSUB 1900:S=S-Z(I,1)*S2/JK @k;65'"Q
342 IF X(I1)=Z(I,3) THEN 350 i|z= WnF$&
345 I1=I1+1: GOTO 337 Dk s n
347 S=S-AK(K)^2*Z(I,1)/(2*E*JK) @yb'h`f]
348 A1=X(I1-1):A2=Z(I,3): GOSUB 1900:S=S-Z(I,1)*S2/JK @&m]:GR
350 NEXT I M 2ex
3m
351 IF LM=3 THEN 355 jj2=|)w$3
352 IF RL<>0 THEN 540 G{6@] 72
354 IF RR<>0 THEN 580 kOo Vqu
355 S=S+SL+SR )jl@ hnA
360 D(K)=-S:Y(K)=S T8\@CV!
370 NEXT K TYD( 6 N
380 S=0:Q=0 mK $E&,OkA
390 FOR I=1 TO R !m:WoQ /
400 T=Z(I,4) j[i*;0) |
410 ON T GOTO 420,430,450 C zw ]5
420 S=S+Z(I,1):Q=Q+Z(I,1)*Z(I,2): GOTO 460 p5E
okh
430 S=S+Z(I,1)*(Z(I,3)-Z(I,2)) :'%|LBc0
440 Q=Q+Z(I,1)*(Z(I,3)-Z(I,2))*(Z(I,2)+(Z(I,3)-Z(I,2))/2): GOTO 460 .S|T{DM Q[
450 Q=Q+Z(I,1) |MKR&%Na
460 NEXT I j;uUM6
470 D(N+1)=-S:D(N+2)=-Q l HM}
E$5
480 FOR K=1 TO N >
"rM\ Q
485 FOR I=1 TO K 0~ nCT&V
490 IF I>K THEN 530 %[KnpJ{\
495 B2=AK(I): GOSUB 1740 Z< >gx m<
500 F(K*(K-1)/2+I)=S1 7.VP7;jys
510 X2(K*(K-1)/2+I)=F(K*(K-1)/2+I) 7r?,wM
520 NEXT I ]tu
OWR
530 NEXT K: GOTO 645 `:7r5}(^
540 T=LM: ON T GOTO 545,560,570 ssQ1u.x9
545 SL=0: FOR I=1 TO RL:X=AL(I)+AK(K): GOSUB 1800 W=A0+t% XC
550 SL=SL+PL(I)*S0*(1-U0^2)/ 3.141592654 /E0/C: NEXT I 3<<wH K;)
555 GOTO 570 1lIs
jBo g
560 SL=0:ST=2* LOG (100*L/C)+2+ LOG (2) 9^g8VlQdT
562 FOR I=1 TO RL:X=AL(I)+AK(K): GOSUB 1850 IY6Ll6OK
565 S0=ST+S0:SL=SL+PL(I)*S0/ 3.141592654 /E0/B: NEXT I sx azl]
570 GOTO 354 O,),0zcYF
580 T=LM: ON T GOTO 585,600,610 !VIxEu^ke
585 SR=0: FOR I=1 TO RR:X=AR(I)+L-AK(K): GOSUB 1800
MOB4t |
590 SR=SR+PR(I)*S0*(1-U0^2)/ 3.141592654 /E0/C: NEXT I <% 3fJt-Ie
595 GOTO 610 ]\K?%z
600 SR=0:ST=2* LOG (100*L/C)+2+ LOG (2) N[O .p]8
602 FOR I=1 TO RR:X=AR(I)+L-AK(K): GOSUB 1850 } 'xGip@W
605 S0=ST+S0:SR=SR+PR(I)*S0/ 3.141592654 /E0/B: NEXT I pD[&,gV$
610 GOTO 355 $/
"+t.ir3
645 T=LM ~SBW`=aP}
650 ON T GOTO 660,730,820 (BtU\ f#d
660 Q=(1-U0^2)/( 3.141592654 *E0*C) 9; XbyA]
665 FOR K=1 TO N eCKm4l'BZ
670 IF K-1<1 THEN 700 MV zj7~+
675 FOR I=1 TO K-1 E h;Ia6}
680 X=AK(K)-AK(I): GOSUB 1800 P =EZ6<c3&
685 S=S0 $:5h5Y#z
690 F(K*(K-1)/2+I)=F(K*(K-1)/2+I)+S*Q Gi-pi=#&cs
695 NEXT I r 6\g #}
700 S= SQR (C^2/B^2+1):S=B* LOG (C/B+S)/C+ LOG (1+S) D3aX\ NGP
705 F(K*(K+1)/2)=F(K*(K+1)/2)+2*C*( LOG (B/C)+S)*Q/B "Cxj_V@\
710 NEXT K KO8vUR*2R
720 GOTO 860 16eP7s
730 Q=1/( 3.141593 *E0*B) 2m*ugBO;
740 FOR K=1 TO N I`XOv SO
745 IF K-1<1 THEN 800 u q 9mq"
750 FOR I=1 TO K-1 G3|23G.~)(
760 X=AK(K)-AK(I): GOSUB 1850 !QAndg{;D
770 S=S0 En7+fQ
780 F(K*(K-1)/2+I)=F(K*(K-1)/2+I)+S*Q cd,'37 pZ
790 NEXT I 0^Ldw)C"
800 NEXT K cHr]{@7Cs
810 GOTO 860 **__&X p1
820 Q=U0/(E0*B*C) YIW9z{rrs
830 FOR K=1 TO N *0,*F ~n
840 F(K*(K+1)/2)=F(K*(K+1)/2)+Q X sJ`x
850 NEXT K "k + :!D
860 FOR I=1 TO N d(t)8k$
870 F((N+1)*N/2+I)=-1 :T$}@& -
880 F((N+2)*(N+1)/2+I)=-AK(I) }N -UlL(
890 NEXT I \mu';[gLd
900 F((N+1)*(N+2)/2)=0 XelFGT E
910 F((N+1)*(N+4)/2)=0 vM5I2C3_>!
920 F((N+2)*(N+3)/2)=0 s[{8:Px
940 N1=N+2 SE(<(w
980 FOR I=1 TO N+2 6 [XaIco=C
990 B1(I)=D(I) * Ib DA
1000 NEXT I {BM:c$3@j
1010 GOSUB 1520 QE"$Lc)
1020 DIM QK(N),MK(N) VB |k
1030 FOR K=1 TO N :| k!hG
1040 QK(K)=0:MK(K)=0 Mz$qe
1050 FOR I=1 TO R +7OE,RoQ
1060 T=Z(I,4) [wSoZBl
1070 ON T GOTO 1080,1100,1180 W:n \,P
1080 IF AK(K)>=Z(I,2) THEN 1090 U7fpaxc-
1085 GOTO 1200 ,B;mG ]_
1090 MK(K)=MK(K)-(AK(K)-Z(I,2))*Z(I,1):QK(K)=QK(K)-Z(I,1): GOTO 1200 hb~d4 J=S
1100 IF AK(K)>Z(I,2) THEN 1120 n%;qIKnIq\
1110 GOTO 1200 =CFg~8 W
1120 IF AK(K)<Z(I,3) THEN 1160 "?k'S{;
1130 QK(K)=QK(K)-Z(I,1)*(Z(I,3)-Z(I,2)) }`
`oojz
1140 MK(K)=MK(K)-Z(I,1)*(Z(I,3)-Z(I,2))*(AK(K)-(Z(I,2)+Z(I,3))/2) +,"[ 0RH
1150 GOTO 1200 PT,*KYF_O"
1160 QK(K)=QK(K)-Z(I,1)*(AK(K)-Z(I,2)) bny@AP(CY+
1170 MK(K)=MK(K)-Z(I,1)*(AK(K)-Z(I,2))^2/2: GOTO 1200 ,e$RvFB
1180 IF AK(K)>Z(I,2) THEN MK(K)=MK(K)+Z(I,1) rkS'OC
1200 NEXT I < hy!B4
1205 NEXT K +Q_xY>ej
1210 FOR K=1 TO N NM9Vi Ym>P
1220 S=0:Q=0 7ojh=i mY
1230 FOR I=1 TO N
Rq| 5%;1
1240 IF I>K THEN Y(K)=X2(I*(I-1)/2+K)*B1(I)+Y(K): GOTO 1260 =3hJti9[
1250 Y(K)=X2(K*(K-1)/2+I)*B1(I)+Y(K) RgFpc*.T
1260 NEXT I }mp`!7?>O
1265 Y(K)=B1(N+1)+AK(K)*B1(N+2)-Y(K) $/Wec ,`&
1267 IF K-1<1 THEN 1300 l]DRJ
1270 FOR I=1 TO K-1
PC@H Nto{
1280 S=S+B1(I):Q=Q+B1(I)*(AK(K)-AK(I)) oIOeX1$V
1290 NEXT I !V #*(_+n
1300 QK(K)=QK(K)+S+.5*B1(K) B> i^ w1
1310 MK(K)=MK(K)+Q+C*B1(K)/8 ?xKiN5q"6
1315 NEXT K #q34>}O< O
1320 FOR K=1 TO N u* iqwm.
1330 D(K)=B1(K)/C 6 T ~+vT
1340 NEXT K b *| ?7
1350 LPRINT TAB(2);"X(m)";TAB(9);"反力D(kN/m?)";TAB(23);"弯距M(kN-m)";TAB(36);"剪力Q(kN)";TAB(49);"变位Y(cm)" ~7}no}7
1360 FOR I=1 TO N |1r y*~
1370 LPRINT TAB(1); INT (100*(I-.5)*C+.5)/100;TAB(8); INT (100*D(I)+.5)/100; sR PQr ?
1372 LPRINT TAB(23); INT (100*MK(I)+.5)/100;TAB(36); INT (100*QK(I)+.5)/100; (*eX'^Q)d
1373 LPRINT TAB(49); INT (100*(Y(I)*100)+.5)/100 _d~GY,WTdO
1375 NEXT I UE_>@_T
1378 GOTO 1710 |:(B I5&S
1380 LPRINT " THE END": BSy4
d>
1420 END ;w% g*S
1500 G1=2*X/C+1 n @@tO#!\
1520 DIM C1(N1-1) q{*[uJ}Xc"
1525 FOR I=1 TO N1 `,pBOh|'
1527 IF I-1<1 THEN 1567 <F_w4!
1530 I1=I*(I-1)/2 fU.hb%m)Q\
1535 FOR J=1 TO I-1 }T902RL0
1540 S=0:J1=J*(J-1)/2:J2=J*(J+1)/2 .6n| hYe
1543 IF J-1<1 THEN 1560 vQXF$/S
1545 FOR K=1 TO J-1 w0js_P-uv
1550 S=S+C1(K)*F(J1+K) myXGMN$ i
1555 NEXT K sdXchVC
1560 C1(J)=F(I1+J)-S:F(I1+J)=C1(J)/F(J2) *URY8 a`bO
1565 NEXT J >@4Ds"Ye"O
1567 S=0:I2=I*(I+1)/2 eWYet2!Q
1568 IF I-1<1 THEN 1585 05 6yhB
1570 FOR K=1 TO I-1 W{ozZuo
1575 S=S+C1(K)*F(I1+K) n$j B "1
1580 NEXT K AS0(NlV
1585 F(I2)=F(I2)-S hHw1<! M
1590 NEXT I _kOuD}_|
1595 FOR I=1 TO N1 8_>: 0(y
1600 I1=I*(I-1)/2:S=0 ZPM7R3%V)z
1603 IF I-1<1 THEN 1625 nC/T$
#G
1605 FOR K=1 TO I-1 .Z QXY%g
1610 S=S+F(I1+K)*B1(K) $+ e (k~
1615 NEXT K FhH*lO&
1620 B1(I)=B1(I)-S {3vm ]
1625 NEXT I cQh{z8Bf?<
1630 FOR I=1 TO N1 Rbm+V{EF&
1635 I2=I*(I+1)/2:B1(I)=B1(I)/F(I2) 3H"F~_H
1640 NEXT I RToX[R;1E
1645 FOR I=N1 TO 1 STEP -1 -{eiV0<^
1647 IF I+1>N1 THEN 1668 ,trh)ZZYW|
1650 S=0
KYcc jX
1655 FOR K=I+1 TO N1 _mdJIa0D6k
1660 S=S+F(K*(K-1)/2+I)*B1(K) b2F1^]p
1663 NEXT K jku Nafp}
1667 B1(I)=B1(I)-S:D(I)=B1(I)/C %E, -dw
1668 NEXT I ) tV]h#4
1700 RETURN m|cRj{xZF
1710 GOTO 1380 $a\X( okx
1740 S1=0:I1=1: IF AK(K)>=B2 THEN 1770 jvd3_L-@E<
1745 M1=B2-X(I1-1):M3=AK(K)-X(I1-1):JK=E*B*H(I1)^3/12 *sfz +8Y
1750 IF X(I1)>AK(K) THEN 1765 0~<t :q!
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 !5m~qet.
1758 IF X(I1)=AK(K) THEN 1790 Vas Q/
1760 I1=I1+1: GOTO 1745 h*P0; V`UX
1765 M2=B2-AK(K):S1=S1+(2*M1*M3+M2*M3)*(M1-M2)/6/JK 3 5-FD{
1768 GOTO 1790 +f]I7e:qp
1770 M1=B2-X(I1-1):M3=AK(K)-X(I1-1):JK=E*B*H(I1)^3/12 *Z"Kvj;> u
1772 IF X(I1)>B2 THEN 1785 Jz;`L3m
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 /Jk.b/t.*S
1780 IF X(I1)=B2 THEN 1790 z SsogAx
1782 I1=I1+1: GOTO 1770 QKz2ONV=)
1785 M4=AK(K)-B2:S1=S1+(2*M1*M3+M1*M4)*(M3-M4)/6/JK (o6 u ^#6
1790 RETURN &U.y):
1800 G1=2*X/C+1:G2=2*X/C-1 Jl"DMUy[kW
1805 G3=2*X/B:G4=G3-C/B:G3=G3+C/B H-5f!>)
1810 S0=2* LOG (B/C)- LOG (G1*G2)-2*X/C* LOG (G1/G2) t@cBuV`9c
1815 G2= SQR (G3^2+1):G1= SQR (G4^2+1) ,h3,& ,
1820 S0=S0+B* LOG ((G3+G2)/(G4+G1))/C :i ?c
1825 S0=S0+2*X* LOG ((1+G2)/(1+G1))/C ;XYfw)
1830 S0=(S0+ LOG ((1+G2)*(1+G1)))/B*C Qw% 0<~<
1835 RETURN 3kJSz-_M
1850 G1=2*X/C+1 QIN# \
1855 S0=2*X/C-1:S0=-2*X* LOG (G1/S0)/C- LOG (G1*S0) \<%FZT_4~
1860 RETURN G r d9yLF
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
&@7|_60
1905 S2=S2+(AK(K)+2*Z(I,3))*(A2^3-A1^3)/6-(A2^4-A1^4)/8 *=E4|> Ul,
1910 RETURN K1<l/
s
