4800--4850编程

合集下载

4800,4850,5800计算程序

CASIO4800-4850计算程序1Lbl 02{LUZIPTKGSH}:3L≤626.303=>Q=0:N=3601886.914:E=547501.524:4J=156-27-10: Goto1:5≠>L≤706.303=>Q=626.303:N=3601312.75:E=547751.706:6J=156-27-20:M=-1:C=80:R=4000:Goto3:7≠>L≤768.044=>Q=706.303:N=3601239.517:E=547783.907:8J=155-52-57.35:R=4000:M=-1:Goto2: 9≠>L＜848.044=>Q=848.044:N=3601111.086:E=547843.863:10J=334-25-30.96: R=4000: M=1:C=80:G=180: Goto 3：11≠>L≤1387.746=>Q=848.044:N=3601111.086:E=547843.863:12J=154-25-30.96: Goto1:13≠>L≤1467.746=>Q=1387.746:N=3600624.263:E=548076.846:14J=154-25-30.96:M=1:C=80:R=4000:Goto3:15≠>L≤1527.249=>Q=1467.746:N=3600551.986:E=548111.14:16J=154-59-53.6:R=4000:M=1:Goto2: 17≠>L＜1607.249=>Q=1607.249:N=3600424.659:E=548168.129:18J=336-25-24.59: R=4000: M=-1:C=80:G=180: Goto 3：19≠>L≤2180.31=> Q=1607.249:N=3601424.659:E=548168.129:20J=156-25-24.59: Goto1:21≠>L≤2300.31=>Q=2180.31:N=3599899.434:E=548397.338:22J=156-25-24.59:M=1:C=120:R=2500:Goto3:23≠>L≤2517.342=>Q=2300.31:N=3599899.434:E=548397.338:24J=156-25-24.59:R=2500:M=1:Goto2:25≠>L＜2637.342=>Q=2637.342:N=3599469.659:E=548551.339:26J=344-08-51.73: R=2500: M=-1:C=120:G=180: Goto 3：27≠>L≤3612.157=>Q=2637.342:N=3599469.659:E=548551.339:28J=164-08-51.73: Goto1:29≠>L≤3712.157=>Q=3612.157:N=3598531.917:E=548817.618:30J=164-08-51.73:M=-1:C=100:R=3000:Goto3:31≠>L≤3984.738=>Q=3712.157:N=3598435.875:E=548845.468:32J=163-11-34:R=3000:M=-1:Goto2:33≠> L ＜4084.738=>Q=4084.738:N=3598086.588:E=548974.532:34J=337-01-95:R=3000:M=1:C=100:G=180:Goto3: 35Lbl 136D=L-Q37P=J38X=N+D*COS J▲39Y=E+D*SIN J▲40U“1,0”:U=0=>Goto0 :≠>Goto 4:41Lbl 242D=L-Q43I=90*D/(兀R)44P=J+2IM:P>360=>P=P-360⊿P≤360=>P=P⊿P<0=>P=P+360⊿45X=N+2R SIN I*COM (J+IM)▲46Y=E+2R SIN I*SIN (J+IM)▲47U“1,0”:U=0=>Goto 0: ≠> Goto 4:48Lbl 349D=Abs (L-Q)50B=D︿3/（6RC）-D︿7/(336(RC) ︿3)+D︿11/(42240(RC) ︿5) 51A=D-D︿5/(40(RC) ︿2)+D︿9/(3456(RC) ︿4)52I=tan-1 (B/A)53P=J+3IM+G: P>360=>P=P-360⊿P≤360=>P=P⊿P<0=>P=P+360⊿54X=N+B/SIN I*COM (J+IM)▲55Y=E+B/SIN I*SIN (J+IM)▲56U“1,0”:U=0=>Goto 0: ≠> Goto 4:57Lbl 458W=X+(Z“JL”)*COM (P+(T“JD”))▲59V=Y+Z*SIN (P+T)▲60S“F,B(0,1)”61S=0=>Goto 0:≠>Goto 5:62Lbl 563O=W+K“JULI” COS (P+T+H“JIAODU”)▲64Q=V+K SIN (P+T+H)▲65Goto 0:66说明：67这个程序是根据你自己管段的曲线要素来编的。

CASIO_fx4800P与fx4850P说明书

次并显示计算结果，按键删除显示后，计算结果不被保存；而 fx-4800P 可以完整地显示用户输
入的表达式，按键显示的计算结果被保存在 Ans 存储器中，按键删除显示后，可以按键调出最近一次表达式的计算结果。还可以在完成表达式的计算后，立即按重演键或调出
并显示最近执行的表达式，再应用或键移动光标、输入新字符覆盖当前光标处的字符、应
第 1 章 CASIO fx-4800P 与 fx-4850P 的操作方法与编程技巧
§1.1 fx-4800P 的特点
众所周知，CASIOfx-4500PA 可编程函数计算器早在 80 年代末就已经在我国工程界得到了广泛的应用， fx-4800P 是比 fx-4500PA 更高一个级别的产品(见图 1-1)，它采用硬翻盖设计，除了具有 fx-4500PA 的全部功能以外，还主要改进了下列特性：
30，但π 2 不可以输入为 2，4 3 显然也不可以输入为 43。
③ 计算器遵循括号内的表达式、A、B 型函数，×、÷、＋、－的运算规则，而在 A、B 型函
数内的表达式，省略了×号的乘法运算级别最高。对于有同级别的运算，遵循由左至右的运算规则。
如输入表达式 2π sin 30 最简单的按键方式为 2
。
③ 2500 加上其 15%(答案为 2875)。按键：2500 15
。
④ 3500 减去其 25%(答案为 2625)。按键：3500 25
。
⑤ 168、98 及 734 的和减去其 20%(答案为 800)。
按键操作：168 98 734
20
。
⑥ 原重量为 543 克，将其重量加上 345 克后，总重量为原重量的百分之几(答案为 163.54%)?

4850 4500 4800正反算程序

