;;; PLT Scheme Science Collection
;;; random-distributions/pareto-graphics.ss
;;; Copyright (c) 2004 M. Douglas Williams
;;;
;;; This library is free software; you can redistribute it and/or
;;; modify it under the terms of the GNU Lesser General Public
;;; License as published by the Free Software Foundation; either
;;; version 2.1 of the License, or (at your option) any later version.
;;;
;;; This library is distributed in the hope that it will be useful,
;;; but WITHOUT ANY WARRANTY; without even the implied warranty of
;;; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
;;; Lesser General Public License for more details.
;;;
;;; You should have received a copy of the GNU Lesser General Public
;;; License along with this library; if not, write to the Free
;;; Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
;;; 02111-1307 USA.
;;;
;;; -------------------------------------------------------------------
;;;
;;; This module implements the graphics for Pareto distributions.
;;;
;;; Version  Date      Description
;;; 1.0.0    09/28/04  Added contracts for functions.  Marked as ready
;;;                    for Release 1.0.  (Doug Williams)

(module pareto-graphics mzscheme
  
  (require (lib "contract.ss"))
  
  (provide/contract
   (pareto-plot
    (-> real? real? any)))

  (require "pareto.ss")
  (require (lib "plot.ss" "plot"))
  
  ;; pareto-plot: real x real -> void
  ;; Plot the pdf and cdf functions for a Pareto distribution with
  ;; exponent a and scale b.
  (define (pareto-plot a b)
    (plot (mix (line (lambda (x) (pareto-pdf x a b))
                     (color 'red))
               (line (lambda (x) (pareto-cdf x a b))
                     (color 'blue)))
          (x-min b) (x-max (+ b 20))
          (y-min 0) (y-max 1)
          (y-label "Density")
          (title "Pareto Distribution")))
  
)