!brkpow
! Eleonora Troja @ NASA/GSFC
! email: eleonora.troja_at_nasa.gov
! My lastmod: 11/09/10
! fit the light curve with a broken power law
! Parameters are: P1 first temporal index, P2 break time,
! P3 second temporal index, P4 normalization
: brkpow
VAR A
VAR B
VAR C
VAR D
X P2 <
IF
X P1 ^ P4 *
ELSE
P1 P3 - A STO
P2 A RCL ^  P4 *  B STO
X P3 ^ B RCL *
THEN
;