1 /* This is the example given and commented on the MPFR web site:
2 * http://www.mpfr.org/sample.html
6 Copyright 1999, 2000, 2001, 2002, 2003, 2004, 2006, 2007, 2008, 2009, 2010, 2011 Free Software Foundation, Inc.
7 Contributed by the Arenaire and Cacao projects, INRIA.
9 This file is part of the GNU MPFR Library.
11 The GNU MPFR Library is free software; you can redistribute it and/or modify
12 it under the terms of the GNU Lesser General Public License as published by
13 the Free Software Foundation; either version 3 of the License, or (at your
14 option) any later version.
16 The GNU MPFR Library is distributed in the hope that it will be useful, but
17 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
18 or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
19 License for more details.
21 You should have received a copy of the GNU Lesser General Public License
22 along with the GNU MPFR Library; see the file COPYING.LESSER. If not, see
23 http://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc.,
24 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA.
38 mpfr_set_d (t
, 1.0, GMP_RNDD
);
40 mpfr_set_d (s
, 1.0, GMP_RNDD
);
42 for (i
= 1; i
<= 100; i
++)
44 mpfr_mul_ui (t
, t
, i
, GMP_RNDU
);
45 mpfr_set_d (u
, 1.0, GMP_RNDD
);
46 mpfr_div (u
, u
, t
, GMP_RNDD
);
47 mpfr_add (s
, s
, u
, GMP_RNDD
);
50 mpfr_out_str (stdout
, 10, 0, s
, GMP_RNDD
);