辛普森公式

辛普森程序

程序名：ZBJS

Lb1 0:M" (ZHI JING)X": N" (ZHI JING)Y": U" (HOU SHI)X" : V" (HOU SHI)Y":I=0:J=0:Pol(U-M,V-N):Z [7]=I: "JU LI": Z [7]◢

W=J:J<0⇒J=J+360Δ" FWJ":J→DMS◢

Lb1 1:｛EG｝: A"XA": B"YA": C"CA"→DMS D"1÷RA": E"1÷RB":F"DK A": G"DKB"↵

Lb1 2: ｛HOR｝: H"DKI": O"DL": R"DR":H>G⇒Goto3Δ↵

P=(E-D) ÷Abs(G-F): Q= Abs(H-F): S= P×Q: T=D+S ↵

L=C+( S+2D) Q ×90÷π: L＜0⇒L = L+360Δ"J=":L→DMS◢

Z [1]=C+( S ÷4+2D) Q×45÷(2π): Z [2]=C+(3 S ÷4+2D) Q ×135÷(2π)↵K=C+( S÷2+2D)Q ×45÷π↵

"X=": X=A+(Q÷12×(cosC+4(cos Z [1] + cos Z [2])+2 cos K + cosL)◢"Y=": Y=B+(Q÷12×(sinC+4(sin Z [1] + sinZ [2]) +2 sin K +sin L)◢

I=0:J=0:Pol(X-M,Y-N):Z=J-W:Z＜0⇒Z=Z+360Δ"Z(JIAO DU)":Z→D MS◢

"Z(JU LI)":I◢

"XL=": Z [3]=X+Ocos(L-90)◢

"YL=": Z [4]= Y+Osin (L-90)◢

I=0:J=0:Pol(Z [3]-M, Z [4]-N):Z=J-W:Z＜0⇒Z=Z+360Δ"L(JIAO DU)": Z→DMS◢

"L(JU LI)":I◢

"XR=": Z [5]=X+Rcos(L+90)◢

"YR=": Z [6]= Y+Rsin (L+90)◢

I=0:J=0:Pol(Z [5]-M, Z [6]-N):Z=J-W:Z＜0⇒Z=Z+360Δ"R(JIAO DU)": Z→DMS◢

"R(JU LI)":I◢

(A= X : B= Y : D=E: C= L : F=H:) Goto2↵

Lb1 3:A= X : B= Y : D=E: C= L : F=G: Goto1↵

注：括号闪烁（）中内容可省略。

半径输1÷半径例如：半径为1000输入形式为1÷1000。右曲为正，左曲为负。

XA……………….起点X YA…………………起点Y CA…………..切线方位角1÷RA……………起点半径

1÷RA…………..起点半径DKA…………….起点里程DKB…………...终点里程DKI…………...测设点里程