Newton's method for quadratics involving continued fractions.
[frac.git] / epow.c
blob4dde0cb638aa0f2bbdf6bc35efd57e97e15649db
1 #include <stdio.h>
2 #include <stdlib.h>
3 #include <string.h>
4 #include <gmp.h>
5 #include "cf.h"
7 int main() {
8 mpz_t z;
9 mpz_init(z);
10 cf_t x, conv;
12 mpz_set_si(z, 3);
13 x = cf_new_tanh(z);
15 printf("tanh 3 = ");
16 conv = cf_new_cf_to_decimal(x);
17 for (int i = 0; i <= 5000; i++) {
18 cf_get(z, conv);
19 gmp_printf("%Zd", z);
20 if (!(i % 5)) putchar(' ');
21 if (!(i % 50)) putchar('\n');
23 cf_free(conv);
24 cf_free(x);
26 mpz_set_ui(z, 2);
27 x = cf_new_epow(z);
29 printf("e^2 = ");
30 conv = cf_new_cf_to_decimal(x);
31 for (int i = 0; i <= 5000; i++) {
32 cf_get(z, conv);
33 gmp_printf("%Zd", z);
34 if (!(i % 5)) putchar(' ');
35 if (!(i % 50)) putchar('\n');
37 cf_free(conv);
38 cf_free(x);
40 mpz_clear(z);
41 return 0;