Mercurial > hg > octave-nkf
annotate scripts/elfun/acot.m @ 9665:1dba57e9d08d
use blas_trans_type for xgemm
author | Jaroslav Hajek <highegg@gmail.com> |
---|---|
date | Sat, 26 Sep 2009 10:41:07 +0200 |
parents | 16f53d29049f |
children | c776f063fefe |
rev | line source |
---|---|
9245 | 1 ## Copyright (C) 1994, 1996, 1997, 1999, 2000, 2004, 2005, 2006, 2007, |
2 ## 2008, 2009 John W. Eaton | |
2313 | 3 ## |
4 ## This file is part of Octave. | |
5 ## | |
6 ## Octave is free software; you can redistribute it and/or modify it | |
7 ## under the terms of the GNU General Public License as published by | |
7016 | 8 ## the Free Software Foundation; either version 3 of the License, or (at |
9 ## your option) any later version. | |
2313 | 10 ## |
11 ## Octave is distributed in the hope that it will be useful, but | |
12 ## WITHOUT ANY WARRANTY; without even the implied warranty of | |
13 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
14 ## General Public License for more details. | |
15 ## | |
16 ## You should have received a copy of the GNU General Public License | |
7016 | 17 ## along with Octave; see the file COPYING. If not, see |
18 ## <http://www.gnu.org/licenses/>. | |
2303 | 19 |
3321 | 20 ## -*- texinfo -*- |
3428 | 21 ## @deftypefn {Mapping Function} {} acot (@var{x}) |
9155
ad20b967e1c9
Update section 17.3 (Trigonometry) of arith.txi
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
22 ## Compute the inverse cotangent in radians for each element of @var{x}. |
ad20b967e1c9
Update section 17.3 (Trigonometry) of arith.txi
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
23 ## @seealso{cot, acotd} |
3321 | 24 ## @end deftypefn |
713 | 25 |
2314 | 26 ## Author: jwe |
27 | |
2311 | 28 function w = acot (z) |
713 | 29 |
30 if (nargin != 1) | |
6046 | 31 print_usage (); |
713 | 32 endif |
33 | |
4692 | 34 w = atan (1./z); |
713 | 35 |
2326 | 36 endfunction |
7385 | 37 |
38 %!test | |
39 %! rt2 = sqrt (2); | |
40 %! rt3 = sqrt (3); | |
41 %! x = [rt3, 1, rt3/3, 0, -rt3/3, -1, -rt3]; | |
42 %! v = [pi/6, pi/4, pi/3, pi/2, -pi/3, -pi/4, -pi/6]; | |
7742 | 43 %! assert(all (abs (acot (x) - v) < sqrt (eps))); |
7385 | 44 |
45 %!error acot (); | |
46 | |
47 %!error acot (1, 2); | |
48 |