; ************************* 
; Math/Utilities Library 17
; *************************
; Romberg Integration 
; (Original TI-59)


#EMULATOR c:\ti-59\mu\17\MU-17.t59

#LISTING c:\ti-59\mu\17\MU-17.LST

#AZEMULATOR c:\ti-59\mu\17\MU-17.TI


lbl a                                     ;000
sto 10 rtn

lbl b                                     ;005
sto 11 rtn

lbl c                                     ;010
cp x%t rtn

lbl d                                     ;015
inv stf 7 rcl 10 pgm 0 a' sto 18
rcl 11 pgm 0 a' sum 18 ( ( rcl 11
- rcl 10 ) / 2 ) prd 18 sto 19
( ( rcl 11 + rcl 10 ) / 2 ) pgm 0
a' prd 19 ( rcl 18 / 2 ) sum 19 ( (
4 * rcl 19 - rcl 18 ) / 3 ) sto 20
sto 18 2 sto 9 

jump4:                                     ;093
op 28 rc* 8 op 19
iff 7 jump1 op 29 ( ( rcl 11 - rcl 10 )
/ 2 yx ( rcl 9 - 1 ) sto 17 ) sto 13
( +/- + rcl 10 ) sto 14 ( 2 yx ( rcl 9
- 2 ) ) dms sto 15 0 sto 16 

jump2:                                     ;151
rcl 13 sum 14 
sum 14 rcl 14 pgm 0 a' sum 16 dsz 15 jump2
1 9 sto 8 ( rcl 19 / 2 + rcl 13 * rcl 16 )
ex* 8 

jump3:                                     ;186
( ( +/- + rc* 8 * ( 4 yx ( rcl 8
- 1 8 ) ) sto 12 ) / ( rcl 12 - 1 ) )
op 28 ex* 8 dsz 17 jump3 ( ( rc* 8 - exc 18
) / rcl 18 ) abs ge jump4 

jump1:                                      ;239
x%t rcl 18 rtn


#end