Mercurial > hg > octave-lyh
view src/DLD-FUNCTIONS/betainc.cc @ 4776:adf8d68d7143 ss-2-1-54
[project @ 2004-02-16 20:32:20 by jwe]
author | jwe |
---|---|
date | Mon, 16 Feb 2004 20:32:20 +0000 |
parents | fa0d84558faa |
children | 9f7ef92b50b0 |
line wrap: on
line source
/* Copyright (C) 1997 John W. Eaton This file is part of Octave. Octave is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2, or (at your option) any later version. Octave 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 General Public License for more details. You should have received a copy of the GNU General Public License along with Octave; see the file COPYING. If not, write to the Free Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ #ifdef HAVE_CONFIG_H #include <config.h> #endif #include "lo-specfun.h" #include "defun-dld.h" #include "error.h" #include "gripes.h" #include "oct-obj.h" #include "utils.h" DEFUN_DLD (betainc, args, , "-*- texinfo -*-\n\ @deftypefn {Mapping Function} {} betainc (@var{x}, @var{a}, @var{b})\n\ Return the incomplete Beta function,\n\ @iftex\n\ @tex\n\ $$\n\ \\beta (a, b, x) = B (a, b)^{-1} \\int_0^x t^{(a-z)} (1-t)^{(b-1)} dt.\n\ $$\n\ @end tex\n\ @end iftex\n\ @ifinfo\n\ \n\ @smallexample\n\ x\n\ /\n\ betai (a, b, x) = beta (a, b)^(-1) | t^(a-1) (1-t)^(b-1) dt.\n\ /\n\ t=0\n\ @end smallexample\n\ @end ifinfo\n\ \n\ If x has more than one component, both @var{a} and @var{b} must be\n\ scalars. If @var{x} is a scalar, @var{a} and @var{b} must be of\n\ compatible dimensions.\n\ @end deftypefn") { octave_value retval; int nargin = args.length (); if (nargin == 3) { octave_value x_arg = args(0); octave_value a_arg = args(1); octave_value b_arg = args(2); if (x_arg.is_scalar_type ()) { double x = x_arg.double_value (); if (a_arg.is_scalar_type ()) { double a = a_arg.double_value (); if (! error_state) { if (b_arg.is_scalar_type ()) { double b = b_arg.double_value (); if (! error_state) retval = betainc (x, a, b); } else { Matrix b = b_arg.matrix_value (); if (! error_state) retval = betainc (x, a, b); } } } else { Matrix a = a_arg.matrix_value (); if (! error_state) { if (b_arg.is_scalar_type ()) { double b = b_arg.double_value (); if (! error_state) retval = betainc (x, a, b); } else { Matrix b = b_arg.matrix_value (); if (! error_state) retval = betainc (x, a, b); } } } } else { Matrix x = x_arg.matrix_value (); if (a_arg.is_scalar_type ()) { double a = a_arg.double_value (); if (! error_state) { if (b_arg.is_scalar_type ()) { double b = b_arg.double_value (); if (! error_state) retval = betainc (x, a, b); } else { Matrix b = b_arg.matrix_value (); if (! error_state) retval = betainc (x, a, b); } } } else { Matrix a = a_arg.matrix_value (); if (! error_state) { if (b_arg.is_scalar_type ()) { double b = b_arg.double_value (); if (! error_state) retval = betainc (x, a, b); } else { Matrix b = b_arg.matrix_value (); if (! error_state) retval = betainc (x, a, b); } } } } } else print_usage ("betainc"); return retval; } /* ;;; Local Variables: *** ;;; mode: C++ *** ;;; End: *** */