financial functions

gregg.irwin (Gregg Irwin) — Fri, 13 Jan 2006 04:05:22 GMT


REBOL [

   Title:   "Financial Library"

   Date:    17-Apr-2002

   Name:    'fin

   Version: 0.0.2

   File:    %fin.r

   ;Home:    http://www.pointillistic.com/

   Author:  "Gregg Irwin"

   Email:   gregg@pointillistic.com

   ;Owner:   "Pointillistic Software"

   ;Rights:  "Copyright © Pointillistic Software 2002. All Rights Reserved."

   ;Tabs:    4

   ;Need:    0.1.4

   ;Language: 'English

   ;Charset:  'ANSI

   ;Category: [script 1]

   Purpose: {Financial module for standard library.}



   Comment: {Prototype stage.}



   History: [

      0.0.1 [13-Mar-2002 "First whack" "Gregg"]

      0.0.2 [17-Apr-2002 {Added NPV and a functional, if ugly, IRR.

      BRACKET-FN supports IRR.} "Gregg"]

   ]



   Example: {do %fin.r}

]



fin: make lib-kernel [



    ; Future Value

    fv: func [

        {Computes future value based on present value.}

        value   [number! money!] "Present value"

        rate    [number!] "Interest Rate"

        periods [number!] "Number of time periods"

    ][

        1 + rate ** periods * value

    ]

    ;lib/fin/fv 1000 .1 12

    ;lib/fin/fv 1000 .1 24

    ;lib/fin/fv 1000 .1 36





    ; Compound Interest

    interest: func [

        {Computes accumulated compound interest.}

        value   [number! money!] "Present value"

        rate    [number!] "Interest Rate"

        periods [number!] "Number of time periods"

    ][

        1 + rate ** periods - 1 * value

    ]

    ;lib/fin/interest 1000 .1 12

    ;lib/fin/interest 1000 .1 24

    ;lib/fin/interest 1000 .1 36





    ; Interest Rate

    rate: func [

        {Computes interest rate for a given future value.}

        value        [number! money!] "Present value"

        future-value [number! money!] "Future value"

        periods      [number!] "Number of time periods"

    ][

        future-value / value ** (1 / periods) - 1

    ]

    ;lib/fin/rate 1000 future-value 1000 .1 12 12

    ;lib/fin/rate 1000 future-value 1000 .1 24 24

    ;lib/fin/rate 1000 future-value 1000 .1 36 36





    ; Present Value

    pv: func [

        {Computes present value based on future value.}

        value   [number! money!] "Future value"

        rate    [number!] "Interest Rate"

        periods [number!] "Number of time periods"

    ][

        1 + rate ** (negate periods) * value

        ; - or -

        ;value / (1 + rate ** periods)

    ]

    ;lib/fin/pv future-value 1000 .1 12 .1 12

    ;lib/fin/pv future-value 1000 .1 24 .1 24

    ;lib/fin/pv future-value 1000 .1 36 .1 36



    npv: func [

        {Computes net present value based on future cash flow values

        using continuous discounting.}

        values  [any-block!] "Cash flow values"

        rate    [number!] "Interest Rate"

        /local result

    ][

        result: $0.0

        repeat i length? values [

            result: result + pv values/:i rate i

        ]

        result

    ]

    ;lib/fin/npv [1000 2000 3000] .1

    ;lib/fin/npv [-10000, 3000, 4200, 6800] .1



    ;[-$70'000 $12'000 $15'000 $18'000 $21'000 $26'000]



    ; Time Periods

    nper: func [

        {Computes the number of time periods from the present value

        to a given future value.}

        value        [number! money!] "Present value"

        future-value [number! money!] "Future value"

        rate         [number!] "Interest Rate"

    ][

        ; 2 versions. Both work.

        (log-e (future-value / value)) / (log-e (rate + 1))

        ;divide log-e divide future-value value  log-e add rate 1

    ]

    ;lib/fin/nper 1000 future-value 1000 .1 12 .1

    ;lib/fin/nper 1000 future-value 1000 .1 24 .1

    ;lib/fin/nper 1000 future-value 1000 .1 36 .1



    irr: func [

        {Computes the internal rate of return based on future cash flow values.

        The internal rate of return of a cash flow is the interest rate that

        makes the present value of a cash flow equal to zero. THIS IS A QUICKLY

        HACKED VERSION AND NEEDS TO BE CLEANED UP! HERE FOR TESTING.}

        values  [any-block!] "Cash flow values"

        /guess

            guess-val

        /local rates rate npv-result result

    ][

        tries: 50

        ACCURACY: $.000001

        range: none

        rtb: dx: x-mid: 0.0

        f-mid: $0



        if not guess [guess-val: .2]



        if range: bracket-fn :npv values 0 guess-val [

            either ((npv values range/1) < $0) [

                rtb: range/1

                dx:  range/2 - range/1

            ][

                rtb: range/2

                dx:  range/1 - range/2

            ]

            loop tries [

            	dx: dx * 0.5

            	x-mid: rtb + dx

            	f-mid: npv values x-mid

            	if (f-mid <= $0) [rtb: x-mid]

            	if any [

                    (abs f-mid < ACCURACY)

                    (abs dx < second ACCURACY)

                ] [return x-mid]

            ]

            return none

        ]

    ]



    bracket-fn: func [

        {Adapted from Numerical Recipes. Given a user supplied function and a

        given 'guess' range (x1 to x2) it expands the range geometrically until

        a root is bracketed by the returned values x1 and x2, in which case

        the range is returned as a block of those two values. If it doesn't

        succeed after 50 tries, it returns false. This was built to support

        IRR so it's designed to call NPV, which takes one arg. It should be

        generalized.}

        fn arg x1 x2

        /local factor tries f1 f2

    ][

        tries: 50

        factor: 1.6

        f1: fn arg x1

        f2: fn arg x2

        loop tries [

            ; We have to use SECOND on f2 because it's a money! value

            ; and the numeric component is returned by SECOND.

            if negative? multiply f1 second f2 [

                return reduce [x1 x2]

            ]

            either (abs f1) < (abs f2) [

                x1: x1 + (factor * (x1 - x2))

                f1: fn arg x1

            ][

                x2: x2 + (factor * (x2 - x1))

                f2: fn arg x2

            ]

        ]

        false

    ]

    ;bracket-fn get in lib/fin 'npv values 0 .2

    ;values: [-$70'000 $12'000 $15'000 $18'000 $21'000 $26'000]

    ;lib/fin/irr values

    ;lib/fin/irr copy/part values 5

    ;lib/fin/irr copy/part values 3

    ;lib/fin/irr/guess copy/part values 3 -.1

    ;values: [-10 1.8 1.8 1.8 1.8 1.8 1.8 1.8 2.8]

    ;lib/fin/irr values

    ;values: [0 -1073 -1459 -1364 -1247 -1110 31789]

    ;lib/fin/irr values



;     irr: func [

;         {Computes the internal rate of return based on future cash flow values.

;         The internal rate of return of a cash flow is the interest rate that

;         makes the present value of a cash flow equal to zero. THIS IS A QUICKLY

;         HACKED VERSION AND NEEDS TO BE CLEANED UP! HERE FOR TESTING.}

;         values  [any-block!] "Cash flow values"

;         /guess

;             guess-val

;         /local rates rate npv-result result

;     ][

;         ACCURACY: $.000001

;         tries: 50

;         range: none

;         rtb: dx: x-mid: 0.0

;         f-mid: $0

;         if not guess [guess-val: .2]

;         if range: bracket-fn :npv values 0 guess-val [

;             f: npv values range/1

;             either (f < $0) [

;                 rtb: range/1

;                 dx: range/2 - range/1

;             ][

;                 rtb: range/2

;                 dx: range/1 - range/2

;             ]

;             loop tries [

;             	dx: dx * 0.5

;             	x-mid: rtb + dx

;             	f-mid: npv values x-mid

;             	if (f-mid <= $0) [rtb: x-mid]

;             	if any [

;                     ((abs f-mid) < ACCURACY)

;                     ((abs dx) < (second ACCURACY))

;                 ] [return x-mid]

;             ]

;             return none

;         ]

;     ]

;

;     bracket-fn: func [

;         {Adapted from Numerical Recipes. Given a user supplied function and a

;         given 'guess' range (x1 to x2) it expands the range geometrically until

;         a root is bracketed by the returned values x1 and x2, in which case

;         the range is returned as a block of those two values. If it doesn't

;         succeed after 50 tries, it returns false. This was built to support

;         IRR so it's designed to call NPV, which takes one arg. It should be

;         generalized.}

;         fn arg x1 x2

;         /local factor tries f1 f2

;     ][

;         tries: 50

;         factor: 1.6

;         f1: fn arg x1

;         f2: fn arg x2

;         loop tries [

;             ; We have to use SECOND on f2 because it's a money! value

;             ; and the numeric component is returned by SECOND.

;             if negative? multiply f1 second f2 [

;                 return reduce [x1 x2]

;             ]

;             either (abs f1) < (abs f2) [

;                 x1: x1 + (factor * (x1 - x2))

;                 f1: fn arg x1

;             ][

;                 x2: x2 + (factor * (x2 - x1))

;                 f2: fn arg x2

;             ]

;         ]

;         false

;     ]



]
DZone Snippets: FV code

financial functions
