第一节使用说明1.数据的填写与穿孔(1)数据表(见附表1)：原始数据标识符含义：NO——题号，可用任意的不大于8位的整数表示；A——断面型号，按程序安排的代号填写；L d——底板半宽，m；H——侧墙高度，m；R1——底拱半径，m；R24——侧拱或底圆角半径，m；R4——顶拱半径，m；F i1——底拱半中心角，°；F i2——侧拱中心角，°；F i——顶拱半中心角，°；D1——底板或底拱厚度，m；D2——侧拱厚度，m；D3——侧墙厚度，m；D4B——顶拱拱脚厚度，m；D4——顶拱拱顶厚度，m；D5——顶板厚度，m；H0——水头，由底板中线算起，内压为正，外压为负，外压应为折减后的有效水头，m；Q n——顶部垂直山岩压力强度，向下为正，t/m；E1——水平山岩压力顶部强度，指向衬砌为正，t/m；E2——水平山岩压力底部强度，指向衬砌为正，t/m；Q d——底部山岩压力强度，向上为正，t/m；K1、K2、K3、K4、K5——底板、侧拱、侧墙、顶拱、顶板等部位，围岩弹性抗力系数，t/m3；M1——衬砌材料容重，t/m3；E——衬砌材料弹性模量，t/m2；B——无底板断面侧墙墙底约束情况(B=1铰支端；B=2固定端；B=3弹性固端)。

(2)数据填写按衬砌结构图形，在附表1中选用。

原始数据的个数与排列顺序，不可随意增减或倒置。

当无某项荷载时，可在相应的位置填零；不计衬砌自重时，材料容重填零；局部部位围岩破碎，不计弹性抗力时，该部位弹性抗力系数可填零。

附表1数据表TQ—16机用(3)数据穿孔使用ALGOL程序(TQ—16)机，数据用‘N’格式分两段穿孔，也可用机器允许的简化格式穿孔。

使用FORTRAN程序(FELIX机)，数据用卡片穿孔，前两个数据穿一张卡片，其余每8个穿一张卡片，每个数据后穿一逗点。

2.上机操作按机器要求的操作命令，先输入源程序，待编译结束，进行标准变量赋值(当需要时)，再输入算题数据，机器开始算题。

最后输出计算成果，计算结束。

本程序在TQ—16机安排标准变量赋值功能为：H00000=1多组数据自动输入，若H00000不赋1或H00000≠l，则无此功能。

H00001=q，q为衬砌组成构件(半跨顶拱、边墙、底板等)的等分计算段数，可取10或30，40，50等分段。

为提高计算精度，大断面隧洞可取较多的计算段数。

当H000001不赋值，或所赋值不等于10或30，40，50时，程序自动按等分20段执行。

H00002=1水头按均匀水头计算。

“本程序的FORTRAN程序，已有多组数据自动输入功能，分段数q值的改变，可用调换q值赋值语句卡片来实现。

等水头计算功能，可用改变编译卡，使程序中的注解成为计算语句来实现。

”3.计算成果程序列表输出各计算点的内力和位移，格式为：单位及符号：M——弯矩，以衬砌内壁受拉为正，t·m。

N——轴向力，拉为正，t。

Q——剪力，逆时针转动为正，t。

V——法向位移，以向外法线方向位移为正，m。

U——切向位移，顺外法线向右为正，m。

W——转角位移，逆时针为正，弧度。

输出是由衬砌顶部的中点开始，以半跨顶拱(或顶板)、边墙、侧拱(或底圆角)、半跨底板(或底拱)为序，每段按等分10小段11点输出结果。

4.算例圆拱直墙形无压隧洞，受顶部垂直山岩压力7t/m，齐墙顶内水压力及衬砌自重等荷载，断面尺寸如附图1。

并已知：围岩弹性抗力系数K=5×104t/m3；混凝土弹性模量E=23×102t/m2；混凝土容重W l=2.4t/m3；题号NO=198306(任选)；断面型号A=4。

附图1采用TQ—16机ALGOL程序原始数据：6*6*198306，4；∧∧……∧6*6*6*6*3，6，60，0.6，0.6，0.6，0.6，6，7，0，0，0，2.4，5410，5410，5410，23510；∧∧……∧6*6*输出成果；原始数据NO：198306A=43.00 6.0060.000.600.600.60 6.00 6.007.000.000.000.00 2.405000050000500002300000计算结果由计算成果正理绘制衬砌内力及抗力分布图，见附图2。

