大地坐標轉(zhuǎn)換成施工坐標公式

大 地 ( 高 斯 平 面 ) 坐 標 系 工 程 坐 標 系 轉(zhuǎn) 換大地坐標系---> 工程坐標系========================待轉(zhuǎn)換點為 P，大地坐標為： Xp、Yp工程坐標系原點 o: 大地坐標：Xo、Yo工程坐標：xo、yo工程坐標系 x軸之大地方位角： adX=Xp-XodY=Yp-YoP點轉(zhuǎn)換后之工程坐標為 xp、yp:xp=dX*COS(a)+dY*SIN(a)+xoyp=-dX*SIN(a)+dY*COS(a)+yo工程坐標系---> 大地坐標系========================待轉(zhuǎn)換點為 P，工程坐標為： xp、yp工程坐標系原點 o: 大地坐標：Xo、Yo工程坐標：xo、yo工程坐標系x軸之大地方位角：adx=xp-xody=yp-yoP點轉(zhuǎn)換后之工程坐標為xp、yp:xp=Xo+dx*COS(a)-dy*SIN(a)yp=Yo+dx*SIN(a)+dy*COS(a)坐標方位角計算程序置鏡點坐標：ZXZY后視點坐標：HXHY方位角：W兩點間距離:SLb10←｛A,B,C,D｝←A〝ZX=〞:B〝ZY=〞:C〝HX=〞:D〝HY=〞:W=tg1((D-B)÷(C-A)):(D-B)>0=>(C-A)>0=>W=W:∟∟(D-B)>0=>(C-A)<0=>W=W+180:∟∟(D-B)<0=>(C-A)<0=>W=W+180:∟∟(D-B)<0=>(C-A)>0=>W=360+W∟∟W=W◢S=√((D-B)2+(C-A)2)◢Goto0←CASIOfx－4500p坐標計算程序根據(jù)坐標計算方位角W＝W＋360△W：“ALF(1~2)＝”L1A“X1＝”：B“Y1＝”：Pol(C“X2－”A，D“Y2”－B：“S＝”▲W＜0直線段坐標計算L1X“X(0)：”Y“Y(0)：”S“S(0)：”A“ALF”L2Lb12L3{L}：L“LX”L4M“X(Z)”＝X＋(L－S)cosA▲L5 N“Y(Z)＝”Y＋(L－S)sinA▲L6 {B} ：B“B(L)：”Q“Q”L7 O“X(L)＝”M＋Bcos(A＋Q＋180)▲L8 P“Y(L)＝”N＋Bsin(A＋Q＋180)▲L9 {C} ：C“B(R)”L10 U“X(R)＝”M＋Ccos(A＋Q)▲L11 V“Y(R)＝”N＋Csin(A＋Q)▲L12Goto2園曲線段坐標計算L1S“S(0)-Km”：X“X(0)：”Y“Y(0)：”A“ALF：”R“R：”K“K(L＝1,R＝2)”L2Lb12L3{L}：L“L(X)”L4V＝180／π×（L－S）／R：W＝V／2L5C＝A＋(-1)K×W：D＝2RsinW：F＝A＋(-1)K×VL6M“X(Z)＝”X＋DcosC▲L7N“Y(Z)＝”Y＋DsinC▲L8{E}：E“B(L)：”Q“Q”L9O“X(L)＝”M＋Ecos(F＋Q＋180)▲L10 P“Y(L)＝”N＋Esin(F＋Q＋180)▲L11{G} ：G“B(R)”L12 T“X(R)＝”M＋Gcos(F＋Q)▲L13 U“Y(R)＝”N＋Gsin(F＋Q)▲L14Goto2正向緩和曲線段坐標計算L1S“ZH-Km”：X“X(ZH)：”Y“Y(ZH)：”A“ALF：”R“R：”H“LS：”K“K(L＝1,R＝2)”L2Lb12L3{L}：L“L(X)”L4D＝30（L－S）2／π／R／H：C＝L－S－（L－S）5／90／（R×H）2：B＝A＋D(-1)K：E＝A＋3D(-1)KL5U“X(Z)＝”X＋CcosB▲L6V“Y(Z)＝”Y＋CsinB▲L7{G}：G“B(L)：”Q“Q”L8F“X(L)＝”U＋Gcos(E＋Q＋180)▲L9I“Y(L)＝”V＋Gsin(E＋Q＋180)▲L10{J}：J“B(R)”L11M“X(R)”＝U＋Jcos(E＋Q)▲L12N“Y(R)＝”V＋jsin(E＋Q)▲L13Goto2卵形曲線坐標計算X＝1,D＝2)”L1 S“Km-YH”：E“X(YH)：”F“Y(YH)：”G“ALF：”B“R1”：D“A：”K“K(L＝1,R＝2)”：Q“R1-R2L2Lb12L3{Z}：Z“L(X)”L4J“L1＝”D2／B：R“RP”＝D2B／(D2＋(-1)Q(Z－S)B)：L“LP＝”D2／RL5M＝(L－J)－(L5－J5)／40／D4＋(L9－J9)／3456／D8L6N＝(L3－J3)／6／D2－(L7－J7)／336／D6＋(L11－J11)／42240／D10L7T＝G－(-1)Q(-1)K×J2×90／D2／πL8X“X(Z)＝”E＋(-1)QMcosT－(-1)KNsinT▲L9Y“Y(Z)＝”F＋(-1)QMsinT＋(-1)KNcosT▲L10A“ALF(P)＝”G＋(-1)K(Z－S)×90×(1／B＋1／R)／πL11{H}：H“B(L)：”U“Q”L12W“X(L)＝”X＋Hcos(A＋U＋180)▲L13V“Y(L)＝”Y＋Hsin(A＋U＋180)▲L14{C}：C“B(R)”L15 I“X(R)＝”X＋Ccos(A＋U)▲L16 P“Y(R)＝”Y＋Csin(A＋U)▲L17Goto2公路逐樁坐標計算4800程序公路逐樁坐標計算程序(可以計算對稱、不對稱緩和曲線)Lb10Z=？V=？W=V+2：Fixm｛K｝Lb11K>Z[W+5Z+4]=>W=W+1:Goto1⊿(判斷樁號在哪個交點范圍，就是該交點曲線起點至下一交點曲線起點)S=K-Z[W+5Z+3]（計算該樁號與曲線起點的距離）R=Z[W+2Z+2]：L=Z[W+3Z+2]:E=Z[W+4Z+2]（讀取該交點曲線要素R、Ls1、Ls2）Pol（Z[W]-Z[W-1],Z[W+Z+2]-Z[W+Z+1]）（計算該交點與下一交點直線方位角）J<0=>J=J+360⊿A=JPol（Z[W-1]-Z[W-2],Z[W+Z+1]-Z[W+Z]）（計算該交點與上一交點直線方位角）J<0=>J=J+360⊿C=A-J:A=J（計算偏角）W=V+2=>Goto2⊿（如果樁號在起點與第一交點曲線起點之間，則轉(zhuǎn)Lb12）I=Abs（tan（c÷2））M=L÷2-L^3÷240R^2:N=E÷2-E^3÷240R^2P=L^2÷6R-L^4÷336R^3-R（1-cos（90L÷πR））Q=E^2÷6R-E^4÷336R^3-R（1-cos（90E÷πR））D=（P-Q）I÷2:F= （P+Q+2R ）I÷2M=F+M-D:Q=F+N+DN=πRAbsC÷180+（L+E）÷2X=Z[W-1]-McosAY=Z[W+Z+1]-MsinAM=Z[W-1]+Qcos（A+C）V=Z[W+Z+1]+Qsin（A+C）Q=AbsC÷CS≤L=>P=0:Goto3⊿（如果樁號在第一緩和曲線內(nèi)，則轉(zhuǎn)Lb13）S≤N-E=>S=S-L:Goto4⊿（如果樁號在圓曲線內(nèi)，則轉(zhuǎn)Lb14）S≤N=>S=N-SQ=-Q:A=A+C-180:X=M:Y=V:L=E：P=180:Goto3⊿（如果樁號在第二緩和曲線內(nèi)，則轉(zhuǎn)Lb13）P=A+C:S=S-N:D=M+ScosP:F=V+SsinPGoto6（如果樁號在直線內(nèi)，則轉(zhuǎn)Lb16）Lb12P=A+CD=Z[W-1]+ScosPF=Z[W+Z+1]+SsinP:Goto6Lb13I=S-S^5÷40R^2÷L^2+S^9÷3456R^4÷L^4J=Q（S^3÷6RL-S^7÷336R^3÷L^3）P=P+A+90QS^2÷πRL:Goto5Lb14M=90（2S+L）÷πRI=RsinM+L÷2-L^3÷240R^2J=Q（L^2÷24R+R（1-cosM））P=A+QMLb15D=X+IcosA- jsinA:F=Y+JcosA+IsinALb16D″X=″◢ （結(jié)果顯示 X坐標）F″Y=″◢ （結(jié)果顯示 Y坐標）P″AT=″◢ （結(jié)果顯示該樁號方位角）{BO} ：B″S″O⊿″″ （輸入邊樁距離，交角）P=P+OL″XB″=D+BcosP◢ （結(jié)果顯示邊樁 X坐標）M″YB″=F+BsinP◢ （結(jié)果顯示邊樁 Y坐標）以上是坐標計算程序，括號內(nèi)是程序計算的大致原理及說明，中間部分為直線、圓曲線、緩和曲線計算的各種公式，大家也知道，書上也有。該程序是最基本的，如再加幾條語句就可以處理斷鏈。該程序還需變量輸入，把整條線路曲線要素一次輸入到 4800 中，就可以計算線路任意中、邊樁坐標。如有需要，需 RMB￥200，附有變量輸入方法及詳細使用說明。聯(lián)系： 曲線任意里程中邊樁坐標正反算 (CASIOfx-4800P 計算器)程序曲線任意里程中邊樁坐標正反算 (CASIOfx-4800P 計算器)程序一、程序功能本程序由一個主程序 (TYQXjs)和兩個子程——正算子程序(SUB1) 、反算子程序(SUB2)序構(gòu)成，可以根據(jù)曲線段——直線、圓曲線、緩和曲線（完整或非完整型）的線元要素（起點坐標、起點里程、起點切線方位角、線元長度、起點曲率半徑、止點曲率半徑）及里程邊距或坐標，對該曲線段范圍內(nèi)任意里程中邊樁坐標進行正反算。另外也可以將本程序中核心算法部分的兩個子程序移植到其它相關(guān)的程序中，用于對曲線任意里程中邊樁坐標進行正反算。本程序也可以在 CASIOfx-4500P 計算器及CASIOfx-4850P 計算器上運行。二、源程序1.主程序(TYQXjs)"1.SZ=>XY "："2.XY=>SZ "：N：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：C=1÷P：D=(P-R) ÷(2HPR)：E=180÷π：N=1=>Goto1 ：≠>Goto 2Δ←┘Lbl1：{SZ} ：SZ：W=Abs(S-O) ：Prog "SUB1"：X"XS"=X◢Y"YS"=Y◢Goto1 ←┘Lbl2：{XY} ：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+W ◢Z"Z"=Z◢Goto22.正算子程序(SUB1)A=0.1739274226 ：B=0.3260725774 ：K=0.0694318442 ：L=0.3300094782 ：F=1-L：M=1-K ：X=U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD))) ：Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))) ：F=G+QEW(C+WD)+90

