double RSProbMax(X, N)
double X;
double N;
{
  double S = 0.0;
  double P = 1.0;
  double PP;
  int I, K;
  K = 1;
  I = 1;
  while (I <= N && I * X < 1) {
    P = P * (N + 1 - I) / I;
    K = -K;
    S -= P * K * pow(1.0 - I * X, N - 1);
    ++I;
  };
  return (S);
}


x86 assembly (clang, -O0):
<RSProbMax>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  0xece(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x1,-0x30(%rbp)
movl   $0x1,-0x2c(%rbp)
cvtsi2sdl -0x2c(%rbp),%xmm1
movsd  -0x10(%rbp),%xmm0
xor    %eax,%eax
ucomisd %xmm1,%xmm0
mov    %al,-0x31(%rbp)
jb     117a <RSProbMax+0x6a>
cvtsi2sdl -0x2c(%rbp),%xmm1
mulsd  -0x8(%rbp),%xmm1
movsd  0xe90(%rip),%xmm0
ucomisd %xmm1,%xmm0
seta   %al
mov    %al,-0x31(%rbp)
mov    -0x31(%rbp),%al
test   $0x1,%al
jne    118a <RSProbMax+0x7a>
jmp    1258 <RSProbMax+0x148>
movsd  -0x20(%rbp),%xmm0
movsd  0xe69(%rip),%xmm1
addsd  -0x10(%rbp),%xmm1
cvtsi2sdl -0x2c(%rbp),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x2c(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
xor    %eax,%eax
sub    -0x30(%rbp),%eax
mov    %eax,-0x30(%rbp)
movsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x30(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
cvtsi2sdl -0x2c(%rbp),%xmm0
movsd  -0x8(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xe01(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  -0x10(%rbp),%xmm1
movsd  0xdf0(%rip),%xmm2
subsd  %xmm2,%xmm1
call   1030 <pow@plt>
movaps %xmm0,%xmm2
movsd  -0x40(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1145 <RSProbMax+0x35>
movsd  -0x18(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<RSProbMax>:
push   %rbp
push   %r14
push   %rbx
sub    $0x40,%rsp
movapd %xmm0,%xmm2
movsd  %xmm1,0x10(%rsp)
ucomisd 0xed6(%rip),%xmm1
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x8(%rsp)
jb     123f <RSProbMax+0x12f>
movsd  0xebe(%rip),%xmm3
ucomisd %xmm2,%xmm3
jbe    123f <RSProbMax+0x12f>
movsd  0x10(%rsp),%xmm0
movapd %xmm0,%xmm1
addsd  %xmm3,%xmm1
movsd  %xmm1,0x28(%rsp)
movsd  0xea0(%rip),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x20(%rsp)
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x8(%rsp)
mov    $0xffffffff,%r14d
mov    $0x1,%eax
mov    $0x2,%ebp
movapd %xmm3,%xmm1
movapd %xmm3,%xmm4
movsd  %xmm2,0x18(%rsp)
nopw   0x0(%rax,%rax,1)
mov    %eax,%ebx
movsd  0x28(%rsp),%xmm0
subsd  %xmm1,%xmm0
mulsd  %xmm0,%xmm4
divsd  %xmm1,%xmm4
movsd  %xmm4,0x38(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %r14d,%xmm0
movsd  %xmm0,0x30(%rsp)
mulsd  %xmm2,%xmm1
movapd %xmm3,%xmm0
subsd  %xmm1,%xmm0
movsd  0x20(%rsp),%xmm1
call   1030 <pow@plt>
movsd  0x38(%rsp),%xmm4
movsd  0x30(%rsp),%xmm1
mulsd  %xmm4,%xmm1
mulsd  %xmm0,%xmm1
movsd  0x8(%rsp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rsp)
xorps  %xmm1,%xmm1
cvtsi2sd %ebp,%xmm1
movsd  0x10(%rsp),%xmm0
ucomisd %xmm1,%xmm0
jb     123f <RSProbMax+0x12f>
add    $0x1,%ebp
movapd %xmm1,%xmm0
movsd  0x18(%rsp),%xmm2
mulsd  %xmm2,%xmm0
movsd  0xdd1(%rip),%xmm3
ucomisd %xmm0,%xmm3
mov    %r14d,%eax
mov    %ebx,%r14d
ja     11a0 <RSProbMax+0x90>
movsd  0x8(%rsp),%xmm0
add    $0x40,%rsp
pop    %rbx
pop    %r14
pop    %rbp
ret

LLVM IR (-O0):
define double @RSProbMax(double %0, double %1) {
%3 = alloca double
%4 = alloca double
%5 = alloca double
%6 = alloca double
%7 = alloca double
%8 = alloca i32
%9 = alloca i32
store double %0, double* %3
store double %1, double* %4
store double 0.000000e+00, double* %5
store double 1.000000e+00, double* %6
store i32 1, i32* %9
store i32 1, i32* %8
br label %10
10: ; preds = %23, %2
%11 = load i32, i32* %8
%12 = sitofp i32 %11 to double
%13 = load double, double* %4
%14 = fcmp ole double %12, %13
br i1 %14, label %15, label %21
15: ; preds = %10
%16 = load i32, i32* %8
%17 = sitofp i32 %16 to double
%18 = load double, double* %3
%19 = fmul double %17, %18
%20 = fcmp olt double %19, 1.000000e+00
br label %21
21: ; preds = %15, %10
%22 = phi i1 [ false, %10 ], [ %20, %15 ]
br i1 %22, label %23, label %53
23: ; preds = %21
%24 = load double, double* %6
%25 = load double, double* %4
%26 = fadd double %25, 1.000000e+00
%27 = load i32, i32* %8
%28 = sitofp i32 %27 to double
%29 = fsub double %26, %28
%30 = fmul double %24, %29
%31 = load i32, i32* %8
%32 = sitofp i32 %31 to double
%33 = fdiv double %30, %32
store double %33, double* %6
%34 = load i32, i32* %9
%35 = sub nsw i32 0, %34
store i32 %35, i32* %9
%36 = load double, double* %6
%37 = load i32, i32* %9
%38 = sitofp i32 %37 to double
%39 = fmul double %36, %38
%40 = load i32, i32* %8
%41 = sitofp i32 %40 to double
%42 = load double, double* %3
%43 = fneg double %41
%44 = call double @llvm.fmuladd.f64(double %43, double %42, double 1.000000e+00)
%45 = load double, double* %4
%46 = fsub double %45, 1.000000e+00
%47 = call double @pow(double %44, double %46)
%48 = load double, double* %5
%49 = fneg double %39
%50 = call double @llvm.fmuladd.f64(double %49, double %47, double %48)
store double %50, double* %5
%51 = load i32, i32* %8
%52 = add nsw i32 %51, 1
store i32 %52, i32* %8
br label %10
53: ; preds = %21
%54 = load double, double* %5
ret double %54
}

LLVM IR (-O2):
define double @RSProbMax(double %0, double %1) {
%3 = fcmp oge double %1, 1.000000e+00
%4 = fcmp olt double %0, 1.000000e+00
%5 = and i1 %3, %4
br i1 %5, label %6, label %32
6: ; preds = %2
%7 = fadd double %1, 1.000000e+00
%8 = fadd double %1, -1.000000e+00
br label %9
9: ; preds = %6, %9
%10 = phi double [ 1.000000e+00, %6 ], [ %27, %9 ]
%11 = phi i32 [ 1, %6 ], [ %12, %9 ]
%12 = phi i32 [ -1, %6 ], [ %11, %9 ]
%13 = phi i32 [ 1, %6 ], [ %26, %9 ]
%14 = phi double [ 0.000000e+00, %6 ], [ %25, %9 ]
%15 = phi double [ 1.000000e+00, %6 ], [ %18, %9 ]
%16 = fsub double %7, %10
%17 = fmul double %15, %16
%18 = fdiv double %17, %10
%19 = sitofp i32 %12 to double
%20 = fneg double %10
%21 = tail call double @llvm.fmuladd.f64(double %20, double %0, double 1.000000e+00)
%22 = tail call double @pow(double %21, double %8)
%23 = fneg double %18
%24 = fmul double %19, %23
%25 = tail call double @llvm.fmuladd.f64(double %24, double %22, double %14)
%26 = add nuw nsw i32 %13, 1
%27 = sitofp i32 %26 to double
%28 = fcmp ole double %27, %1
%29 = fmul double %27, %0
%30 = fcmp olt double %29, 1.000000e+00
%31 = and i1 %28, %30
br i1 %31, label %9, label %32
32: ; preds = %9, %2
%33 = phi double [ 0.000000e+00, %2 ], [ %25, %9 ]
ret double %33
}