!   eclipse_speed_brake.mac
!
mod_assemb_prep
!
!----------------------------------------------------------------------
!
ndiv=64
!
Ls=5.646
Lsmin=4.228
ra=15.402
rb=11.787
x1=1.25
x3=6.071
y3=1.207-.593
!
Lsmax=Lsmin+Ls
Lstep=Ls/ndiv
!
!----------------------------------------------------------------------
!
csys,0
!
view,'f'
!
*do,i,0,ndiv,1
    y1=Lsmin+i*Lstep
    y1%i%=y1
    c=sqrt((x3-x1)**2+(y3-y1)**2)
    c%i%=c
    !
    a_=acos(((c**2)+(rb**2)-(ra**2))/(2*c*rb))
    b_=acos(((c**2)+(ra**2)-(rb**2))/(2*c*ra))
    c_=acos(((ra**2)+(rb**2)-(c**2))/(2*ra*rb))
    e_=atan((y3-y1)/(x3-x1))
    !
    a_%i%=a_*dconv
    b_%i%=b_*dconv
    c_%i%=c_*dconv
    e_%i%=e_*dconv
    !
    phi%i%=(a_+b_+c_)*dconv
    theta%i%=(a_+e_)*dconv
    !
    x2=rb*cos(a_+e_)+x1
    y2=rb*sin(a_+e_)+y1
    !
    x2%i%=x2
    y2%i%=y2
    !
    k,,x1,y1,0
    k,,x2,y2,0
    k,,x3,y3,0
    !
    create_area
    !
    !voffst,am,-(i+1)
    !init_entity_num
    !k,,20,ky(km-2),0
    !lstr,km-2,km+1
    !init_entity_num
    !
    !e_msr%i%=(anglek(km-3,km-1,km))*dconv
    lplot
*enddo