1．中边桩计算程序<1>主程序名—XYLbl 1:{KM}：K：M：Prog”3”:X=C+Mcos(W-90:Y=D+Msin(W-90:”X=”:X:Pause 0:”Y=”:Y▲Goto1<2>子程序—3PrOg”QX”：Prog”20”<3>曲线要素—QX数据实例：某高速公路中线设计如下图所示K<2000=>B:E:F:O:P:Q:R: K≤4570=>B=25:E=1000:F=1000:O=100P=66:Q=4000:R=800: K≤5040=>B=-17:E=1086.4378:F=1555.1451:O=90:P=91:Q=4570:R=800: K≤5670=>B=35:E=1168.2556:F=2016.0909:O=120:P=74:Q=5040:R=1000: B=-25:E=1180.9927:F=2636.2334:O=150:P=109:Q=5670:R=1500注：无缓和曲线时，O=0.00001当输入里程小于2000时，为新曲线要素输入，可修改值域范围<4> 计算过程—20Z=O：Prog”8”：A=B/AbsB:L=ABπR/180+O:T=(D+RcosN)tan(AB/2)+C-RsinN：K≤Q-T=>C=(K-Q)cosP+E：D=(K-Q)sinP+F:W=P： K≤Q-T+O=>Z=K-Q+T：Prog”9”：K≤Q-T+L-O=>Z=O: Prog”9”:V=90(K-Q+T-O)/ πR:S=2R sinV:N=W+AV:C= C+S cosN:D=D+S sinN:W=W+2AV: Z=Q-T+L-K: Prog”8”:V=P+B-180-AV:W=P+B-AN :C=S cosV+E+T cos(P+B:D=S sinV+F+T sin(P+B: Z=K-Q+2T-L:W=P+B:C=E+ZcosW:D=F+Z sinW△△△△N=Z[3]<5> 缓和曲线方程式—8V=RO:C=Z-Z^5/40V2+Z^9/3456V^4:D=Z^3/6V-Z^7/336V^3+Z^11/42240V^5:N= 90Z^2/πV:V=tan-1(D/C:S=(D^2+C^2)^0.5<6> 子程序—9Prog”8”:C=Scos(P+AV)+E-TcosP:D=Ssin(P+AV)+F-TsinP:W=P+AN2.桩号、边桩距返算程序—KM方式一：试算法（速度慢，程序简单）Fixm:K=Q-T+L-O:Prog”20”Lbl 1:{XY}Lbl 2:POL(X-C,Y-D:K=Isin(J-W+90)+K:Abs sin(J-W+90)<0.0001＝>”K=”:K:Pause0:”M=”:M=Icos(J-W+90▲≠>Prog”20”:Goto2△方式二：解算法（程序复杂，计算速度快）Defm5: Lbl 1:Fixm:{XY}:Pol(X-E,Y-F:Z[1]=I:Z[2]=J:cos(J-P<0＝>K=Q+Icos(J-P: ≠>K=Q+Icos(J-P-B)+2T+L△K≤Q-T＝>M=Icos(J-P+90:Goto5≠>Goto2△△Lbl 2:Z=O:Porg”8”:U=(D+RcosN)+cos(B/2):Z=P+B/2+90A:C=UcosZ+E:D=UsinZ+F:Pol(X-C,Y-D:Asin(J-P-AN+90 A<0＝>K=K-AZ[1]cos(Z[2]-P+90)tan(90 (K-Q+T)^2/πRO):Goto3:≠>Asin(J-P-B+AN+90A>0 ＝>K=K-AZ[1]cos(Z[2]-P-B-90 )tan(90(K-Q+T-L)^2/πRO):Goto4:≠>K=Q-T+O+Asin-1(sin(J-P-AN+90 A))×πR/180 :M=A(I-R:Goto5△△Lbl 3:Z=K-Q+T:Prog”9”Pol(X-C,Y-D:K=K+Isin(J-W+90 :Abs sin(J-W+90<.001＝>M=Icos(J-W+90 :Goto5: ≠>Goto3△Lbl 4:Z=Q-T+L-K:Prog”8”:V=P+B-180-AV:W=P+B-AN:C=cosV+E+Tcos(P+B:Pol(X-C,Y-D:K=K+Isin(J-W+90 :Abs sin(J-W+90 <.001＝>M=Icos(J-W+90:Goto5: ≠>Goto4△Lbl 5:”K=”:K:Pause0:”M=”:M▲注：程序“20”为中桩计算程序Q—直缓点里程；T—切线长，L—曲线长；O—缓和曲线长；W—切线方位角Pol—极坐标转换，为涵数特别注意数字0和字母O的区别一、程序功能本程序由一个主程序(TYQXjs)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

测量程序4800和4850

1.主程序(ZCX1)"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=>Goto 1：≠>Goto 2Δ换行Lbl 1：{SZ}：SZ：W=Abs(S-O)：Prog "SUB1"：X"XS"=X◢Y"YS"=Y◢F"FS"=F-90◢Prog“ZDS”2.Lbl 2：{XY}：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+W◢Z"Z"=Z◢Prog“ZDS”2. 正算子程序(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+ZsinF换行3. 反算子程序(SUB2)T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog "SUB1"：L=T+QEW(C+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ换行Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF换行Lbl 0：S≤终点桩号=>O=起点桩号：U=起点X：V=起点Y：G=起点在曲线上的方位角：H=曲线或直线长：P=起点半径（直线段和缓和曲线为无穷大）：R=终点半径（直线段和缓和曲线为无穷大）：Q=路线转向（左“—”右“+”直线为“0”）：Prog “ZCX1”：Goto 0⊿↙Lbl 1：S≤终点桩号=>O=起点桩号：U=起点X：V=起点Y：G=起点在曲线上的方位角：H=曲线或直线长：P=起点半径（直线段和缓和曲线为无穷大）：R=终点半径（直线段和缓和曲线为无穷大）：Q=路线转向（左“—”右“+”直线为“0”）：Prog “ZCX1”：Goto 0⊿↙。

卡西欧4800和4850计算程序

卡西欧4800\4850程序一：坐标反算：(ZBFS)①程序步骤：C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C,F-D:I”S1-2=”▲J≤Ｏ＝＞J=J+360△J”A1-2=”②操作过程：ZBFS→EXE→输入X1值→EXE→输入Y1值→EXE→输入X2值→EXE→输入Y2值→EXE→EXE→EXE注：①：S1-2:计算得出的距离。

②：A1-2:计算得出的角度。

（按shift°′″转换为60进制的角度）③：此程序可循环计算。

③案例：已知：D1坐标（7811.23，606.136），D2坐标（7805.915，431.910）求解： a D1-D2(A1-2)=268°15′09.56 ″S D1-D2(S1-2)=174.3071二：坐标正算：(ZBZS)①程序步骤：C”X1”:D”Y1”:S”S1-2”:R”A1-2”:Fixm:X”X2”=C+Rec(S,R) ▲Y”Y2”=D+J②操作过程：ZBZS→EXE→输入X1值→EXE→输入Y1值→EXE→输入S1-2距离值→EXE→输入A1-2角度（例268°15′09.56 ″）值→EXE→EXE注：①：X2计算得出坐标，Y2计算得出坐标。

②：此程序可循环计算。

③案例：已知：D1坐标（7811.23，606.136），两点距离：S D1-D2(S1-2)=174.3071方位角：a D1-D2(A1-2)=268°15′09.56 ″求解：D2坐标（7805.915，431.910三：线路曲线计算程序:（XLQXJS）(万能公式)①程序步骤：1.Lbl 0:{E G}：A“XA”：B“YA”：C“CA”：D“1/RA”：E“1/RB”：F“DKA”：G“DKA”2.Lbl 1:{H O R}:H“DKI”：O“DL”：R“DR”：H＞G=﹥Goto 3⊿3.P=(E-D)/Abs(G-F):Q=Abs(H-F):I=P×Q:T=D+I4.J=C+(I+2D)×Q×90/π ▲5.M=C+(I/4+2D)Q×45/(2π):N=C+(3I/4+2D)Q×135/(2π):6.K=C+(I/2+2D)Q×45/π7.X=A+Q(cosC+4(cosM+cosN)+2cosK+cosJ)/12 ▲8.Y=B+Q(sinC+4(sinM+sinN)+2sinK+sinJ)/12 ▲9.U“XL”=X+Ocos(J-90) ▲10.V“YL”=Y+Osin(J-90) ▲11.W“XR”=X+Rcos(J+90) ▲12.Z“YR”=Y+ Rsin(J+90) ▲13. A=X：B=Y：C=J：D=T：F=H：Goto 1Lbl 3: A=X：B=Y：D=E：F=G：C=J：Goto 0程序说明：A- 曲线元起点A的坐标；B- 曲线元起点B的坐标；C- 曲线元起点A的切线坐标方位角；F- 曲线元起点A的里程；G- 曲线元起点B的里程；H- 曲线上待求点i的里程；D- 曲线元起点A的曲率；E- 曲线元终点B的曲率；XL-左边线点位X坐标；YL-左边线点位Y坐标；XR-右边线点位X坐标；YR-右边线点位Y坐标；X- 中线点位纵坐标；Y- 中线点位横坐标；DL-左边线距中线平距；DR-右边线距中线平距；该程序需要输入的数据为：（1）.曲线元起点A的坐标及切线坐标方位角，计算器上用“XA”，“YA”，“CA”显示；（2）.曲线元起点A和B的曲率，计算器上用I÷RA，I÷RB显示（曲线左偏时取“-”）；（3）.曲线元起点A和终点B的里程，计算器上用“DKA”，“DKB”显示；（4）.输入待求点里程和该点距左右的水平距离，计算器上用“DKI”，“DL”，“DR”显示；每算完一个待求点的中线及边线坐标，程序又让输入下一点的“DKI”，“DL”，“DR”当输入的“DKI”大于“DKB”时程序中显示“1÷RB”和“DKB”，此时输入下一个曲线元起点的曲率和里程，然后重复步骤（4），即可计算下一个曲线中线及边线点位坐标。

对卡西欧4800计算器编程

在计算能力飞速发展的今天，各种计算工具应运而生代替了以往我们大量、繁琐的分析计算工作，简化了计算步骤节省了计算时间。

在种类繁多的计算工具中适合测量工作的计算器也大有所在。

比如有一款SHAPPC-E500计算器，它相当于是基于BASIC语言的微型电脑，功能也相当齐全能编写较复杂的程序、内存也相当大有32KB、64 KB、128 KB、256 KB、几种型号，同时具有输入输出通讯设备。

但价格不菲，现在市场价格大约为2000元到3000元，个人经济承受能力有限。

通过几年的实际测量工作我个人为认为真正适合测量人员较为专业的计算器应首推基于科学计算器具有简单程序语言的casio４x ００计算器。

它有很多优点：一、使用简单，最大的一个优点就是只要具有初中水平会运用数学公式对函数有一般的了解的人，能基本理解casio４x００的内装函数即可进行简单的编程。

非常适合测量初学者和数学程序爱好者学习，也适合专业的测量人员的使用和能力提高。

不像其它的专业一定要具有较高的专业技术水平才能进行电脑编程。

二、成本低廉，价格在400元左右。

三、携带方便，体积很小可随时放在口袋里随拿随用。

casio系列较好的编程型号有casio４５００（以下简称４５００）、casio４８００（以下简称４８００），（好象近来还推出了一款casio4850）前者较内存小，只有１１０３个字节，能应付一些较为简单的公式计算和科学计算，但由于内存有限，对一些较复杂或子程序过多的程序就力不从心了，不能出色的完成测量任务。

４８００就比４５００有较大的改进，４８００内存达到的４５００个字节，而且显示屏是４５００的几倍大，能更准确的显示数据，内装函数字符一目了然。

且具有简单的人机对话功能，出现了菜单子菜单。

４８００还在４５００增加了啊佛加德罗常数、万有引力长常数、详见《操作说明书》。

４８００有很多很实用的功能如他有一个公式解答功能，其原理是用牛顿法解方程。

举个例子：有一方程式：a=2b-c求当a=2、c=5时的b值。

卡西欧fx-5800P程序编制与卡西欧fx-4800P或4850P程序编制的异同点

Ifend表达的是本条件判断结束。如果条件判断成立，本条件判断结束后不执行下一条件判断。如果条件判断不成立，本条件判断结束后执行下一条件判断。最后一条件判断中，可以使用“If <条件>：Then < 命令1> ： <命令2> ： <命令n> :Ifend：Else <命令1> ： <命令2> ： <命令n> ：Ifend”结构结束。 b、For语句的使用 For <表达式，始值>→ <变量> To <表达式，终值>：Next 控制变量中可加步长（从始值开始）
b、扩充变量组 5800：在程序编制前使用“shift” DimZ命令具体操作：n → DimZ 4800：具体操作：Defm n c、中间计算值赋值 5800操作为：X-M →U Y-N →V 4800操作为： U= X-M V =Y-N d、赋值显示 5800操作为：X-M →U：“U”：U▲ 4800操作为： U= X-M ▲ e、度分秒的显示问题 ”FWJ=”:J ▲DMS ▲ 如果不使用▲DMS ，则需要按 “shift” ，再按“°′″”来显示 4、括号使用
卡西欧fx-5800P程序编制与卡西欧fx-4800P 或4850P程序编制的异同点
1、5800与4800程序结构编制基本一致：指的是运行结构和运行顺序。 2、二者调用子程序：a、5800主程序调用不许使用 RETURN命令；b、子程序编制完后，要用RETURN命令结束，与4800不同点之一。 3、变量赋值： a、单变量赋值，有两点须注意：基本语句：“A”？A或“A”？→A 两者区别： “A”？A表达了提示符显示该变量的当前值，不输入任何内容时，按“EXE”执行的是当前显示值继续。而“A”？→A表达了不显示当前值，必须输入一个值才能继续。 4800的单变量赋10 Step 0.5：A×A →B ： B ▲ :Next

4850chengxu

一、程序简介本套程序是共有2个主程序，7个子程序。

包括了路线坐标正反算、竖曲线、超高加宽、边坡放样及断面计算等程序。

适用于CASIO4800/4850，4850用户需要修改清单内结果显示的字段（例如：X“X=”◢改为“X =”：X◢或G“Z”=H+G◢改为“Z=”：G=H+G◢）。

二、程序功能本着用最懒的方法做最繁杂的事（当然认真更重要！），根据不同需要，本程序共设置了6种计算模式J-JSMS：1-ZS，正算模式：适用于中、边桩计算，里程桩号及边距=>平面坐标。

2-FI，放样1模式：适用于路面施工放样，里程桩号及边距=>平面坐标、路面高程。

3-FS，反算模式：适用于任意点里程、边距计算，任意点坐标=>里程、边距。

4-F2，放样2模式：适用于路基施工放样，任意点坐标=>里程、边距、填挖高度，边坡超欠挖宽度，最近平台高度（<1.5）。

5-Z0，设计高程计算，即竖曲线计算，里程桩号=>设计高程。

6-DM，断面计算，里程桩号、边距=>设计高程、加宽值、指定边桩的高程。

三、计算范围平曲线：直线-缓和曲线（或无）-圆曲线-缓和曲线（或无）-直线中任意桩号，允许两缓和曲线不等长，即可以算卵形曲线，回头曲线竖曲线：前交点SYZ到后交点SZY中的任意桩号四、数据输入及各计算模式运行过程1、正算模式：1-ZS，运行主程序J-PQX（平曲线），计算过程如下：显示输入计算结果备注DJ？交点桩号JDX？交点X值JDY？交点Y值FJ？计算方位角ZJ：Z-，Y+？交点转角左转为负值，右转正值R？圆曲线半径LS1？第一缓和曲线LS2？第二缓和曲线T1= 前切线T2= 后切线L= 曲线长度E= 外距ZH= ZH桩号若不计算平曲线要素桩号，程序清单J-PQX内有下划线部分可以不写入计算器内。

HY= HY桩号QZ= QZ桩号YH= YH桩号HZ= HZ桩号1-ZS,2-F1,3-FS4-F2,5-Z0,6-DMMS？计算模式选择，此模式请输入：1 总共有6种计算模式供选择，符号说明见“二、程序功能”。

卡西欧4800和4850放样坐标计算程序

已知任意点的里程、半宽，计算该点平面坐标的程序前言1、本计算程序适用范围仅限于某ZH点至下一相邻ZH点间的线路范围。

2、使用本程序前需进行导线点坐标计算并求出各个曲线的曲线要素如第一切线方位角、转向角、曲线半径、缓和曲线长度、切线长度、曲线长度等。

3、使用过程中要注意长短链对里程传递的影响关系，否则将导致计算结果错误。

4、往计算器中输入程序符号时将⊿改为。

5、本程序仅供参考,实际使用时请先进行验证。

6、本程序适用的计算器型号为卡西欧4800/4850。

“wangDongping” (程序名称)Defm 13↙XYVWRLZTC：Fixm↙Z[4]=0.5×L－L3÷(240×R2)↙Z[2]=Z[4]+(R+L2÷24÷R)Abs tan(T÷2)：G=Z＋Ans：Fix3：G＂JDDK=＂◢Z[9]=R Abs Tπ÷180＋L：Z[9]＂QXL=＂◢Z[13]=Z＋Ans：Z[13]＂HZDK=＂◢D=C＋T：D＜0⇒D=D＋3600⊿D＞3600⇒D=D－3600⊿D◢Z[11]＂HZX＂=V＋Z[2]×(CosC＋CosD)◢Z[12]＂HZY＂=W＋Z[2]×(SinC＋SinD)◢Norm：Prog＂WEIXONG2＂(程序“WaangDongping”运行到此时，自动转到程序“WangDangping2”进行下一步运行)“WangDongping2”（程序名称）Lb1 0：{K,S,Q,E,F,N,U,H}：P=Q＋Abs(E－F)×N：K≤Z⇒Goto1⊿Z＜K⇒K≤Z＋L⇒Goto2⊿Z＋L＜K⇒K≤Z＋Z[9]÷2⇒Goto3⊿Z＋Z[9]÷2＜K⇒K＜Z[13]－L⇒Goto4⊿Z[13]－L≤K⇒K＜Z[13]⇒Goto5⊿ K≥Z[13]⇒Goto6⊿Goto0↙Lb1 1：Z[1]=V＋(Z－K)×cos(C＋1800)＋Pcos(C＋S)＋Ucos(C＋H)：Z[3]=W＋(Z －K)×sin(C＋1800)＋Psin(C＋S)+Usin(C+H)：Goto7⊿Goto0↙Lb1 2：Z[5]=K－Z－(K－Z)5÷(40R2L2)：Z[6]=（K－Z）3÷（6RL）：A=T÷AbsT：Z[7]=(K－Z)2×180÷(6πRL)：Z[1]=V＋√((Z[5])2＋(Z[6])2)×cos(C＋A×Z[7])＋Pcos(C＋3×A×Z[7]＋S)＋Ucos(C＋3×A×Z[7]＋H)：Z[3]=W＋√((Z[5])2＋(Z[6])2)×sin(C＋A×Z[7])＋Psin(C＋3×A×Z[7]＋S)＋Usin(C＋3×A×Z[7]＋H)：Goto7⊿Goto0↙Lb1 3：B=L×180÷(2πR)：Z[8]=B＋(K－Z－L)×180÷(πR)：Z[5]=Z[4]+R×sinZ[8]：Z[6]=R+L2÷(24R)－RcosZ[8]：Z[7]=tan-1(Z[6]÷Z[5])：Z[10]=(2×(K－Z)－L)×180÷(2πR)：A=T÷AbsT：Z[1]=V+ √((Z[5])2+(Z[6])2)×cos(C+A×Z[7])+Pcos(C+A×Z[10]+S)＋Ucos(C＋A×Z[10]＋H)：Z[3]=W+√((Z[5])2＋(Z[6])2)×sin(C+A×Z[7])+Psin(C+A×Z[10]+S)＋Usin(D＋A×Z[10]＋H)：Goto7⊿Goto0↙Lbl 4：B=L×180÷(2πR)：Z[8]=B+(Z[13]－K－L)×180(πR)：Z[5]=Z[4]+R×sinZ[8]：Z[6]=R+L2÷(24R)－RcosZ[8]：Z[7]=tan-1(Z[6]÷Z[5])：Z[10]=(2×(Z[13]－K)－L)×180÷(2πR)：A=－T÷AbsT：Z[1]=Z[11]+ √((Z[5])2+(Z[6])2)×cos(D+180+A×Z[7])+Pcos(D+A×Z[10]+S)＋Ucos(D＋A×Z[10]＋H)：Z[3]=Z[12]+√((Z[5])2+(Z[6])2)×sin(D+180+A×Z[7])+Psin(D+A ×Z[10]+S)＋Usin(D＋A×Z[10]＋H)：Goto7⊿Goto0↙Lbl 5：Z[5]=Z[13]－K－(Z[13]－K) 5÷(40R2L2)：Z[6]=(Z[13]－K) 3÷(6RL)：A=－T÷AbsT：Z[7]=(Z[13]－K)2×180÷(6πRL)：Z[1]=Z[11]+ √((Z[5])2+(Z[6])2)×cos(D+180+A×Z[7])+Pcos(D+3×A×Z[7]+S)+Ucos（D+3×A ×Z[7]+H）：Z[3]=Z[12]+√((Z[5])2+(Z[6])2)×sin(D+180+A×Z[7])+Psin(D+3×A×Z[7]+S)＋Usin(D+3×A×Z[7]+H)：Goto7⊿Goto0↙Lb1 6：Z[1]=Z[11]+(K－Z[13])×CosD+Pcos(D+S)＋Ucos(D+H)：Z[3]=Z[12]+(K －Z[13])×sinD+Psin(D+S)＋Usin(D+H)：Goto7⊿Goto0↙Lb1 7：M=Z[1]－X：O=Z[3]－Y：Pol(M,O)：J＜0 J=J+3600⊿Fix3：J◢I◢P◢E－F◢M+X◢O+Y◢Norm：Goto0↙坐标程序符号说明1、X，Y：测站坐标2、V，W：ZH点X,Y坐标3、R：曲线半径4、L：缓和曲线长度5、Z：ZH点里程6、T：曲线转向角（左负右正）7、C：第一切线方位角8、K：放样点里程9、Q：放样半宽（平行于S角非切线边在里程点左右侧的宽度，若已知实际半宽P，则Q=P；若放路基边桩则输入变坡点处的半宽P）10、S：放样点与其对应的中线点连线和该中线点切线的夹角（左负右正，即放样点在线路左侧S为负，同理在右侧为正）11、F：放样点地面标高12、N：路基边坡坡度13、E：变坡点设计标高14、P=Q+Abs（E－F）×N（Abs为绝对值），若P为已知，则Q的值直接输入P的值即可，E、F、N均为零。

卡西欧计算器4850使用说明

举个例子：有一方程式：a=2b-c求当a=2、c=5时的b值。

将该方程式存入公式存贮器中，：先按按解答键“SOLVE” a输入2 、c输入5，再按解答键“SOLVE”计算器就会显示：b=3.5。

此项功能被称为自动解答功能。

同时它也是非常实用的，在实际工作中通常要有一个经常使用的小公式，可借助它来完成。

4800使用的程序语言可以算做简单的BASIC语言，4850基本上跟4800一样，有的命令如GOTO（转移到）、PAUSE（暂停）就与BASIC语言的一模一样。

现在4800的程序语言来说说。

其主要命令有：1、=>…… 条件转移成立符号，其用法相当于BASIC中的IF……THEN（假设语句相当于假如……然后，IF相当于条件……THEN相当于结果）语句2、≠>……条件转移不成立符号，其用法相当于BASIC中的IF……ELSE语句通常二者连用，相当于BASIC中的IF……THEN……ELSE语句（它的英语形式一般为if a＞b then c＞d else if b＞a the……）3、⊿…… 条件转移结束符号，与=>和≠>配合使用，放在条件语句最后面。

4、LbI……标记命令。

用于将一段语句作转换标记。

后可接字母、数字、符号，但不能超过两个字节，如不能用≥10的数字作行标，否则会出现出错信息。

5、Goto…… （条件）转移命令。

前面可加条件语句，与BASIC中的GOTO作用相同。

通常与LbI一起用，如果所转移的行号无效，则会显示：GO ERROR（详见说明书）出错信息6、Dsz……减量循环命令。

可减少未知数的数量。

7、Isz……增量循环命令。

8、Pause……暂停命令。

后可接0~9之间的整数n，可使某一数据显示n/2秒钟，然后继续运行下面的程序。

它被认为是一个语句。

9、Fixm……变量锁定命令。

该命令能使其所有变量值（A~Z）均当成定数处理。

当程序运行时，将不需要输入变数（“{}”内的变数除外），而是将存贮器中原有的数值来完成计算。

4800和4850及5800编程区别

FX-4500PA,FX-4800P,FX-4850P,FX-5800P计算器区别在FX-4500PA,FX-4800P,FX-4850P,FX-5800P等计算器的背面都有复原孔（RESET或P）,每次更换新电池或新计算器拔掉隔电膜后,都必须先点击该小孔进行初始化处理。

CASIO FX-4850 4800工程计算器在2008年3月停产现在的升级产品FX-5800PFx-5800价格和fx-4850相当，但程序不能共享，Fx-5800的编程语言更接近BASIC程序语言。

2、5800变量输入语句改为?，赋值语句改为—＞。

程序中和1850一样，除可以使用26个字母变量A~Z存储数据和独立于字母变量存储器的额外变量存储器Z[1]，Z[2]，Z[3]，Z[4]……外，另外可以使用统计串列存储器List X，List Y，List Freq。

3、程序命令由4850的16个增加到了34个，取消了⊿、≠=>、｛｝命令，使得由⊿、≠=>、｛｝等重要程序命令构成的fx-4850程序无法在5800上运行。

同时，5800的条件语句修改为更易于使用的If~Then~Else~IfEnd形式，新增For~To~Step~Next循环、条件循环Do~LpWhile与While~WhileEnd语句。

4、5800增加了数据通讯功能，使用SB-62数据线连接好两台机器的3Pin通讯口可方便地传输程序与数据，避免了4850只能以手工输入的方式录入程序，容易出现程序输入错误的弊端。

5、5800有矩阵计算功能，可以编写严密平差程序时组成与解算法方程，能够满足测量专业测量平差编程的计算要求。

同时字母变量除了可以输入26个大写英文字母外，还能用小写英文字母、希腊字母。

现在的升级产品FX-5800PFx-5800价格和fx-4850相当，但程序不能共享，Fx-5800的编程语言更接近BASIC程序语言。

2、5800变量输入语句改为?，赋值语句改为—＞。

Casio fx4800任意曲线程序

Casio fx-4800(4850)计算器程序任意曲线坐标计算程序(HHQX)LbI 0 ：｛QD｝←┘G=I=>Goto 1ΔG=-1=>Goto 2△G=0=>Goto 4 △←┘LbI 1 ：L=C÷R ：M=L+Q-K ：F=A+90(M2-L2)÷（πC）：F◢T=A-90L2/πC：Goto 3△←┘LbI 2 ：L=C÷R ：M=L+K-Q ：:T=A+90L2÷πC+180 ：F=A-90(M2-L2)÷πC:F◢Goto 3△←┘LbI 3 ：N=(M-L)-(M^5-L^5)÷40C2+(M^9-L^9)÷3456C^4-(M^13-L^13)÷599040C^6：E=(M^3-L^3)/6C-(M^7-L^7)÷336C^3+(M^11-L^11)÷42240C^5-(M^15-L^15)÷9676800 C^7 ：X=I+NcosT-GEsinT+Dcos(F-90)◢Y=J+NsinT+GEcosT+Dsin(F-90)◢Goto 0←┘LbI 4 ：L=Q-K ：N=Rsin(180L÷πR)：E=R(1-cos(180L÷πR))：F=A+180L÷πR：F◢X=I+NcosA-EsinA+Dcos(F-90)◢Y=J+NsinA+EcosA+Dsin(F-90)◢Goto 0说明:当G=1时,计算第一缓和曲线或由大圆到小圆的卵型曲线；G=-1时,计算第二缓和曲线或由小圆到大圆的卵型曲线；G=0时,计算直线或圆曲线。

R为起算点半径,线路左转时R、C取负值，右转时R、C取正值，计算第一缓和曲线时R取10^45。

计算第二缓和曲线时R取圆曲线半径，计算直线时R取10^45.此计算思路是沿路线的前进方向依次进行分段计算。

4800,4850编写程序详解

14
综合上述关于直线，圆曲线，缓和曲线缓和段程序的介绍，道路测量就是这三种线型的组合，计算方法就是这些。希望做的同事在回去工作中多看关于这方面的书籍，了解其原理，能够自己写出更多更好实用性更强的程序。
2020/8/27
15
道路测量直线，缓和曲线。圆曲线综合程序
路线程序（LX） D“JD”：A：﹛R﹜：C“LS” P=C2÷24÷R-C4÷2688÷R3+C6÷506880÷R5-
V=F+ZA：圆曲线上任意点切线方位角
2020/8/27
9
第4步： X=M+Hcos(F+ZI) ▲ 所求任意点X坐标 Y=N+Hsin(F+ZI) ▲ 所求任意点Y坐标第5步： J”J=” 距中桩距离。左-，右+。
X=X+Jcos(V+90) ▲
Y=Y+Jsin(V+90) ▲ 边桩X,Y坐标结果。 Goto 0 同上
经多年的实践经验，自己的总结，编程出了一个通过偏角法进行计算坐标的程序，适用于各种圆曲线，包括匝道的卵形曲线中的圆曲线，希望大家在以后的工作中能够用上。
2020/8/27
7
第1步： M “X0=”:N”Y0=”:R”R=“: F”F=”: 0”0=”:Z”Z=” 同直线赋予值一样，只介绍 R—半径，Z—转角，左-右+
5
第4步
X=X+J*cos(F+90) ▲ Y=Y+J*sin(F+90) ▲
所求的XY坐标就是左右边桩的坐标，90度为正交。
第5步 Goto 0 Goto为无条件转移命令， Goto 0即为无条件转移到
标识符为0的程序运行段落中来

4800及4850计算器程序

4800及4850计算器常用程序ZBZS（坐标正算）E“X1=”:F“Y1=”:G“T1=LbI 0：Fixm：{AS}：G>360=>G=G-360 “T2=”:G DMSE=E+ScosG：“X2=”：F=F+SsinG: “Y2=”X1：起算点的X坐标Y1：起算点的Y坐标T1：起算边指向计算前进方向的坐标方位角A：计算前进方向的左角ZBFS（坐标反算）A“X1=”：B“Y1=”：C“X2=”：D“Y2=”：POL（C-A，D-B）：Fixm：J>0=>W=J:≠“W（1-2）=”“S（1-2）=”：XLYH（线路优化程序）A“XA=”：B“YA=”：C“TA=”：E“1÷RA，Z-，Y+”：P“KA=LbI1: ｛FQ｝:F“1÷RB，Z-，Y+”：Q“KB=LbI2:｛VDZ｝:V“K=”:Z“JIA-JIAO=”:D“WAIYI=(Z-,Y+)”: O=V:O<P=>Goto 3: ≠LbI3: “K-SMALL”LbI4: O>Q=>Goto 5: ≠LbI5: “JISUAN:A-TXY”LbI6: T=C+90+3.1415926×(2E+(F-E)(V-P)÷(Q-P))(V-P):Rad:W=A+∫(cos(C0+0.5×(2E+(F-E)(X-P)÷(Q-P))(X-P), P,V,4)+D×cos(T0+ Z0)Y=B+∫(sin(C0+0.5×(2E+(F-E)(X-P)÷(Q-P))(X-P)),P,V,4)+D×sin(T0+ Z0):LbI 7: O≤Q=>Goto 9:≠LbI 9: W“X=Y“Y=”T“T=”XA：起算点的X坐标YA：起算点的Y坐标TA：起算点的切线方位角1÷RA，Z-，Y+：起算点的曲率（1/R），左偏曲线为正，右偏曲线为负KA：起算点的里程1÷RB，Z-，Y+：终算点的曲率（1/R），左偏曲线为正，右偏曲线为负KB：终算点的里程JIA-JIAO：边桩计算时所要输的与线路切向的夹角，计算点与线路正交时输90度，斜交时则输与线路切向右侧所夹的斜交角度WAIYI=(Z-,Y+)：边桩计算时所要输的计算点到中线点的直线距离（即边桩的外移距离），计算点位于线路左侧时为负，位于线路右侧时为正程序说明：本程序计算线路坐标时要将线路分段（一般分直线段、第一缓和曲线段、圆曲线段、第二缓和曲线段、直线段以此循环）进行计算，例如：缓和曲线加圆曲线计算时，必须将曲线分三段进行计算，先输（起点要素）ZH点的坐标XA、YA，ZH到JD的方位角（即ZH点的切线方位角）TA，ZH 点的曲率1÷RA、里程KA，再输（终点要素）HY点的曲率1÷RB、里程KB，然后输入ZH到HY间任意点的里程便可计算；当所输计算点里程大于HY点里程时，程序显示“JISUAN:A-TXY”表示所要计算点超过第一缓和曲线计算范围，此时程序自动默认HY点的曲线要素为圆曲线计算段的起点要素，用户只需再输入（终点要素）YH点的曲率1÷RB、里程KB，就可计算圆曲线上任意点；当计算点里程大于YH 点里程时，程序显示“JISUAN:A-TXY”表示所要计算点超过圆曲线计算范围，程序自动默认YH点的曲线要素为第二缓和曲线计算段的起点要素，用户只需再输入（终点要素）HZ点的曲率1÷RB、里程KB，就可计算第二缓和曲线上任意点。

CASIO4800测量程序源代码

曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序一、程序功能本程序由一个主程序(TYQXJS)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。

本程序也可以在CASIO fx-4500P计算器及CASIO fx-4850P计算器上运行。

特别申明：(1). 适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm；(2). 在引用该核心计算部分时，请注明来源。

二、源程序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=>Goto 1：≠>Goto 2Δ←┘Lbl 1：{SZ}：SZ：W=Abs(S-O)：Prog "SUB1"：X"XS"=X◢Y"YS"=Y◢F"FS"=F-90◢Goto 1←┘Lbl 2：{XY}：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+W◢Z"Z"=Z◢Goto 22. 正算子程序(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：Lbl 0：Prog "SUB1"：L=T+QEW(C+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

4800计算器程序

4800-4850直线、曲线、螺旋线、匝道中边桩坐标计算程序直线、曲线、螺旋线、匝道中边桩坐标计算程序主程序名称（4800-4850）A“XA”：B“YB”：C“CA”：D“RA”：E“RB”：F“DKA”：G“DKB”换行D≠0=>D=1÷D：换行E≠0=>E=1÷E：换行Lb1 1:｛H，O，R｝换行H“DKI”：O“DL”：R“DR” 换行Lb1 0 换行P=（E－D）÷Abs（G－F）换行Abs（H－F）＞Abs（G－F）=>Q= Abs（G－F）：≠=>Q= Abs（H－F）：换行I=P×Q 换行J=C＋（I＋2D）×Q×90÷π 换行M=C＋（I÷8＋2D）×Q×45÷（4π）换行N=C＋（3I÷8＋2D）×Q×135÷（4π）换行U=C＋（5I÷8＋2D）×Q×225÷（4π）换行V=C＋（7I÷8＋2D）×Q×315÷（4π）换行K=C＋（I÷4＋2D）×Q×45÷（2π）换行W=C＋（I÷2＋2D）×Q×45÷π 换行Z=C＋（3I÷4＋2D）×Q×135÷（2π）换行X=A＋Q÷24×（cosC＋4×（cosM＋cosN＋cosU＋cosV）＋2×（cosK＋cosW＋cosZ）＋cosJ）换行Y=B＋Q÷24×（sinC＋4×（sinM＋sinN＋sinU＋sinV）＋2×（sinK＋sinW＋sinZ）＋sinJ）换行U=X＋O×cos（J－90）换行V=Y＋O×sin（J－90）换行W=X＋R×cos（J＋90）换行Z=Y＋R×sin（J＋90）换行Abs（H－F）＞Abs（G－F）=> Goto 3：换行H“DKI=” ◢ fx4850输入改为“DKI=”：H ◢X“X=” ◢ fx4850输入改为“X=”：X ◢Y“Y=” ◢ fx4850输入改为“Y=”：Y ◢U“XL=” ◢fx4850输入改为“XL=”：U◢V“YL=” ◢fx4850输入改为“YL=”：V◢W“XR=” ◢fx4850输入改为“XR=”：W◢Z“YR=” ◢fx4850输入改为“YR=”：Z◢Goto 1 换行Lb1 3 换行A=X：B=Y：D=E：F=G：C=J 换行｛E，G｝换行E“RB”：G“DKB”换行E≠0=>E=1÷E：换行Goto 0计算器显示及输入方法：XA？线型起算点的X坐标。

4850线路计算程序

4850线路计算程序1：已知线外任意点坐标，求对应线路里程在缓和曲线上，要计算任意里程的法线方向及任意宽度的边线坐标，非常简单。

但要计算任意一个已知坐标点，是对应哪一个里程法线方向上的点，就有一些困难。

很难推导一个这样的计算公式。

唯一的方法“渐进”，如果手工计算这可不是一个好方法。

但在有CASIO系列可编程计算器，如：FX-4500的情况下就变的非常简单了。

亦可用于直线和圆曲线的计算。

首先在缓和曲线上任选一点A为起始点，计算该点的坐标和切线方位角，通过坐标反算求起始点A与计算点B的方位角和距离，B点肯定对应A点切线方向上有一个垂足C点，把三点看成一个直角三角形，通过解直角三角形计算AC的距离，当该距离大于某一数值，如0。

001m，A点里程加AC的距离等于C点的里程，回到开始重新进入新一轮的计算，如果AC的距离小于某一规定值，则计算C点的里程与BC的距离即可。

求对应线路里程程序：主程序QLC (已知坐标求里程)Lb1 0：｛LＤE｝：Prog XH：Goto 0子程序：XH (循环)L1 Lb1 1L2 Norm：Prog LYYD：L3 PO1(D-X,E-Y)：W≤0=> W=W+360⊿L4 Z=W-I：A=V×cos Z：L=L+AL5 Abs A≥0.001=>Goto 1：≠=>B=V×sinZ：Fix 3：“FXJL=” ◢L6 L:Fix3:“DYLC=” ◢程序中字母代表D 任意点X坐标，E 任意点Y坐标，DYLC 对应里程，FXJL 中线法线距离。

程序中有坐标反算功能。

使用方法：只需输入计算点坐标、和较为接近的桩号。

桩号越接近计算速度越快2：逐桩坐标计算2.1编制方法：线路坐标程序是按照平曲线为单元，直线部分归属在曲线两端的方法，把整段路线分段装进数据库，根据桩号判断采用数据通过共用程序，进行任意点的坐标计算，在坐标转换示意土，第一直线段，是通过方位角和距离直接计算大地坐标，第一缓和曲线和圆曲线段，是先计算任意点切线支距和方位角然后转换大地坐标，第二缓和曲线段和直线段是先计算任意点切线支距和方位角。

CASIO(4850-4800)交点法计算平曲线要数各线元起点坐标及切线方位角(更新)1

平曲线要数、线元各起点坐标及切线方位角计算程序CASIO（4850-4800）原创版单位：中铁二十局集团第二工程有限公司编程：辜建君Defm12:LbI0:{ABCFΟZREK}:A”JD”:B”JDX”:C”JDY”:F”FJ”:Ο”ZJ:Z-,Y+”:Z:R:E”LS1”:K”LS2”:E<1=>E=1E-9⊿K<1=>K=1E-9⊿Z[2]=E÷2-EEE÷240RR:Z[6]=K÷2-KKK÷240RR: Z[7]= E-E^3÷40R^2+E^5÷3456R^4-E^7÷599040R^6+E^9÷17542600R^8:Z[8]= E^2÷6R-E^4÷336R^3+E^6÷42240R^5-E^8÷9676800R^7+E^10÷3530096640R^9: Z[9]= K-K^3÷40R^2+K^5÷3456R^4-K^7÷599040R^6+K^9÷17542600R^8:Z[10]= K^2÷6R-K^4÷336R^3+K^6÷42240R^5-K^8÷9676800R^7+K^10÷3530096640R^9: G=sin-1((Z[7]-Z[2])÷R): U= sin-1((Z[9]-Z[6])÷R): Z[1]=RcosG+Z[8]-R: Z[5]=RcosU+Z[10]-R ：V=(EE-KK)÷24R÷SinAbsΟ: ”T1”:Z[3]=(R+Z[1])tan(AbsΟ÷2)+Z[2]-V◢”T2”:Z[4]=(R+Z[5])tan(AbsΟ÷2)+Z[6]+V◢I=tan-1((R+Z[1])÷(Z[3]-Z[2])):“L”:L=AbsΟлR÷180+(E+K)÷2◢”E0”:U=(R+Z[1])÷sinI-R◢X=A-Z[3]:Y=X+E:E<1=>GOTO 1: ≠> GOTO 2 ⊿LbI 1:”ZY”:X◢“ZYX=”:N=B+Z[3]Cos(F+180) ◢“ZYY=”:Q=C+Z[3]sin(F+180) ◢“FWJ=”:F->DMS◢GOTO 3LbI 2:”ZH”:X◢“ZHX=”:N=B+Z[3]cos(F+180) ◢“ZHY=”:Q=C+Z[3]sin(F+180) ◢“FWJ=”:F->DMS◢“HY”:Y◢”HYX=”:N=B+(Z[3]-Z[7])cos(F+180)+Z[8]cos(F+90Z) ◢“HYY=”:Q=C+(Z[3]-Z[7])sin(F+180)+Z[8]sin(F+90Z) ◢J=F+ZG：J>360=>J=F+ZG-360:⊿J<0=>J=F+ZG+360: ⊿“FWJ=”:J->DMS◢LbI 3:”QZ”:S=X+(L-K-E)÷2+E◢D=Rsin(90 * (L-K-E)÷ЛR+sin-1((Z[7]-Z[2])÷R))+Z[2]:H=R(1-cos(*90 (L-K-E)÷ЛR+sin-1((Z[7]-Z[2])÷R)))+Z[1]“QZX=”:N= B+(Z[3]-D)cos(F+180)+Hcos(F+90Z) ◢“QZY=”:Q= C+(Z[3]-D)sin(F+180)+Hsin(F+90Z) ◢J=F+ZG+90Z(L-K-E)÷ЛR:J>360=>J= F+ZG+90Z(L-K-E)÷ЛR -360: ⊿J<0=>J=F+ZG+90Z(L-K-E)÷ЛR+360:⊿“FWJ=”:J >DMS◢LbI 4:M=X+L-K:P=X+L:K<1=>GOTO 5: ≠> GOTO 6⊿LbI 5:”YZ”:P◢“YZX=”:N=B+Z[4]com(F+Ο) ◢“YZY=”:Q=C+Z[4]sin(F+Ο) ◢W=F+Ο:W>360=>W=F+Ο-360: ⊿W<0=>W=F+Ο+360: ⊿”FWJ=”:W->DMS◢G0TO 0LbI 6:”YH”:M◢”YHX=”:N=B+(Z[4]-Z[9])cos(F+Ο)+Z[10]cos(F+Ο+90Z) ◢“YHY=”:Q=C+(Z[4]-Z[9])sin(F+Ο)+Z[10]sin(F+Ο+90Z) ◢J=F+ZG+180Z(L-K-E)÷ЛR:J>360=>J= F+ZG+180Z(L-K-E)÷ЛR -360: ⊿J<0=>J=F+ZG+180Z(L-K-E)÷ЛR+360:⊿“FWJ=”:J >DMS◢“HZ”:P◢“HZX=”:N=B+Z[4]com(F+Ο) ◢“HZY=”:Q=C+Z[4]sin(F+Ο) ◢W=F+Ο:W>360=>W=F+Ο-360: ⊿W<0=>W=F+Ο+360: ⊿”FWJ=”:W->DMS◢GOTO 0说明：一、便于区分数字零于字母O，在程序中用Ο代替字母O二、本程序适合CASIO4850、4800、4500计算器使用。

CASIO4800程序

O=514354.0958:Prog"HQZB":←
LbI D:L=430:R=10000:T=1953.99:A=173°27′41.99″:V=270:Q=-1:U=2121683.69:M=2604622.035:
O=514354.0958:Prog"YQZB":←
LbI E:A=333°43′58.99″:Q=1:U=2125126.98:R=10000:L=430:T=1953.99:M=2604622.035:
Л/R←
"B=":B▲
X=R COS B+H COS(A+180-QF)+M"JD(X)ZB"←
"X=":X▲
Y=R sin B+H sin(A+180-QF)+O"JD(Y)ZB"←
"Y=":Y▲
POI(X-K,Y-Z)▲
S=J+360←
"S=":S▲
Q=1=>B=B+180:△←
O=514354.0958:Prog"HQZB":←
YQZB
P=L"HQC"＾2/(24R)-L^4/(2688R^3):H=√((P+R)^2+(T"QXC"-L/2+L^3/(240R^2))^2):F=tan-1((P+R)
/(T-L/2+L^3/(240R^2))):LbI 1:{D,C,K,Z,W}:B=A+180+V+90LQ"-1,1"/(ЛR)+180Q(C"CD"-U"HY")/