卡西欧FX5800---辛普森公式(万能公式)
1.Lbl 0：“XA=”?A： “YA=”?B: “CA=”?C：“1/RA=”?D：“1/RB=”?E：“DKA=”?F:“DKB=”?G
2.Lbl 1:“DKI=”?H：“DL=”?O：“DR=”?R：IF H＞G:THEN Goto0 IFEND
3.(E-D)/Abs(G-F)→P:Abs(H-F) →Q: P×Q→I:D+I→T
4.C+(I+2D)×Q×90/π→J
5. C+(I/4+2D)Q×45/(2π) →M: C+(3I/4+2D)Q×135/(2π) →N
6. C+(I/2+2D)Q×45/π→K
7. A+Q(cosC+4(cosM+cosN)+2cosK+cosJ)/12 →X
8.B+Q(sinC+4(sinM+sinN)+2sinK+sinJ)/12 →Y
9.”FW=”:J▲DMS ▲
10. “X=”:X▲
11. “Y=”:Y▲
12.“XL=”: X+Ocos(J-90) →U ▲
13.“YL=”: Y+Osin(J-90) →V ▲
14.“XR=”:X+Rcos(J+90) → W▲
15.“YR=”:Y+Rsin(J+90) → Z▲
16.Goto 1
程序结束
程序说明：
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”时，此时输入下一个曲线元起点的曲率和里程，即可计算下一个曲线中线及边线点位坐标。
使用该程序应注意事项；
该程序以前进方向为有意识，不可倒退计算；缓和段和圆曲线段应分开计算在计算圆曲线时应记下缓和
曲线尾的坐标方位角即“J”的角度；在计算第一段缓和曲线时曲率“1÷RA”输入0；在计算第二段缓和曲
线时“1÷RB” 输入0 。
（5）第一个0为零。