：X=X+ZcosF

：Y=Y+ZsinF3.反算子程序

(SUB2)T=G-90 ：W=Abs((Y-V)cosT-(X-U)sinT) ：Z=0：Lbl0：Prog "SUB1"：L=T+QEW(C+WD)：Z=(J-Y)cosL-(I-X)sinL ：AbsZ<1E-6=>Goto1 ：≠>W=W+Z ：Goto 0Δ←┘Lbl1：Z=0：Prog"SUB1"：Z=(J-Y) ÷sinF三、使用說明、規(guī)定以道路中線的前進方向（即里程增大的方向）區(qū)分左右；當線元往左偏時，Q=-1；當線元往右偏時， Q=1；當線元為直線時， Q=0。當所求點位于中線時，Z=0；當位于中線左鍘時，Z取負值；當位于中線中線右側(cè)時，Z取正值。(3) 當線元為直線時，其起點、止點的曲率半徑為無窮大，以 10的45次代替。當線元為圓曲線時，無論其起點、止點與什么線元相接，其曲率半徑均等于圓弧的半徑。(5)當線元為完整緩和曲線時，起點與直線相接時，曲率半徑為無窮大，以10的45次代替；與圓曲線相接時，曲率半徑等于圓曲線的半徑。止點與直線相接時，曲率半徑為無窮大，以 10的45次代替；與圓曲線相接時，曲率半徑等于圓曲線的半徑。當線元為非完整緩和曲線時，起點與直線相接時，曲率半徑等于設(shè)計規(guī)定的值；與圓曲線相接時，曲率半徑等于圓曲線的半徑。止點與直線相接時，曲率半徑等于設(shè)計規(guī)定的值；與圓曲線相接時，曲率半徑等于圓曲線的半徑。、輸入與顯示說明輸入部分：SZ=>XYXY=>SZN? 選擇計算方式，輸入 1表示進行由里程、邊距計算坐標 ；輸入2表示由坐標反算里程和邊距。X0？線元起點的X坐標Y0？線元起點的Y坐標S0？線元起點里程F0？線元起點切線方位角LS？線元長度R0？線元起點曲率半徑RN？線元止點曲率半徑？線元左右偏標志(左偏Q=-1，右偏Q=1，直線段Q=0）？正算時所求點的里程Z ？正算時所求點距中線的邊距 (左側(cè)取負，值右側(cè)取正值，在中線上取零 )？反算時所求點的X坐標？反算時所求點的Y坐標顯示部分：XS=×××正算時，計算得出的所求點的X坐標YS=×××正算時，計算得出的所求點的Y坐標S=×××反算時，計算得出的所求點的里程Z=×××反算時，計算得出的所求點的邊距四、算例某匝道的由五段線元（直線+完整緩和曲線+圓曲線+非完整緩和曲線+直線）組成，各段線元的要素（起點里程S0、起點坐標X0Y0、起點切線方位角F0、線元長度LS、起點曲率半徑R0、止點曲率半徑RN、線元左右偏標志Q）如下：S0X0Y0F0LSR0RNQ500.00019942.83728343.5611251631.00269.2561E451E450769.25619787.34028563.3781251631.0037.4921E45221.75-1806.74819766.56628594.5741202554.07112.779221.75221.75-1919.52719736.07228701.893911730.6380.285221.759579.228-1999.81219744.03828781.659804050.00100.0001E451E450、正算（注意：略去計算方式及線元要素輸入，請自行根據(jù)所求點所在的線元輸入線元要素）S=700Z=-5計算得XS=19831.41785YS=28509.72590S=700Z=0計算得XS=19827.33592YS=28506.83837S=700Z=5計算得XS=19823.25398YS=28503.95084S=780Z=-5計算得XS=19785.25749YS=28575.02270S=780Z=0計算得XS=19781.15561YS=28572.16358S=780Z=5計算得XS=19777.05373YS=28569.30446S=870Z=-5計算得XS=19747.53609YS=28654.13091S=870Z=0計算得XS=19742.68648YS=28652.91379S=870Z=5計算得XS=19737.83688YS=28651.69668S=940Z=-5.123 計算得 XS=19741.59118YS=28722.05802S=940Z=0計算得XS=19736.47687YS=28722.35642S=940Z=3.009計算得XS=19733.47298YS=28722.531682、反算X=19831.418Y=28509.726計算得S=699.9999974Z=-5.00018164X=19827.336Y=28506.838計算得S=699.9996493Z=0.000145136X=19823.25398Y=28503.95084計算得S=699.9999985Z=5.000003137X=19785.25749Y=28575.02270計算得S=780.0000035Z=-5.000001663X=19781.15561Y=28572.16358計算得S=780.0000025Z=-0.000002979X=19777.05373Y=28569.30446計算得S=780.0000016Z=4.99999578X=19747.536Y=28654.131計算得S=870.0001137Z=-4.99941049X=19742.686Y=28652.914計算得S=870.0003175Z=-0.00041814X=19737.837Y=28651.697計算得S=870.0002748Z=4.999808656X=19741.5912Y=28722.0580計算得S=939.9999786Z=-5.123024937X=19736.4769Y=28722.3564計算得S=939.9999862Z=-0.000027710X=19733.4730Y=28722.5317計算得S=940.0000238Z=3.00898694[point=10][/point]自己寫的 fx-4800 小程序一、ZBFS(坐標反算)二、NEXY,三、Fixm 。四、F,五、Pol(N-X,E-Y) 。六、J≥0→J=J。七、≠→J≤0→J=J+360 。八、九、說明：“，”表示回車符；“?！北硎窘K止符。十、十一、二、ZBFS-2(坐標反算)十二、NEXY,十三、Pol(N-X,E-Y)。十四、J≥0→J=J。十五、≠→J≤0→J=J+360。十六、S=I。十七、R=J。十八、十九、說明：“，”表示回車符；“。”表示終止符。二十、二十一、三、ZBZS(坐標正算）二十二、二十三、SRXY,二十四、 V=X+Rec(S,R),二十五、 W=Y+J,二十六、 N=V。二十七、 E=W用全站儀測設(shè)公路中樁和邊樁程序程序說明：本程序用于用全站儀測設(shè)中樁和邊樁，采用坐標變換的方法編制而成，在使用 fx-4500 計算器時，可直接輸入以下程序便可使用；而在使用 fx-4800 計算器時，見 fx-4800 程序清單。注意：因程序的 0(零)和O易混淆,所以字母 O用較大的全角符號，零仍用標準字。程序輸入時應連續(xù)輸入，無空格。而且輸入時文件名需和給出的一樣，否則會出現(xiàn)調(diào)用文件錯誤。變量說明： <程序運算符定義 >X輸入擺站點的 x坐標C輸入第一個曲線角樁 x坐標I輸入中間曲線角樁的 x坐標YEF

輸入擺站點的 y坐標D輸入第一個曲線角樁 y坐標J輸入中間曲線角樁的 y坐標最后一個曲線角樁 x坐標R中間曲線半徑 L中間曲線緩和曲線長最后一個曲線角樁 y坐標K第一個K為中間角樁樁號 S中間角樁曲線總長，檢驗用K第二個K為要測設(shè)的樁號 B測中樁時輸入 0，測邊樁時輸入寬度（左輸負值，右輸正值）程序清單注釋文件:F1文件名：1L1L2L3Defm8XYCDIJEFRLK:Fixm:Z[2]=0:ProgBZ[5]=0.5L-LXy3/240R2:Z[3]=Z[5]+(R+L2/24R)AbstanT

┛

2:Z=K-Ans:S=RAbs

Tπ/180+L:S

◢H=Z+Ans:G=RL:Prog2文件:F2

文件名：

2L1Lbl0:Fixm:{K}L2Lbl6:{B}:Norm:K(-1)XyZ[2]>(H-L)(-1)XyZ[2]A=90T/AbsTL3B<0＝>N=0:≠>N=1 ⊿B=AbsBL4

＝

>Goto4

⊿

FracZ[2]

┛

2=0=>Z[6]=Z:

≠>Z[6]

=H

⊿０=Abs(K-Z[6]:K(-1)XyZ[2]>Z[6](-1)XyZ[2]

＝>Goto2

⊿Q=Z[3]+

０:Z[4]=BL5B≠0＝>ProgC

⊿ProgD:P=0L6ProgE:Goto0L7Lbl2:

０≥L＝>Goto3

⊿L8Q=Z[3]-

０+０

Xy5/40G2:Z[4]=

０

Xy3/6G-

０

Xy7/336GXy3:ProgD:P=BL9B≠0＝>ProgF

⊿ProgE:Goto0L10Lbl3:Q=Z[3]-Z[5]:Z[4]=R+L2/24R:ProgDL11Z[8]=U+A(2(０-L)/πR-1+L/πR):B≠0＝>ProgG:≠>P=R⊿ProgE:Goto0L12Lbl4:M=C:０=D:C=E:D=F:E=M:F=０:Z[2]=Z[2]+1:ProgB:Goto6文件:F3文件名：0L1Pol(M,０):W<0＝>W=W+360文件:F4文件名：BL1L2M=E-I:０=F-J:Prog0:T=W:M=I-C:０=J-D:Prog0:U=W:T=T-WT<-180＝>T=T+360⊿T>180＝>T=T-360文件:F5文件名：CL1A=-90(-1)Xy(N+Z[2])文件:F6文件名：DL1Z[7]=I-QcosU+Z[4]cos(U+A):Z[1]=J-QsinU+Z[4]sin(U+A)文件:F7文件名：EL1L2M=Z[7]+PcosZ[8]-X: ０=Z[1]+PsinZ[8]-Y:Prog0W◢Fix3:V ◢Norm文件:F8 文件名：FL1Z[8]=U+A+A ０2/πG+90+A(-1)Xy(N+Z[2])文件:F9 文件名：GL1P=R+AB(-1)Xy(N+Z[2])/90此處L2句加入：M+X◢０+Y◢可顯示測點坐標值CASIO4800 超高及高程（源程序）超高及高程（源程序）L1: R:L“LS”:Z“ZH”:H“HZ”:V“F)÷2:X“X(R+1,L-1)”:P“P(R+,L-1)”

RS”:F

“I1

”:G“I2

”:O“BPH”:N“BP”:T=Vabs(G

－L2:Lbl6: ｛S｝L3:S+T －N<0D=O+(S －N)F:S－N<0D=O+(S－S+N)2÷ 2÷V:D=O+(S －N)G D“H=”◢L4:R<E=:M=Goto1:R<E=:M=Goto1:

－N)F

－U(S

－N+T)2

÷2÷V:S

－N≤TD=O+(S

－N)G

－U(TL5:Lbl1:M ≥LC=L:Goto2:C=ML6:Lbl2:PX>0Goto3:Goto4L7:Lbl4:S≤ZI=0.02:Goto5:S<Z+CI=0.02－(0.02－EPX)÷C×(S－Z):Goto5:S<H－CI=EPX:Goto5:S<HI=EPX+(0.02－EPX)÷C×(S－H+C):Goto5:I=0.02:Goto5L8:Lbl3:Z[1]=0.04(EPX+0÷.02)×C:S≤Z+Z[1]I=0.02:Goto5:S<Z+CI=0.02+(EPX－0.02)÷(C－Z[1])×(S－Z－Z[1]):Goto5:S<H－CI=EPX:Goto5:S<H－Z[1]I=EPX－(EPX－0.02)÷(C－Z[1])×(S－H+C):Goto5:I=0.02:Goto5L9:Lbl5:{B}:W“HB”=D－IB◢L10:Goto6說明：R:平曲線半徑L“LS”:緩和曲線長P：偏轉(zhuǎn)，Z=－1,Y=1X“XL：”線路（Z,Y）左線＝－1,右線＝1Z“ZH”:直緩點里程H“HZ”：緩直點里程V“RS”：豎曲線半徑F“I1”：前坡堵（帶符號）G“I2：”后坡度值（帶符號）O“BPH”：變坡點高程N“BP”：變坡點里程T：切線長E：最大超高值C：漸變段長度（LC）I：橫坡值M：最大漸變長度（Lmax）B：距中心距離W“HB”：邊部高程公路中線坐標計算程序[CASIO4500]公路中線坐標計算程序[casio4500]ZXCLCXDefm16Lbl1X=1:{X}:X“ZJF-1ZBF-2TT-3SJD-4P2L-5B2A-6”X=1＝>Prog"ZJF"△X=2＝>Prog“ZBF”△X=3＝>Prog“TT"△X=4＝>Prog“SJD"△X=5＝>Prog“P2L"△X=6＝>Prog“B2A"△Goto1ZBF“-----ZXZBjs-----”:{YD}:Y“I1”“I2”Y≠D=>Z=0rog“FP”:A=B:B=0:Goto3△X=0:{TPXVWU}:T“HJD(X+Yi)”Lbl1X“XJA=1”“JD(X+Yi) ”:U“JDLC”:Y=Arg(P-T)Y<0=>Y=Y+360 △Y“I◢”C“L”=Abs(P-T)◢X=0=>V=P:Goto2 △V“JDB(X+Yi) ”=Arg(V-P)D<0=>D=D+360 △D“I◢”Z“L”=Abs(V-P)◢Prog"FP":A=BLbl2W“QJD(X+Yi) ”=Arg(W-V)D<0=>D=D+360 △D“I◢”C“L”=Abs(W-V)◢Lbl4Prog

“FP”X=1=>B=B-A

△X=0=>A=B:B=0:Z=0:E=RepP:F=ImpP:Goto3

△M=Rep(P-T):N=Imp(P-T):I=Rep(W-V):J=Imp(W--I÷J) ◢ E“X”=M(F-ImpT) ÷N+RepTU=U+C:C“L2”=√((E-RepV)2+(F-ImpV)2) ◢Lbl3X=0:{X}:X “ZBF(0)ZJF(