附图2算例内力及抗力分布图第二节源程序6*6*beginreal r1，r24，r4，fi1，fi2，fi，d1，d2，d3，d4，d5，d4b，e，ml，k1，k2，k3，k4，k5，h3，h，ld，ho，qn，e1，e2，qd，fo2，fn2，go3，gn3，eo4，en4，ds1，ds2，ds3，ds4，ds5，ab，ab3，ab4，hf，hf3，hf4，sf2，cf4，ho3，ho4，ho5，en2，eo2；in teger m1，m2，m3，m4，m5，m，a，b，no，q，t，s；switch sw1：=la1，la2，la3，la4，la5，la6，la7，la8，la9，la10，la11；m24：#read(0，‘N‘，no，a)；#print(0，‘80s-，2/，28X，4Hyuan，4X，3Hshi，4X，3Hshu，4X，2H u，/，28X，24s*，2/，28X，3Hno=，I8，9X，2Ha=，I2，/‘，no，a，)；m1∶=m2∶=m3∶=m4∶=m5∶=0；m∶=0.01745329；q∶=20；for s：10，20，30，40，50doif H00001=s then q∶=H00001；goto sw1[a]；la1：#read(0，‘N‘，r24，d2，ho，qn，e1，e2，qd，mL，k2，k4，e)；#print(0，‘8F10.2，3I10‘，r24，d2，ho，qn，e1，e2，qd，ml，k2，k4，e)；d4b∶=d4∶=d2；r4∶=r24；fi∶=fi2∶=90；fn2∶=e2；en4∶=fo2∶=(e1+e2)/2；ab∶=#abs(ho)；ho4∶=if ab＞r24then#sign(ho)*(ab﹣r24)else0；goto Loo；la2：#read(o，‘N‘，ld，fi，d4，d4b，ho，qn，e1，e2，ml，k4，e)；#print(0，‘9F10.2，2I10‘，ld，fi，d4，d4b，ho，qn，e1，e2，ml，k4，e)；en4∶=e2；ho4∶=ho；r4∶=Ld/#sin(fi*m)；goto Loola3：#read(0，‘N‘，b，ld，h，fi，d3，d4b，d4，ho，qn，e1，e2，ml，k3，k4，e)；#print(0，‘I10，11F10.2，3I10‘，b，ld，h，fi，d3，d4b，d4，ho，qn，e1，e2，ml，k3，k4，e)；m3∶=q；ds3：=h/m3；h3：=h；la31：r4：=ld/#sin(fi*m)；ho3：=ho；gn3：=e2；ab3：=#abs(ho3)；hf3：=#sign(ho3)；ab：=#abs(ho)；en4：=go3：=e2-(e2-e1)/(h+r4*(1-#cφs(fi*m)))*h；ho4：=if ab＞h then#sign(ho)*(ab-h)else o；goto Loo；la4：#read(0，‘N’，ld，h，fi，d1，d3，d4b，d4，ho，qn，e1，e2，qd，ml，k1，k3，k4，e)；#print(0，‘13F10.2，4I10‘，ld，h，fi，d1，d3，d4b，d4，ho，qn，e1，e2，qd，ml，k1，k3，k4，e)；goto La31；la5：#read(0，‘N‘，ld，h，fi，r24，d1，d3，d4b，d4，ho，qn，e1，e2，qd，ml，k1，k3，k4，e)；#print(o，‘14F10.2，4I10‘，ld，h，fi，r24，d1，d3，d4b，d4，ho，qn，e1，e2，qd，ml，k1，k3，k4，e)；d2：=(d1+d3)/2；k2：=(k1+k3)/2；r4：=(ld+r24)/#sin(fi*m)；en4：=go3：=e2-(e2-e1)/(r24+h+r4*(1-#cφs(fi*m)))*(r24+h)；gn3：=fo2：=e2-(e2-e1)/(r24+h+r4*(1-#cφs(fi*m)))*r24；fn2：=e2；ab：=#abs(ho)；hf：=#sign(ho)；fi2：=90；ho3：=if ab＞r24then hf*(ab-r24)else o；ho4：=if ab＞h+r24then hf*(ab-h-r24)else o；goto Loo；la6：#read(0，‘N‘，r24，fi2，r4，fi，d1，d2，d4，ho，qn，e1，e2，qd，ml，k1，k2，k4，e)；#print(0，‘13F10.2，4I10‘，r24，fi2，r4，fi，d1，d2，d4，ho，qn，e1，e2，qd，ml，k1，k2，k4，e)；m1：=q；d4b：=(d2+d4)/2；ld：=r4*(1﹣#cφs(fi*m))﹣r24*(1﹣#cφs(fi2*m))；ds1：=ld/m1；fi1：=o；la61：cf4：=#cφs(fi*m)；sf2：=#sin(fi2*m)；en4：=fo2：=e2﹣(e2﹣e1)/(r24*sf2+r4*(1﹣cf4)*r24*sf2；fn2：=e2；la62∶ab：=#abs(ho)；ho4∶=if ab＞r24*sf2then#sign(ho)*(ab﹣r24*sf2)else0；goto Loo；la7：#read(0，‘N‘，r1，fi1，r24，fi2，r4，fi，d1，d2，d4，ho，qn，e1，e2，qd，ml，k1，k2，k4，e)；#print(0，‘15F10.2，4I10‘r1，fi1，r24，fi2，r4，fi，d1，d2，d4，ho，qn，e1，e2，qd，ml，k1，k2，k4，e)；m1∶=q；fi1：=fi1*m；d4b：=(d2+d4)/2；ds1：=r1*fi1/*m1；goto La61；La8：#read(0，‘N‘，b，r24，fi2，r4，fi，d2，d4b，d4，ho，qn，e1，e2，ml，k2，k4，e)；#print(0，‘I10，12F10.2，3I10‘，b，r24，fi2，r4，fi2，d2，d4b，d4，ho，qn，e1，e2，ml，k2，k4，e)；if a=8then goto la61；sf2：=#sin(fi2*m)；en4：=eo2：=e2﹣(e2﹣e1)/((r24﹣r4)*sf2+r4)*r24*sf2；en2：=e2；goto la62；la10：#read(0，‘N‘，ld，r24，h，d1，d3，d5，ho，qn，e1，e2，qd，ml，k1，k3，k5，e)；#print(0，‘12F10.2，4I10‘，lb，r24，h，d1，d3，d5，ho，qn，e1，e2，qd，nl，k1，k3，k5，e)；m1：=m3：=m5：=q；d2：=(d1+d3)/2；d4b：=d3；d4：=d5；r4：=r24；fi：=fi2：=90；ds1：=ds5：=ld/m5；h3：=h；ds3：=h/m3；k2：=(k1+k3)/2；k4：=(k3+k5)/2；if ld=0then m1：=m5：=0；if h=0then m3：0，en4：=go3：=e1+(e2﹣e1)/(h+2*r24)*r24；fn2：=e2；gn3：=fo2：=e2﹣(e2﹣e1)/(h+2*r24)*r24；ab：=#abs(ho)；hf：=#sign(ho)；ho3：=if ab＞r24then hf*(ab﹣r24)else0；ho4：=if ab＞h+r24then hf(ab﹣r24﹣h)else0；ho5：=if ab＞h+2*r24then hf(ab﹣h﹣2*r24)else0；ab3：=#abs(ho3)；hf3：=#sign(ho3)；goto Loo；la11：#read(0，‘N‘，ld，h，d1，d3，d5，ho，qn，e1，e2，qd，ml，k1，k3，k5，e)；#print(0，‘11F10.2，4I10‘，lb，h，d1，d3，d5，ho，qn，e1，e2，qd，ml，k1，k3，k5，e)；ab：=#abs(ho)；m5：=q；ds5：=ld/m5；go3：=e1，gn3：=e2；ho3：=ho；ho5：=if ab＞h then#sign(ho)*(ab﹣h)else0；loo：if a≤10thenbegin m4：=q；fi：=fi*m；ds4：=fi*r4/m4；eo4：=e1；cf4：=#cφs(fi)；ab4：=#abs(ho4)；hf4：=#sign(ho4)；ena；if a=1∨(a≥5∧a≤10)thenbegin m2：=q；fi2：=fi2*m；ds2：=fi2*r24/m2；sf2：=#sin(fi2)；ab：=#abs(ho)；hf：=#sign(ho)；end；if a=4∨a=5∨a=11thenbegin m1：=m3：=q；ds1：=ld/m1；ds3：=h/m3；h3：=h；ab3：=#abs(ho3)；hf3：=#sign(ho3)；end；if a=10∧r24=0then m2：=m4：=0；m：=m1+m2+m3+m4+m5；beginreal s，so，ke，w，we，eps；integer i，j，k，mx；array c，d[1∶3，1∶6]；bt，a1，dl[1∶4]，c1，c2，c3，co，cl，g，gi，ae，an，aa，de[1∶6，1∶6]，xo，fe，p，cp，h，hi，b1，b2[1∶6]w1，wo[1∶251]，x[1∶1506]，y[1∶330]；procedu re mat1(a，b，c)；array a，b，c；begin integer i，j，k；for i：=1step1until6dofor k：=1step1until6dobegin c[i，k]：=0；for j：=1step1until6doc[i，k]：=c[i，k]+a[i，j]*b[j，k]；end；end mat1；procedure mat2(a，f，h，d，b)；v alue d；real d；array a，f，h，b；begin integer i，j；for i：=1step1until6dobegin b[i]：=0；for j：=1step1until6dob[i]：=b[i]+a[i，j]*f[j]；b[i]：=b[i]+h[i]*d；end；end mat2；procedure gg4(x，y)；array x，y；begin integer i，j；for i：=1step1until6dofor j：=1step1until6doy[i，j]：=x[i，j]；end gg4；procedure gg5(x，y)；array x，y；begin integer i；for i：=1step1until6doy[i]：=x[i]；end gg5；procedure gnhn；begin integer i，j，l；s：=so；for i：=1step1until6dofor j：=1step1until6dobeging[i，j]：=if i=j then1else0；gi[i，j]：=hi[i]：=h[i]：=0；end；for j：=1step1until4dobegin so：=s+dl[j]；begin real gx，gy，d，b1，b2，bo，hm，cf；switch sw2：=lm1，lm2，lm3，lm4，lm5；goto sw2[mx]；lm1：d：=d1；p[1]：=0；p[2]：=qd﹣ho﹣d*ml；goto ll；lm2：if a≠9thenbegin d：=d2；bo：=fi2﹣so/r24；b1：=#sin(bo)；b2：=#cφs(bo)；gx：=qd﹣d*ml/(if bo=0then10﹣4else b1)；gy：=fo2+b1*(fn2﹣fo2)/sf2；cf：=r24*(sf2﹣b1)；hm：=if ab＞cf then hf*(ab﹣cf)else0；lm21：if H00002=1then hm：=ho：p[1]：=(gy﹣gx)*b2*b1；p[2]：=gy*b2*b2+gx*b1*b1﹣hm；goto ll；and；bo：=so/r24；b1：=#sin(bo)；b2：=#cφs(bo)；d：=d4b+(fi2﹣bo)/fi2)↑2*(d2﹣d4b)；gx：=eo2+(sf2﹣b1)/sf2*(en2﹣eo2)；gy：=qn+d*ml/(if bo=0then10﹣4else b1)；cf：=r24*b1；hm：=if ab＞cf then hf*(ab﹣cf)else0；if H00002=1then hm：=ho：p[1]：=(gy﹣gx)*b1*b2；p[2]：=gy*b1*b1+gx*b2*b2﹣hm；goto ll；lm3：d：=d3；hm：=if ab3＞so then hf3*(ab3﹣so)else0；if H00002=1then hm：=ho；p[1]：=d*ml；p[2]：=(h3﹣so)/h3*(gn3﹣go3)+go3﹣hm；goto ll；lm4：bo：=fi﹣so/r4；b1：=#sin(ho)；b2：=cφs(bo)；d：=d4+(bo/fi)↑2*(d4b﹣d4)；gx：=eo4+(1﹣b2)*(en4﹣eo4)/(1﹣cf4)；gy：=qn+d*ml/(if bo＞1.57then10﹣4else b2)；cf：=r4*(b2﹣cf4)；hm：=if ab4＞cf then hf4*(ab4﹣cf)else0；goto lm21；lm5：d：=d5；p[1]：=0；if H00002=1then ho5：=ho；p[2]：=qn+d*ml﹣ho5；ll：an[2，5]：=if so=s then w*ke else we*ke；an[4，1]：=1/d；an[6，3]：=12/d↑3；end；for i：=1step1until6dofor l：=1step1until6dobegingi[i，l]：=if i=l then1+a1[j]*gi[i，l]else a1[j]*gi[i，l]；ae[i，l]：=dl[4]*an[i，l]；end；for i：=1step1until6dohi[i]：=a1[j]*hi[i]；mat1(ae，gi，cl)；gg4(cl，gi)；mat2(ae，hi，p，1，cp)；gg5(cp，hi)；for i：=1step1until6dobegin h[i]：=h[i]+bt[j]*hi[i]；for l：=1step1until6dog[i，l]：=g[i，l]，+bt[j]*gi[i，l]；end；end；end gnhn；procedure gj(n，a，b，ep，c1，c2)；value n，ep；real n，ep；arrag a，b，c1，c2；begininteger i，j，k，io，jo；real c，t；for i：=1step1until n doc1[i]：=i；for k：=1step1until n dobegin c：=0；for i：=k step1until n dofor j：=k step1until n doif#abs(a[i，j])＞#abs(c)thenbegin c：=a[i，j]；io：=i；jo：=jend；if#abs(c)≤ep thenbegin#print(0，‘10‘，c)；stopend；if jo≠k thenbeginfor i：=1step1until n dobegin t：=a[i，jo]；a[i，jo]：=a[i，k]；a[i，k]：=tend；j：=c1[k]；c1[k]：=c1[jo]；c1[jo]：=jend；if io≠k thenbeginfor j：=k step1until n dobegin t：=a[io，j]；a[io，j]：=a[k，j]；a[k，j]：=t end；t：=b[io]；b[io]：=b[k]；b[k]：=tend：c：=1/c；for j：=k+1step1until n doa[k，j]；=a[k，j]*c；b[k]：=b[k]*c；for i：=k+1step1until n dobeginfor j：=k+1step1until n doa[i，j]：=a[i，j]﹣a[i，k]*a[k，j]；b[i]：=b[i]﹣a[i，k]*b[k]end；end；for i：=n﹣1step﹣1until1dofor j：=i+1step1until n dob[i]：=b[i]﹣a[i，j]*b[j]；for k：=1step1until n doc2[c1[k]]：=b[k]；for k：=1step1until n dob[k]：=c2[k]；end gj；a1[1]：=a1[2]：=a1[3]：=0.5；a1[4]：=1；bt[1]：=bt[4]：=1/6；bt[2]：=bt[3]：=1/3；for j：=1step1until6dobegin p[j]：=0；for j：=1step1until3doc[i.j]：=d[i.j]：=0；for i：=1step1until6doc1[i，j]：=c2[i，j]：=c3[i，j]：=an[i，j]：=0；end；if a=2thenbegin c[1，1]：=c[2，3]：=c[3，5]：=1；c[1，4]：=﹣k4/e*d4b；c[2，6]：=c[1，4]*d4b↑2/12；goto r；and；if a=3∨a=8∨a=9thenbeginif b=1then c[1，3]：=c[2，4]：=c[3，5]：=1；if b=2then c[1，4]：=c[2，5]：=c[3，6]：=1；if b=3thenbegin c[1，1]：=c[2，3]：=c[3，5]：=1；c[1，4]：=﹣1.2*(if a=3then k3/e*d3else k2/e*d2)；c[2，6]：=c[1，4]*(if a=3then d3eles d2)↑2/12；end；end else c[1，2]：=c[2，4]：=c[3，6]=1；r：d[1，2]：=d[2，4]：=d[3，6]：=1；c1[1，2]：=c1[4，5]：=﹣1；c1[2，1]：=c1[3，3]：=c1[5，4]：=c1[6，6]：=1；if a≤2∨a≥9then goto l01；c2[1，1]：=c2[2，2]：=c2[4，4]：=c2[5，5]：=#sin(fi)；c2[1，2]：=c2[4，5]：=﹣#cφs(fi)；c2[2，1]：=c2[5，4]：=#cφs(fi)；c2[3，3]：=c2[6，6]：=1；if a=6∨a=7thenbegin c3[1，1]：=c3[2，2]：=c3[4，4]：=c3[5，5]：=#sin(fi2+fi1)；c3[1，2]：=c3[4，5]：=﹣#cφs(fi1+fi2)；c3[2，1]：=c3[5，4]：=#cφs(fi1+fi2)；c3[3，3]：=c3[6，6]：=1；end；lo1：an[3，2]：=﹣1；an[5，6]：=1；for i：=1step1until m+1dowo[i]：=1；b：=0；llo：s：=0；b：=b+1；for i：=1step1until6dobegin fe[i]：=0；for j：=1step1until6dode[i，j]：=if i=j then1etse0；end；if m1≠0then goto lo else goto l1；lo：so=0；dl[1]：=0；dl[4]：=ds1；dl[2]：=dl[3]：=dl[4]/2；ke：=k1/e；mx：=1；if a=7thenbegin an[2，1]：=an[5，4]：=1/r1；an[1，2]：=an[4，5]：=﹣1/r1；end；l1：for i：=1step1until m dobegin w：=wo[i]；we：=wo[i+1]；if i=m1+1∧m2≠0then goto l2；if i=m1+m2+1∧m3≠0then goto l3；if i=m1+m2+m3+1∧m4≠0then goto l4；if i=m1+m2+m3+m4+1∧m5≠0then goto l5；l11：gnhn；mat1(g，de，cl)；gg4(cl，de)；mat2(g，fe，h，dl[4]，cp)；gg5(cp，fe)；goto l6；l2：so：=0；dl[1]：=0；dl[4]：=ds2；ke：=k2/e；dl[2]：=dl[3]：=dl[4]/2；an[2，1]：=an[5，4]：=1/r24；an[1，2]：=an[4，5]：=﹣1/r24；mx：=2；if a＜6∨a＞7then goto l11；gnhn；mat1(g，c3，cl)；mat1(cl，de，co)；gg4(co，de)；mat2(cl，fe，h，dl[4]，cp)；gg5(cp，fe)；goto l6；l3：so：=0；dl[1]：=0；dl[4]：=ds3；mx：=3；dl[2]：=dl[3]：=dl[4]/2；ke：=k3/e；an[2，1]：=an[5，4]：=an[1，2]：=an[4，5]：=0；if m2=0∧m1≠0thenbeginl31：gnhn；mat1(g，c1，cl)；mat1(cl，de，co)；gg4(co，de)：mat2(cl，fe，h，dl[4]，cp)；gg5(cp，fe)；goto l6；end else goto l11；l4：so：=0；dl[4]：=ds4；dl[2]：=dl[3]：=dl[4]/2；dl[1]：=0，ke：=k4/e；mx：=4；an[2，1]：=an[5，4]：=1/r4；an[1，2]：=an[4，5]：=﹣1/r4；if a≤2∨a≥9then goto l11；gnhn；mat1(g，c2，cl)；mat1(cl，de，co)；gg4(co，de)；mat2(cl，fe，h，dl[4]，cp)；gg5(cp，fe)；goto l6；l5：so：=0；dl[4]：=ds5；dl[2]：=dl[3]：=dl[4]/2；ke：=k5/e；mx：=5；an[1，2]：=an[2，1]：=an[4，5]：=an[5，4]：=0；if m4=0then goto l31else goto l11；l6：end；for i：=1step1until3dobegin xo[i]：=xo[i+3]：=0；for k：=1step1until6dobegin aa[i，k]：=c[i，k]；aa[i+3，k]：=0；xo[i+3]：=xo[i+3]﹣d[i，k]*fe[k]；for j：=1step1until6doaa[i+3，k]：=aa[i+3，k]+d[i，j]*de[j，k]；end；end；gj(6，aa，xo，10﹣5，b1，b2)；for i：=1step1until6dox[i]：=xo[i]；w1[1]：=if xo[5]＞0then1else0；if m1≠0then goto l7else goto l8；l7：so：=0；dl[4]：=ds1；dl[2]：=dl[3]：=dl[4]/2；ke：=k1/e；an[1，2]：=an[2，1]：=an[4，5]：=an[5，4]：=0；if a=7thenbegin an[2，1]：=an[5，4]：=1/r1；an[1，2]：=an[4，5]：=﹣1/r1；end；mx：=1；l8：for i：=1step1until m dobegin w：=wo[i]；we：=wo[i+1]；if i=m1+1∧m2≠0then goto l12；if i=m1+m2+1∧m3≠0then goto l13；if i=m1+m2+m3+1∧m4≠0then goto l14；if i=m1+m2+m3+m4+1∧m5≠0then goto l15；l111：gnhn；mat2(g，xo，h，dl[4]，cp)；for j：=1step1until6dox[6*i+j]：=xo[j]：=cp[j]；goto l16；l12：so：=0；dl[4]：=ds2；dl[2]：=dl[3]：=dl[4]/2；an[2，1]：=an[5，4]：=1/r24；ke：=k2/e；an[1，2]：=an[4，5]：=﹣1/r24；mx：=2；if a＜6∨a＞7then goto l111；gnhn；mat1(g，c3，cl)；goto l132；l13：so：=0；dl[1]：=0；dl[4]：=ds3；dl[2]：=dl[3]：=dl[4]/2，ke：=k3/e；mx：=3；an[1，2]：=an[2，1]：=an[4，5]：=an[5，4]：=0；if m2=0∧bm1≠0thenbeginl131：gnhn；mat1(g，c1，cl)l132：mat2(cl，xo，h，dl[4]，cp)；for j：=1step1until6dox[6*i+j]：=xo[j]：=cp[j]；gotol16；end else goto l111；l14：so=0；dl[4]：=ds4；dl[2]：=dl[3]：=dl[4]/2；ke：=k4/e；mx：=4；an[2，1]：=an[5，4]：=1/r4；an[1，2]：=an[4，5]：=﹣1/r4；if a≤2∨a≥9then goto l111；gnhn；mat1(g，c2，cl)；mat2(cl，xo，h，dl[4]，cp)；for j：=1step1until6dox[6*i+j]：=xo[j]：=cp[j]；goto l16；l15：so：=0；dl[4]：=ds5；dl[2]：=dl[3]：=dl[4]/2；an[1，2]：=an[2，1]：=an[4，5]：=an[5，4]：=0；ke：=k5/e；mx：=5；if m4=0hen goto l131elsegoto l111；l16：w1[i+1]：=if xo[5]＞0then1else0；end；eps：=0；for j：=1step1until m+1doeps：=eps+(w1[i]﹣wo[i])↑2；if eps≥2then goto l9else goto l10；l9：for i：=1step1until m+1dowo[i]：=w1[i]；an[1，2]：=an[2，1]，an[4，5]：=an[5，4]：=0；goto ll0；l10：for i：=0step1until m dofor j：=4step1until6dox[6*i+j]：=x[6*i+j]/e；k：=q；t：=q/10；mx：=0；for i：=0step1until m dobeginfor j：=1step1until6dobegin mx：=mx+1；y[mx]：=x[6*i+j]；end；if i=k∧i≠m thenbegin k：=k+q；if i=m1∧m2≠0∧(a=6∨a=7)thenbegin gg4(c3，co)；goto ll；end；if i=m1+m2∧m3≠0∧(a=4∨a=11)thenlc1：begin gg4(c1，co)；goto ll；end；if i=m1+m2+m3∧m4≠0∧a≥3∧a≤8thenbegin gg4(c2，co)；goto ll；end；if i=m1+m2+m3+m4∧m5≠0∧a=11thengoto lc1；for s：=1step1until6dofor j：=1step1until6doco[s，j]：=if s=j then1else0；ll：for j：=1step1until6doxo[j]：=x[6*i+j]；mat2(co，xo，b1，0，b2)；for j：=1step1until6dobegin mx：=mx+1；y[mx]：=b2[j]；end；end；end；#print(0，‘2/，28X，2Hji，4X，4Hsuan；4X，3Hjie，3X，4Hgua0，/，28X，24s*，2/，2X，2Hmx，9X，Hm，11X，Hn，11X，Hq，11X，Hv，11X，Hu，11X，Hw，/‘)；begin array x1，x2[1∶3]；k：=0；q=10；for i：=mx/6step﹣1until1dobegin t：=6*(i﹣1)；x1[1]：=y[t+3]；x1[2]：=y[t+1]；x1[3]：=y[t+2]；x2[1]：=y[t+5]；x2[2]：=y[t+4]；x2[3]：=y[t+6]；#print(0，‘I4，F12.3，2X，F12.7‘.k，x1，x2)；if k=q thenbegin q：=q+10；#print(0，‘/‘)；goto lpp；end；k：=k+1；lpp：end；end；#print(0，‘I6‘，b)；#print(0，‘2/，80s﹣‘)；end；if H00000=1then goto m24；end∧∧……∧6*6*。