这是个用渐开线精确画直齿圆柱齿轮的例子

xcb

这个例子里面包含几个子程序，子程序这两天放上来。

原帖地址：http://www.ugufun.com/?p=77

$$ 定义参数：-----------------------------------------------------------------

NUMBER/RESP,R,RB,RA,RF,B,I,HA,HF,MAT(12),COL,CEN(3),P,FD,ZMIN
ENTITY/RFC,SPLN1,SLL1,SL1,SPLN2,SLL2,SL2,LN1,ARCC,SOLIDZ,SOLID(200),CSYS0
! G1 l9 X6 N, oSTRING/DIR(50),FNAME(50),TEMP(100),COLO(10),JTYPES(10)

$$ ---------------------------------------------------------------------------
$$ 取文件目录和文件名---------------------------------------------------------

DIR=&PSPEC $$ 取出文件目录
9 S  y, O6 C7 V0 C0 E4 z. \) bANS=LENF(DIR)
7 l! \, A3 Y; R# i3 W, aFNAME=SUBSTR(DIR,1,ANS-4,)+'.txt' $$ 取出文件名

$$ ---------------------------------------------------------------------------
$$ 默认值---------------------------------------------------------------------

M=2.5 $$ 模数
; V; L+ l, [+ F- B+ S/ v  WZ=20 $$ 齿数
A=20 $$ 压力角
TB=15 $$ 齿宽
+ r" d' Y( ^$ z: b" A$ \8 GHAC=1 $$ 齿顶高系数
9 ^% |( [# s8 T3 }+ P2 t0 ?CC=0.25 $$ 顶隙系数
XX=0 $$ 变位系数

$$ ---------------------------------------------------------------------------
$$ 定义界面-------------------------------------------------------------------
0 I9 L' z2 k$ tL10:
+ L3 c: Q) l  gPARAM/'请输入渐开线直齿圆柱齿轮的各项基本参数:','模数',M,'齿数',INT,Z,$
'压力角',A,'齿宽',TB,'齿顶高系数',HAC,'齿侧间隙系数',CC,'变位系数',$
XX,RESP
JUMP/L10:,TERM:,,,RESP

L20:
CHOOSE/'请选择齿轮键的形式:','平键','花键','无',RESP
JUMP/L10:,TERM:,,,L30:,L40:,L50:,RESP
! [* a8 v: X+ G* HL30:
JTYPE=1
JUMP/L60:
6 m' n- ^6 M  B' WL40:
JTYPE=2
; H7 K+ H: ~  G: _" IJUMP/L60:
/ M; D- z  X# [8 K( sL50:
JTYPE=0
" G& i2 f" u( ^7 D0 I% x% \, J1 gL60:

$$ ---------------------------------------------------------------------------
$$ 基本几何参数计算-----------------------------------------------------------

PI=3.1415926 $$ 圆周率
5 w# \3 l$ u7 ]  j: L: t! jR=M*Z/2 $$ 分度圆半径
5 [$ R4 r7 K: N: YRB=R*COSF(A) $$ 基圆半径
' e7 W8 m" R9 N1 ?" xHA=M*(HAC+XX) $$ 齿顶高
2 Z4 X8 B5 z& Y: W, J* L, YHF=M*(HAC+CC-XX) $$ 齿根高
2 K6 \. c/ Z+ @! WRA=R+HA $$ 齿顶圆半径
RF=R-HF $$ 齿根圆半径
5 D1 N' d' \$ l; X0 [' p/ KP=PI*M $$ 周节
$ \. R, z$ F1 s5 x" D4 J1 IPB=PI*M*COSF(A) $$ 基节
+ H" u. \1 t0 }3 r+ Y5 S+ dSS=PI*M/2+2*XX*M*SINF(A)/COSF(A) $$ 分度圆弧齿厚
. E+ f9 Q2 s, o; b( @' yFD=HA+HF $$ 全齿高
ZMIN=INTF(2*HAC/(SINF(A)*SINF(A)))+1 $$ 最小齿数
XMIN=HAC-Z*SINF(A)*SINF(A)/2 $$ 最小变位系数

$$ ---------------------------------------------------------------------------
4 L7 g! }7 m) n4 Y( {& S/ X$$ 画出齿轮-------------------------------------------------------------------

CALL/'SUB1',A,XX,M,R,RB,HA,HF,RA,RF,Z,RFC,SPLN1,SLL1,SL1,ARCC,SPLN2,SLL2,$
SL2,LN1 $$ 调用渐开线子程序

SOLIDZ=SOLEXT/RFC,HEIGHT,TB,AXIS,0,0,1 $$ 拉伸齿根圆柱
! ^* x: E9 o) }# bSOLID(1)=SOLEXT/SPLN1,ARCC,SPLN2,SLL2,SL2,LN1,SL1,SLL1,HEIGHT,TB,AXIS,$
0 |' z' W' `) y8 T! P1 E0,0,1 $$ 拉伸轮齿
MAT=MATRIX/XYROT,360/Z
/ l  I$ d( w7 _7 |I=1
LLL:
IFTHEN/I<Z
SOLID(I+1)=TRANSF/MAT,SOLID(I) $$ 旋转拷贝轮齿
( j9 F. a+ n. ZI=I+1
JUMP/LLL:
ENDIF

$$ ---------------------------------------------------------------------------
9 n2 b: d) ?& {4 u( W5 X5 E2 X$$ 画出键---------------------------------------------------------------------

IFTHEN/JTYPE==1
7 F7 V$ h2 V( p7 wCALL/'SUB5',TB,SOLIDZ $$ 调用平键子程序
8 I# m: p+ {' U" ]/ [; qELSEIF/JTYPE==2
0 }7 v: Y5 P3 p! Z  G0 T. \/ d. GCALL/'SUB11',TB,SOLIDZ $$ 调用花键子程序
ELSE
5 H- n$ _( s4 V) X- o" ZENDIF

$$ ---------------------------------------------------------------------------
$$ 建立齿轮信息文件-----------------------------------------------------------

COL=&COLOR(SOLIDZ) $$ 取出实体颜色
. P' U: [% P) r5 C& m5 y5 CCSYS0=&WCS
- e, g& V0 U. m$ j; ~# @- n&WCS=&ABS
CEN=&CENTER(RFC) $$ 取出齿轮位置中心
8 e7 c" u1 ^( y) W( ~  C; E&WCS=CSYS0

IFTHEN/JTYPE==0
( _! `: l, F% C# ^& AJTYPES='无'
+ D& ~) j/ {' r; U* q4 iELSEIF/JTYPE==1
JTYPES='平键'
ELSE
& R1 r: F) _1 t) s5 M! v5 `+ oJTYPES='花键'
ENDIF

IFTHEN/COL==1
COLO='BLUE'
: ]# \8 X0 ^4 M) `ELSEIF/COL==2
+ b) q; l8 o& s" cCOLO='GREEN'
+ g. |- T6 W/ ~. L) Q# K( jELSEIF/COL==3
COLO='CYAN'
ELSEIF/COL==4
4 Z/ V* u8 H' x4 t% h- LCOLO='RED'
ELSEIF/COL==5
COLO='MAGENT'
ELSEIF/COL==6
COLO='YELLOW'
ELSEIF/COL==7
COLO='WHITE'
' E1 s) a. b  qELSEIF/COL==8
COLO='OLIVE'
$ j. ?, ~0 l% T5 w' [ELSEIF/COL==9
COLO='PINK'
% @' N' ~% K: `% lELSEIF/COL==10
% C; r, F' v- B  {COLO='BROWN'
5 H* K- ?& P/ YELSEIF/COL==11
COLO='ORANGE'
! Y2 }6 |+ A& \2 WELSEIF/COL==12
& ~! x3 `6 z/ ^* |9 S# SCOLO='PURPLE'
0 V: I9 Z- ^' |# K* j8 w- ~ELSEIF/COL==13
3 n1 D( f/ ^7 o. PCOLO='DKRED'
ELSEIF/COL==14
COLO='AQUAMR'
. \" Z9 [0 x& U" j) z: `, aELSE
% S9 O, i+ i, x$ ]5 iCOLO='GRAY'
( S- C  [8 `5 \3 R/ tENDIF

FETCH/TXT,1,FNAME,IFERR,FL10:
1 ?' f1 q. E. I0 V4 w! {: jRESET/1
" W7 h# g% u- g# mNUM1=1
5 w8 G, D0 E: @! ^/ e. F) P' aFL40:
READ/1,IFEND,FL20:,TEMP
IFTHEN/TEMP=='#END#'
NUM1=NUM1+1
ENDIF
, ^* W2 ~! s6 b0 RJUMP/FL40:

FL10:
; J  x, B) r' H0 x  ~+ iCREATE/TXT,1,FNAME
$ n7 v5 H' D* |4 K4 N- Q# iNUM1=1

FL20:
WRITE/1,'----------------------------------------------------------------- '
5 C. [8 F5 n2 _0 t" E2 WWRITE/1,'有问题，请到http://ugapi.com 论坛上提问! '
WRITE/1,' '
4 N1 F2 ]' E8 |# NWRITE/1,'齿轮编号： ('+ISTR(NUM1)+') :'
$ r' W1 Q& H  w2 M% r* s: zWRITE/1,' '
7 v1 P0 f2 r9 ^9 ZWRITE/1,' 齿轮形式 : '+'渐开线直齿圆柱齿轮'
: S% k+ [" b. ZWRITE/1,' 齿轮信息文件 : '+FNAME
( a( q) q7 f9 J: M- {; PWRITE/1,' 齿轮颜色 : '+COLO
7 A. |; c7 {; C) `* NWRITE/1,' 齿轮键的形式 : '+JTYPES
WRITE/1,' 齿轮中心坐标 ( X,Y,Z) : '
WRITE/1,' '+FSTR(CEN(1))+' ,'+FSTR(CEN(2))+$
2 Z7 s% H) f6 t( H' {$ G' ,'+FSTR(CEN(3))
3 W( r* y( l) W% }) o4 F2 b: [$ IWRITE/1,' '
WRITE/1,' 齿轮几何参数 :'
) l& Y# ~$ m6 z2 \$ YWRITE/1,' 模数 : '+FSTR(M)
WRITE/1,' 齿数 : '+ISTR(Z)
- ]- A2 Z/ |% _) ?0 EWRITE/1,' 压力角 : '+FSTR(A)
WRITE/1,' 齿宽 : '+FSTR(TB)
9 b& g; q' F3 h4 @3 ^* z, _/ X  h5 N* mWRITE/1,' 齿顶高系数 : '+FSTR(HAC)
WRITE/1,' 齿侧间隙系数 : '+FSTR(CC)
4 E5 R- e6 j2 yWRITE/1,' 变位系数 : '+FSTR(XX)
WRITE/1,' '
* z8 i" ?2 u8 ~! B- IWRITE/1,' 分度圆直径 : '+FSTR(2*R)
WRITE/1,' 基圆直径 : '+FSTR(2*RB)
. [1 I0 S/ L: b/ Q! mWRITE/1,' 齿顶高 : '+FSTR(HA)
2 Z- b4 i! K2 QWRITE/1,' 齿根高 : '+FSTR(HF)
WRITE/1,' 齿顶圆直径 : '+FSTR(2*RA)
2 W  q! O# k  ?8 a3 G3 L! PWRITE/1,' 齿根圆直径 : '+FSTR(2*RF)
" v% q( R8 p2 lWRITE/1,' 周节 : '+FSTR(P)
6 b( d1 ?9 @8 r) d( c0 _WRITE/1,' 基节 : '+FSTR(PB)
; c; S; O; Z# S: i2 T# aWRITE/1,' 分度圆弧齿厚 : '+FSTR(SS)
" f, t$ E  W% ^7 ]6 w' S9 I  hWRITE/1,' 全齿高 : '+FSTR(FD)
WRITE/1,' 最小齿数 : '+ISTR(ZMIN)
WRITE/1,' 最小变位系数 : '+FSTR(XMIN)
+ D  K0 ^4 f" |9 U; qWRITE/1,' '
WRITE/1,'#END#'

FILE/TXT,1
RESET/1
FPRINT/1 $$ 将齿轮信息显示在窗口上
+ w! a( z4 o2 K4 d$ q: {8 UFTERM/TXT,1

$$ ---------------------------------------------------------------------------
&WCS=CSYS0
TERM:
& K2 |6 h/ Z3 FHALT
, I+ J8 N' W8 y9 o0 y5 R9 e$$ 建模结束-------------------------------------------------------------------

ab1025950132

楼主直接把二次开发的程序给大家共享吧！

xcb

$$ 渐开线子程序_sub1

7 X4 r$ ]0 ^7 a# a1 m6 O6 b9 Q- j
$$ 传递参数-------------------------------------------------------------------

4 R% M: }  I" |8 B$ a1 D7 X0 hPROC/A,XX,M,R,RB,HA,HF,RA,RF,Z,RFC,SPLN1,SLL1,SL1,ARCC,SPLN2,SLL2,SL2,LN1
$ q& K8 @8 b* Z- t& z8 E
$ y* T& Q/ h- i+ j# e$$ ---------------------------------------------------------------------------


$ g" D5 t* r& R) y$$ 定义参数-------------------------------------------------------------------
5 p7 Q# }# J, [8 L. h, c/ l
NUMBER/A,XX,M,R,RB,RA,RF,B,I,HA,HF,Z,X(200),Y(200),MAT(12),NU,K1,K2
" {3 j+ M& p" B9 t5 C4 j6 S4 ]ENTITY/RFC,SPLN1,SL1,SPLN2,SL2,LN1,ARCC,P(200),PT(2),PTT(2),SLL1,SLL2

$$ ---------------------------------------------------------------------------

; \& r' R9 P+ r2 M
$$ 画出渐开线-----------------------------------------------------------------
+ |* T& v3 P$ q! E) {7 \
PI=3.1415926
  n: i" U  \  xB=0
I=1
: f7 ~$ k6 w- aRFC=CIRCLE/0,0,RF      $$ 画齿顶圆
) m' s" r' T( h- W7 m7 {L10:
6 y3 l# P; {9 \! RX(I)=RB*COSF((SINF(B)/COSF(B))*180/PI-B)/COSF(B)  
Y(I)=RB*SINF((SINF(B)/COSF(B))*180/PI-B)/COSF(B)  

IFTHEN/(RB/COSF(B))>RA
. P5 ~0 I4 ]* u( C6 v" i1 W& ~9 SJUMP/L20:
ELSE
: X9 N6 _- M3 F; v3 {! GP(I)=POINT/X(I),Y(I)   $$ 渐开线上的点
I=I+1
, h' P/ H" E9 ~8 A& yB=B+0.5
7 A: q4 m7 ]& {8 @1 g0 d9 uENDIF
JUMP/L10:

L20:
B=ACOSF(RB/RA)
X(I)=RB*COSF((SINF(B)/COSF(B))*180/PI-B)/COSF(B)
  i8 U: T4 M# S, ?+ W+ E4 OY(I)=RB*SINF((SINF(B)/COSF(B))*180/PI-B)/COSF(B)
# h# x( o8 k0 N2 k* x+ `* aP(I)=POINT/X(I),Y(I)
SPLN1=SPLINE/P(1..I)   $$ 轮齿一侧渐开线
( `7 G/ M) |4 W$ I, P7 M. U
IFTHEN/RB>RF
     PT(1)=POINT/RB,-(RB*RB-RF*RF)/(2*RF)
     SL1=FILLET/SPLN1,RFC,CENTER,PT(1),RADIUS,$
               (RB*RB-RF*RF)/(2*RF),NOTRIM
ELSE
     NU1=SQRTF((RF+0.38*M)**2-RB**2)
     NU2=ATANF((NU1-0.38*M)/RB)
     NU3=(NU1-0.38*M)/RB*180/PI-NU2
/ H, \* |: _  T2 `     NU4=NU2+NU3-ACOSF(RB/(RF+0.38*M))
     PT(1)=POINT/(RF+0.38*M)*COSF(NU4),(RF+0.38*M)*SINF(NU4)
     SL1=FILLET/SPLN1,RFC,CENTER,PT(1),RADIUS,0.38*M
     DELETE/RFC
     RFC=CIRCLE/0,0,RF
+ i! f) F' A) |, C) S* p/ w9 aENDIF                  $$ 轮齿一侧倒角圆
2 I+ g' F6 [! h6 v
& \0 Y1 X2 {; G8 u3 yPTT(1)=POINT/ENDOF,XSMALL,SPLN1
$ }( B. o: T8 z1 |3 o% kPTT(2)=POINT/ENDOF,XLARGE,SL1
, `7 y7 w; K/ G1 W  XSLL1=LINE/PTT(1),PTT(2)$$ 轮齿一侧倒角圆与渐开线的连线
! B6 f% T, l5 D+ U
B=(SINF(ACOSF(RB/R))/COSF(ACOSF(RB/R)))*180/PI-ACOSF(RB/R)+$
# |; Y% R  }$ _! J5 T- o  90/Z+(2*XX*SINF(A)/COSF(A))*180/PI/Z
DELETE/PT(1)
5 h3 ?0 l3 q$ jPT(1)=POINT/0,0
! P* T. a- M' i0 b: T( RLN1=LINE/PT(1),ATANGL,B
B=2*B-((SINF(ACOSF(RB/RA))/COSF(ACOSF(RB/RA)))*180/PI-ACOSF(RB/RA))
K1=ACOSF(RB/RA)
0 {$ e( T* k8 qK2=(SINF(K1)/COSF(K1))*180/PI-K1
ARCC=CIRCLE/0,0,RA,START,K2,END,B
$ y. k  z; m& ^' }0 |; q8 KMAT=MATRIX/MIRROR,LN1
9 I% e5 Z  d1 O1 d# JSPLN2=TRANSF/MAT,SPLN1 $$ 轮齿另一侧渐开线
7 W' Q( Z6 L' C3 U' c4 lSL2=TRANSF/MAT,SL1     $$ 轮齿另一侧倒角圆
SLL2=TRANSF/MAT,SLL1   $$ 轮齿另一侧倒角圆与渐开线的连线
DELETE/LN1,P(1..I),PT(1)
$ O! q, n4 ?, I2 a  V2 d( |; x
0 V! K7 y1 C& z$ p; U8 k1 M: BPT(1)=POINT/ENDOF,YLARGE,SL2
& R6 I. n$ D  N- Z3 V8 b* Q% ZPT(2)=POINT/ENDOF,YSMALL,SL1
LN1=LINE/PT(1),PT(2)   $$ 齿根连线
! y4 z: Y1 L) }& M8 O( ODELETE/PT(1..2)
3 A0 f1 o: v# ^) U! M9 }3 y+ l$ HDELETE/PTT(1..2)
4 h" f3 C1 J: {% o( f2 \* TRETURN

$$ 返回-----------------------------------------------------------------------

xcb

$$ 渐开线子程序_sub5

. n4 v( |2 a# j7 K' y" s1 @! A/ d$$ 传递参数-------------------------------------------------------------------

PROC/TB,SOLIDZ

$$ ---------------------------------------------------------------------------

; y% P3 @) o/ c& v3 f5 B
( B* ?/ w9 {) k9 u: i3 ~$$ 定义参数-------------------------------------------------------------------

0 E+ w0 B! y" t* O/ XNUMBER/G1(3),MAT(12),RESP
+ M( M. g$ Y  P: P# PENTITY/CRH,OBJJ,OBJJ1,PTA,PTB,PTC,LNA,LNB,LNC,SOLIDD,SOLIDZ,SOLIDE
- d3 L" ^) f% n/ s* ^$ X# T* h
) {/ |; d5 @: k/ }* G: x* v$$ ---------------------------------------------------------------------------
( e: p$ `/ f. N( W' U7 z. C7 H

$$ 输入界面参数---------------------------------------------------------------
/ E. y# e1 E5 O9 ~4 t9 S' H1 K
0 W/ w  L3 Y" a' u# p: U5 FL70:
PARAM/'请输入键参数.','圆半径',RH,$
       '键高',HH,'键宽',HR,RESP
' S0 O& O+ G- |  l( |$ \JUMP/L70:,L80:,,,RESP
* C, |1 q( M6 e! T+ E+ D6 w( B% P1 g
$$ ---------------------------------------------------------------------------


$$ 画出键---------------------------------------------------------------------
( B8 Q: q8 W2 h$ D
0 i- x% E, {# K" d3 Q( r# R9 \BT=ASINF(HR/(2*RH))
: K1 P4 M4 w  A0 sCRH=CIRCLE/0,0,RH,START,180,END,360-BT
PTA=POINT/ENDOF,XLARGE,CRH
# l6 n/ E/ [3 Y. `, {G1=&POINT(PTA)
PTB=POINT/G1(1)+HH,G1(2)
LNA=LINE/PTA,PTB
PTC=POINT/G1(1)+HH,0
LNB=LINE/PTB,PTC
OBJJ=SPLINE/APPROX,DELETE,CRH,LNA,LNB
LNC=LINE/0,0,10,0
MAT=MATRIX/MIRROR,LNC
OBJJ1=TRANSF/MAT,OBJJ
  Z; _6 |$ }* |: {; l% ]2 D7 z0 tSOLIDD=SOLEXT/OBJJ,OBJJ1,HEIGHT,TB,AXIS,0,0,1
SOLIDE=SUBTRA/SOLIDZ,WITH,SOLIDD
3 s5 [  r1 \/ v+ R" \
L80:
RETURN
$ S- C0 b) ^$ X0 U2 e, V% X1 [
$$ 返回-----------------------------------------------------------------------

xcb

$$ 渐开线子程序_sub11
5 l6 j/ |$ T! N. \  r5 g! w! P# P8 B
$$ 传递参数-------------------------------------------------------------------

1 ?% {& X; R. V  K3 L. ?0 d2 ~& [/ aPROC/TB,SOLIDZ

9 K: G% A9 R: v$$ ---------------------------------------------------------------------------
- }: F9 g6 V+ a# I/ a

$$ ---------------------------------------------------------------------------

NUMBER/RESP,R,RB,RA,RF,B,I,HA,HF,MAT(12),BC,COL,CEN(3),P,FD,ZMIN
ENTITY/RFC,SPLN1,SLL1,SL1,SPLN2,SLL2,SL2,LN1,ARCC,SOLIDZ,SOLID(200),CSYS0,$
       SOLIDK,SOLIDL,SOLIDM

; T; m- A% V8 ]% Y$$ ---------------------------------------------------------------------------


# @7 E# A6 k3 |4 M$ c$$ 输入界面参数---------------------------------------------------------------

) Y& I3 x  _; G5 Y9 X9 LL10:
8 X# ~; a  E% z. c  ~PARAM/'请输入花键参数 ','模数',M,'齿数',INT,Z,'压力角',A,'齿顶高系数',HAC,$
2 f' x+ }2 l. w: \( x0 D      '齿侧间隙系数',CC,'变位系数',XX,RESP
JUMP/L10:,TERM:,,,RESP
! r6 ]- a- n8 Y7 }+ t( Y
. X6 z- M: {% }1 ^$$ ---------------------------------------------------------------------------


$$ 画出花键-------------------------------------------------------------------
! a" t1 h* m# X& Z3 w8 {# E. v0 C
7 }+ @0 F9 v& N" j: J1 VPI=3.1415926
R=M*Z/2
( `2 J5 |6 I5 S( h% \8 Z3 [# kRB=R*COSF(A)
HA=M*(HAC+XX)
HF=M*(HAC+CC-XX)
$ l, c$ E1 [: R8 \6 k3 f' i& [8 `& z' RRA=R+HA
# {7 _! ]2 T5 m( oRF=R-HF
% p# O' \: u" W$ P9 {# wBC=0
P=PI*M
FD=HA+HF
6 h! _4 y# R- H: E* \4 cZMIN=INTF(2*HAC/(SINF(A)*SINF(A)))+1
0 I( D% e5 P) s5 Q, H0 h3 `
CALL/'SUB1',A,XX,M,R,RB,HA,HF,RA,RF,Z,RFC,SPLN1,SLL1,SL1,ARCC,SPLN2,SLL2,$
$ ]- s! M. q4 l  W0 P4 M     SL2,LN1          $$ 调用渐开线子程序
( ~  ~# w1 I# S7 B) {
, |- i, V6 B# G3 x/ ISOLIDM=SOLEXT/RFC,HEIGHT,TB,AXIS,0,0,1
SOLID(1)=SOLEXT/SPLN1,ARCC,SPLN2,SLL2,SL2,LN1,SL1,SLL1,HEIGHT,TB,AXIS,0,0,1
MAT=MATRIX/XYROT,360/Z
! R8 o  {7 X& v$ lJ=1
LLL:
IFTHEN/J<Z
* I& j" J& P6 v* t( K/ \; vSOLID(J+1)=TRANSF/MAT,SOLID(J)
1 r: ^' U  K2 b! v) r  W! hJ=J+1
JUMP/LLL:
ENDIF
. T+ Q% a; v4 h
& z7 k2 _1 x+ e0 Z5 xSOLIDK=SUBTRA/SOLIDZ,WITH,SOLIDM
3 j( f& s4 t" h- f! T2 z( Q( TSOLIDL=SUBTRA/SOLIDK,WITH,SOLID(1..J)
* g! C! x; l" ]- ?' G) G) @# fTERM:
: Y0 F" @) j6 ]% j: Q% e1 x; Z& |9 SRETURN
8 A3 E! R2 A. ^: b0 H2 G
# h6 k+ G0 c5 o$$ 返回-----------------------------------------------------------------------

xcb

这个是编译好的程序：大家可以试一下。

" f" K% T( ^1 y7 Y( E. `; u