#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
scalar_t__ EOF;
double *dvector(int /*<<< orphan*/, int);
double fn(double);
int /*<<< orphan*/ free_dvector(double *, int /*<<< orphan*/, int);
int /*<<< orphan*/ pade(double *, int, float *);
int /*<<< orphan*/ printf(char *, ...);
double ratval(double, double *, int, int);
scalar_t__ scanf(char *, int *);
int main(void) {
  int j, k, n;
  float resid;
  double b, d, fac, x, *c, *cc;
  c = dvector(0, 100);
  cc = dvector(0, 100);
  for (;;) {
    printf("Enter n for PADE routine:\n");
    if (scanf("%d", &n) == EOF)
      break;
    fac = 1;
    for (j = 1; j <= 2 * n + 1; j++) {
      c[j - 1] = fac / ((double)j);
      cc[j - 1] = c[j - 1];
      fac = -fac;
    }
    pade(c, n, &resid);
    printf("Norm of residual vector= %16.8e\n", resid);
    printf("point, func. value, pade series, power series\n");
    for (j = 1; j <= 21; j++) {
      x = (j - 1) * 0.25;
      for (b = 0.0, k = 2 * n + 1; k >= 1; k--) {
        b *= x;
        b += cc[k - 1];
      }
      d = ratval(x, c, n, n);
      printf("%16.8f %16.8f %16.8f %16.8f\n", x, fn(x), d, b);
    }
  }
  free_dvector(cc, 0, 100);
  free_dvector(c, 0, 100);
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
movl   $0x0,-0x4(%rbp)
xor    %edi,%edi
mov    $0x64,%esi
call   1050 <dvector@plt>
mov    %rax,-0x40(%rbp)
xor    %edi,%edi
mov    $0x64,%esi
call   1050 <dvector@plt>
mov    %rax,-0x48(%rbp)
lea    0xe6a(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    0xe77(%rip),%rdi
lea    -0x10(%rbp),%rsi
call   1070 <scanf@plt>
mov    0x2e24(%rip),%rcx
cmp    (%rcx),%rax
jne    11d2 <main+0x62>
jmp    1385 <main+0x215>
movsd  0xe26(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x1,-0x8(%rbp)
mov    -0x8(%rbp),%eax
mov    -0x10(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     126d <main+0xfd>
movsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0x8(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x40(%rbp),%rax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%rax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8(%rbp)
jmp    11e6 <main+0x76>
mov    -0x40(%rbp),%rdi
mov    -0x10(%rbp),%esi
lea    -0x14(%rbp),%rdx
call   1090 <pade@plt>
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xda1(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
lea    0xdb4(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
movl   $0x1,-0x8(%rbp)
cmpl   $0x15,-0x8(%rbp)
jg     1380 <main+0x210>
mov    -0x8(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  0xd43(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
cmpl   $0x1,-0xc(%rbp)
jl     1325 <main+0x1b5>
movsd  -0x38(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x48(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0xc(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0xc(%rbp)
jmp    12e2 <main+0x172>
movsd  -0x38(%rbp),%xmm0
mov    -0x40(%rbp),%rdi
mov    -0x10(%rbp),%esi
mov    -0x10(%rbp),%edx
call   1060 <ratval@plt>
movsd  %xmm0,-0x28(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1030 <fn@plt>
movaps %xmm0,%xmm1
movsd  -0x50(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm2
movsd  -0x20(%rbp),%xmm3
lea    0xd13(%rip),%rdi
mov    $0x4,%al
call   1040 <printf@plt>
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8(%rbp)
jmp    12a9 <main+0x139>
jmp    119f <main+0x2f>
mov    -0x48(%rbp),%rdi
xor    %esi,%esi
mov    $0x64,%edx
call   1080 <free_dvector@plt>
mov    -0x40(%rbp),%rdi
xor    %esi,%esi
mov    $0x64,%edx
call   1080 <free_dvector@plt>
xor    %eax,%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
xor    %edi,%edi
mov    $0x64,%esi
call   1060 <dvector@plt>
mov    %rax,%r14
xor    %edi,%edi
mov    $0x64,%esi
call   1060 <dvector@plt>
mov    %rax,%rbx
lea    0xea6(%rip),%rdi
call   1030 <puts@plt>
lea    0xe59(%rip),%rdi
lea    0x8(%rsp),%rsi
call   1080 <scanf@plt>
mov    0x2e18(%rip),%rcx
cmp    (%rcx),%rax
jne    1204 <main+0x84>
mov    %rbx,%rdi
xor    %esi,%esi
mov    $0x64,%edx
call   1090 <free_dvector@plt>
mov    %r14,%rdi
xor    %esi,%esi
mov    $0x64,%edx
call   1090 <free_dvector@plt>
xor    %eax,%eax
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    $0xffffffff,%r12d
lea    0xe2b(%rip),%r15
movsd  0xde7(%rip),%xmm2
movsd  0xde7(%rip),%xmm3
jmp    126d <main+0xed>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
lea    0xe22(%rip),%rdi
call   1030 <puts@plt>
lea    0xdd5(%rip),%rdi
lea    0x8(%rsp),%rsi
call   1080 <scanf@plt>
mov    0x2d94(%rip),%rcx
cmp    (%rcx),%rax
movsd  0xda1(%rip),%xmm2
movsd  0xda1(%rip),%xmm3
je     11d5 <main+0x55>
mov    0x8(%rsp),%esi
lea    (%rsi,%rsi,1),%eax
add    $0x1,%eax
test   %eax,%eax
jle    1311 <main+0x191>
mov    %esi,%ecx
add    %ecx,%ecx
je     12f0 <main+0x170>
add    $0x2,%ecx
add    $0xfffffffffffffffe,%rcx
mov    $0x1,%eax
mov    $0x2,%edx
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
lea    -0x1(%rdx),%ebp
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
movapd %xmm2,%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,-0x8(%r14,%rax,8)
xorps  %xmm0,%xmm0
cvtsi2sd %edx,%xmm0
movsd  %xmm1,-0x8(%rbx,%rax,8)
movapd %xmm3,%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r14,%rax,8)
movsd  %xmm1,(%rbx,%rax,8)
add    $0x2,%rax
add    $0x2,%rdi
add    $0x2,%edx
cmp    %rcx,%rdi
jne    12a0 <main+0x120>
jmp    12f5 <main+0x175>
nopl   0x0(%rax,%rax,1)
mov    $0x1,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movapd %xmm2,%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,-0x8(%r14,%rax,8)
movsd  %xmm1,-0x8(%rbx,%rax,8)
mov    %r14,%rdi
lea    0xc(%rsp),%rdx
call   10a0 <pade@plt>
movss  0xc(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xcec(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
lea    0xd36(%rip),%rdi
call   1030 <puts@plt>
mov    $0x1,%r13d
jmp    13a9 <main+0x229>
nopw   0x0(%rax,%rax,1)
xorpd  %xmm1,%xmm1
movsd  %xmm1,0x20(%rsp)
movsd  %xmm0,0x10(%rsp)
mov    %r14,%rdi
mov    %edx,%esi
call   1070 <ratval@plt>
movsd  %xmm0,0x18(%rsp)
movsd  0x10(%rsp),%xmm0
call   1040 <fn@plt>
movapd %xmm0,%xmm1
mov    %r15,%rdi
movsd  0x10(%rsp),%xmm0
movsd  0x18(%rsp),%xmm2
movsd  0x20(%rsp),%xmm3
mov    $0x4,%al
call   1050 <printf@plt>
add    $0x1,%r13d
cmp    $0x16,%r13d
je     1230 <main+0xb0>
lea    -0x1(%r13),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0xc54(%rip),%xmm0
mov    0x8(%rsp),%edx
lea    (%rdx,%rdx,1),%eax
add    $0x1,%eax
test   %eax,%eax
jle    1350 <main+0x1d0>
lea    (%rdx,%rdx,1),%esi
or     $0x1,%esi
lea    -0x1(%rsi),%rbp
mov    %esi,%edi
and    $0x3,%edi
add    $0xfffffffffffffffe,%rsi
xorpd  %xmm1,%xmm1
mov    %rbp,%rcx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mulsd  %xmm0,%xmm1
mov    %ecx,%eax
addsd  (%rbx,%rax,8),%xmm1
mov    %rsi,%rax
add    $0xffffffffffffffff,%rcx
add    $0xffffffffffffffff,%rsi
add    $0xffffffffffffffff,%rdi
jne    13f0 <main+0x270>
cmp    $0x3,%ebp
jb     1354 <main+0x1d4>
cs nopw 0x0(%rax,%rax,1)
nop
mulsd  %xmm0,%xmm1
lea    (%rax,%r12,1),%ecx
add    $0x1,%ecx
addsd  (%rbx,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
lea    (%rax,%r12,1),%ecx
addsd  (%rbx,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
lea    -0x2(%rax),%rcx
mov    %ecx,%esi
addsd  (%rbx,%rsi,8),%xmm1
mulsd  %xmm0,%xmm1
lea    -0x3(%rax),%esi
addsd  (%rbx,%rsi,8),%xmm1
add    $0xfffffffffffffffc,%rax
cmp    $0x1,%rcx
ja     1420 <main+0x2a0>
jmp    1354 <main+0x1d4>

LLVM IR (-O0):
@.str = [27 x i8] c"Enter n for PADE routine:\0A\00"
@.str.1 = [3 x i8] c"%d\00"
@EOF = global i64 0
@.str.2 = [33 x i8] c"Norm of residual vector= %16.8e\0A\00"
@.str.3 = [47 x i8] c"point, func. value, pade series, power series\0A\00"
@.str.4 = [29 x i8] c"%16.8f %16.8f %16.8f %16.8f\0A\00"
define i32 @main() {
%1 = alloca i32
%2 = alloca i32
%3 = alloca i32
%4 = alloca i32
%5 = alloca float
%6 = alloca double
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca double*
%11 = alloca double*
store i32 0, i32* %1
%12 = call double* @dvector(i32 0, i32 100)
store double* %12, double** %10
%13 = call double* @dvector(i32 0, i32 100)
store double* %13, double** %11
br label %14
14: ; preds = %105, %0
%15 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([27 x i8], [27 x i8]* @.str, i64 0, i64 0))
%16 = call i64 @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.1, i64 0, i64 0), i32* %4)
%17 = load i64, i64* @EOF
%18 = icmp eq i64 %16, %17
br i1 %18, label %19, label %20
19: ; preds = %14
br label %106
20: ; preds = %14
store double 1.000000e+00, double* %8
store i32 1, i32* %2
br label %21
21: ; preds = %50, %20
%22 = load i32, i32* %2
%23 = load i32, i32* %4
%24 = mul nsw i32 2, %23
%25 = add nsw i32 %24, 1
%26 = icmp sle i32 %22, %25
br i1 %26, label %27, label %53
27: ; preds = %21
%28 = load double, double* %8
%29 = load i32, i32* %2
%30 = sitofp i32 %29 to double
%31 = fdiv double %28, %30
%32 = load double*, double** %10
%33 = load i32, i32* %2
%34 = sub nsw i32 %33, 1
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %32, i64 %35
store double %31, double* %36
%37 = load double*, double** %10
%38 = load i32, i32* %2
%39 = sub nsw i32 %38, 1
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %37, i64 %40
%42 = load double, double* %41
%43 = load double*, double** %11
%44 = load i32, i32* %2
%45 = sub nsw i32 %44, 1
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double, double* %43, i64 %46
store double %42, double* %47
%48 = load double, double* %8
%49 = fneg double %48
store double %49, double* %8
br label %50
50: ; preds = %27
%51 = load i32, i32* %2
%52 = add nsw i32 %51, 1
store i32 %52, i32* %2
br label %21
53: ; preds = %21
%54 = load double*, double** %10
%55 = load i32, i32* %4
%56 = call i32 @pade(double* %54, i32 %55, float* %5)
%57 = load float, float* %5
%58 = fpext float %57 to double
%59 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([33 x i8], [33 x i8]* @.str.2, i64 0, i64 0), double %58)
%60 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([47 x i8], [47 x i8]* @.str.3, i64 0, i64 0))
store i32 1, i32* %2
br label %61
61: ; preds = %102, %53
%62 = load i32, i32* %2
%63 = icmp sle i32 %62, 21
br i1 %63, label %64, label %105
64: ; preds = %61
%65 = load i32, i32* %2
%66 = sub nsw i32 %65, 1
%67 = sitofp i32 %66 to double
%68 = fmul double %67, 2.500000e-01
store double %68, double* %9
store double 0.000000e+00, double* %6
%69 = load i32, i32* %4
%70 = mul nsw i32 2, %69
%71 = add nsw i32 %70, 1
store i32 %71, i32* %3
br label %72
72: ; preds = %87, %64
%73 = load i32, i32* %3
%74 = icmp sge i32 %73, 1
br i1 %74, label %75, label %90
75: ; preds = %72
%76 = load double, double* %9
%77 = load double, double* %6
%78 = fmul double %77, %76
store double %78, double* %6
%79 = load double*, double** %11
%80 = load i32, i32* %3
%81 = sub nsw i32 %80, 1
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds double, double* %79, i64 %82
%84 = load double, double* %83
%85 = load double, double* %6
%86 = fadd double %85, %84
store double %86, double* %6
br label %87
87: ; preds = %75
%88 = load i32, i32* %3
%89 = add nsw i32 %88, -1
store i32 %89, i32* %3
br label %72
90: ; preds = %72
%91 = load double, double* %9
%92 = load double*, double** %10
%93 = load i32, i32* %4
%94 = load i32, i32* %4
%95 = call double @ratval(double %91, double* %92, i32 %93, i32 %94)
store double %95, double* %7
%96 = load double, double* %9
%97 = load double, double* %9
%98 = call double @fn(double %97)
%99 = load double, double* %7
%100 = load double, double* %6
%101 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([29 x i8], [29 x i8]* @.str.4, i64 0, i64 0), double %96, double %98, double %99, double %100)
br label %102
102: ; preds = %90
%103 = load i32, i32* %2
%104 = add nsw i32 %103, 1
store i32 %104, i32* %2
br label %61
105: ; preds = %61
br label %14
106: ; preds = %19
%107 = load double*, double** %11
%108 = call i32 @free_dvector(double* %107, i32 0, i32 100)
%109 = load double*, double** %10
%110 = call i32 @free_dvector(double* %109, i32 0, i32 100)
ret i32 0
}

LLVM IR (-O2):
@.str.1 = [3 x i8] c"%d\00"
@EOF = global i64 0
@.str.2 = [33 x i8] c"Norm of residual vector= %16.8e\0A\00"
@.str.4 = [29 x i8] c"%16.8f %16.8f %16.8f %16.8f\0A\00"
@str = [26 x i8] c"Enter n for PADE routine:\00"
@str.5 = [46 x i8] c"point, func. value, pade series, power series\00"
define i32 @main() {
%1 = alloca i32
%2 = alloca float
%3 = bitcast i32* %1 to i8*
%4 = bitcast float* %2 to i8*
%5 = tail call double* @dvector(i32 0, i32 100)
%6 = tail call double* @dvector(i32 0, i32 100)
%7 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([26 x i8], [26 x i8]* @str, i64 0, i64 0))
%8 = call i64 @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.1, i64 0, i64 0), i32* %1)
%9 = load i64, i64* @EOF
%10 = icmp eq i64 %8, %9
br i1 %10, label %122, label %16
11: ; preds = %115
%12 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([26 x i8], [26 x i8]* @str, i64 0, i64 0))
%13 = call i64 @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.1, i64 0, i64 0), i32* %1)
%14 = load i64, i64* @EOF
%15 = icmp eq i64 %13, %14
br i1 %15, label %122, label %16
16: ; preds = %0, %11
%17 = load i32, i32* %1
%18 = shl i32 %17, 1
%19 = or i32 %18, 1
%20 = icmp slt i32 %19, 1
br i1 %20, label %53, label %21
21: ; preds = %16
%22 = icmp eq i32 %18, 0
br i1 %22, label %45, label %23
23: ; preds = %21
%24 = add i32 %18, 2
%25 = zext i32 %24 to i64
%26 = add nsw i64 %25, -2
br label %27
27: ; preds = %27, %23
%28 = phi i64 [ 1, %23 ], [ %42, %27 ]
%29 = phi i64 [ 0, %23 ], [ %43, %27 ]
%30 = trunc i64 %28 to i32
%31 = sitofp i32 %30 to double
%32 = fdiv double 1.000000e+00, %31
%33 = add nsw i64 %28, -1
%34 = getelementptr inbounds double, double* %5, i64 %33
store double %32, double* %34
%35 = getelementptr inbounds double, double* %6, i64 %33
store double %32, double* %35
%36 = trunc i64 %28 to i32
%37 = add i32 %36, 1
%38 = sitofp i32 %37 to double
%39 = fdiv double -1.000000e+00, %38
%40 = getelementptr inbounds double, double* %5, i64 %28
store double %39, double* %40
%41 = getelementptr inbounds double, double* %6, i64 %28
store double %39, double* %41
%42 = add nuw nsw i64 %28, 2
%43 = add i64 %29, 2
%44 = icmp eq i64 %43, %26
br i1 %44, label %45, label %27
45: ; preds = %27, %21
%46 = phi i64 [ 1, %21 ], [ %42, %27 ]
%47 = trunc i64 %46 to i32
%48 = sitofp i32 %47 to double
%49 = fdiv double 1.000000e+00, %48
%50 = add nsw i64 %46, -1
%51 = getelementptr inbounds double, double* %5, i64 %50
store double %49, double* %51
%52 = getelementptr inbounds double, double* %6, i64 %50
store double %49, double* %52
br label %53
53: ; preds = %45, %16
%54 = call i32 @pade(double* %5, i32 %17, float* %2)
%55 = load float, float* %2
%56 = fpext float %55 to double
%57 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([33 x i8], [33 x i8]* @.str.2, i64 0, i64 0), double %56)
%58 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([46 x i8], [46 x i8]* @str.5, i64 0, i64 0))
br label %59
59: ; preds = %53, %115
%60 = phi i32 [ 1, %53 ], [ %120, %115 ]
%61 = add nsw i32 %60, -1
%62 = sitofp i32 %61 to double
%63 = fmul double %62, 2.500000e-01
%64 = load i32, i32* %1
%65 = shl i32 %64, 1
%66 = or i32 %65, 1
%67 = icmp sgt i32 %66, 0
br i1 %67, label %68, label %115
68: ; preds = %59
%69 = or i32 %65, 1
%70 = zext i32 %69 to i64
%71 = add nsw i64 %70, -1
%72 = and i64 %70, 3
br label %73
73: ; preds = %73, %68
%74 = phi i64 [ %70, %68 ], [ %78, %73 ]
%75 = phi double [ 0.000000e+00, %68 ], [ %82, %73 ]
%76 = phi i64 [ 0, %68 ], [ %83, %73 ]
%77 = fmul double %63, %75
%78 = add nsw i64 %74, -1
%79 = and i64 %78, 4294967295
%80 = getelementptr inbounds double, double* %6, i64 %79
%81 = load double, double* %80
%82 = fadd double %77, %81
%83 = add i64 %76, 1
%84 = icmp eq i64 %83, %72
br i1 %84, label %85, label %73
85: ; preds = %73
%86 = icmp ult i64 %71, 3
br i1 %86, label %115, label %87
87: ; preds = %85, %87
%88 = phi i64 [ %109, %87 ], [ %78, %85 ]
%89 = phi double [ %113, %87 ], [ %82, %85 ]
%90 = fmul double %63, %89
%91 = add i64 %88, 4294967295
%92 = and i64 %91, 4294967295
%93 = getelementptr inbounds double, double* %6, i64 %92
%94 = load double, double* %93
%95 = fadd double %90, %94
%96 = fmul double %63, %95
%97 = add i64 %88, 4294967294
%98 = and i64 %97, 4294967295
%99 = getelementptr inbounds double, double* %6, i64 %98
%100 = load double, double* %99
%101 = fadd double %96, %100
%102 = fmul double %63, %101
%103 = add nsw i64 %88, -3
%104 = and i64 %103, 4294967295
%105 = getelementptr inbounds double, double* %6, i64 %104
%106 = load double, double* %105
%107 = fadd double %102, %106
%108 = fmul double %63, %107
%109 = add nsw i64 %88, -4
%110 = and i64 %109, 4294967295
%111 = getelementptr inbounds double, double* %6, i64 %110
%112 = load double, double* %111
%113 = fadd double %108, %112
%114 = icmp ugt i64 %103, 1
br i1 %114, label %87, label %115
115: ; preds = %85, %87, %59
%116 = phi double [ 0.000000e+00, %59 ], [ %82, %85 ], [ %113, %87 ]
%117 = call double @ratval(double %63, double* %5, i32 %64, i32 %64)
%118 = call double @fn(double %63)
%119 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([29 x i8], [29 x i8]* @.str.4, i64 0, i64 0), double %63, double %118, double %117, double %116)
%120 = add nuw nsw i32 %60, 1
%121 = icmp eq i32 %120, 22
br i1 %121, label %11, label %59
122: ; preds = %11, %0
%123 = call i32 @free_dvector(double* %6, i32 0, i32 100)
%124 = call i32 @free_dvector(double* %5, i32 0, i32 100)
ret i32 0
}