#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 */
int /*<<< orphan*/ back_substitution(double *, double *, int, int);
int /*<<< orphan*/ gauss_jordan(double *, int, int);
int /*<<< orphan*/ gaussian_elimination(double *, int, int);
int /*<<< orphan*/ printf(char *, ...);
int main() {
  double a[3][4] = {
      {3.0, 2.0, -4.0, 3.0}, {2.0, 3.0, 3.0, 15.0}, {5.0, -3.0, 1.0, 14.0}};
  gaussian_elimination((double *)a, 3, 4);
  printf("Gaussian elimination:\n");
  for (size_t i = 0; i < 3; ++i) {
    printf("[");
    for (size_t j = 0; j < 4; ++j) {
      printf("%f ", a[i][j]);
    }
    printf("]\n");
  }
  printf("\nGauss-Jordan:\n");
  gauss_jordan((double *)a, 3, 4);
  for (size_t i = 0; i < 3; ++i) {
    printf("[");
    for (size_t j = 0; j < 4; ++j) {
      printf("%f ", a[i][j]);
    }
    printf("]\n");
  }
  printf("\nSolutions are:\n");
  double x[3] = {0, 0, 0};
  back_substitution((double *)a, x, 3, 4);
  printf("(%f,%f,%f)\n", x[0], x[1], x[2]);
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0xb0,%rsp
movl   $0x0,-0x4(%rbp)
lea    -0x70(%rbp),%rdi
lea    0xe83(%rip),%rsi
mov    $0x60,%edx
call   1070 <memcpy@plt>
lea    -0x70(%rbp),%rdi
mov    $0x3,%esi
mov    $0x4,%edx
call   1050 <gaussian_elimination@plt>
lea    0xebf(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
movq   $0x0,-0x78(%rbp)
cmpq   $0x3,-0x78(%rbp)
jae    1232 <main+0xd2>
lea    0xeb5(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
movq   $0x0,-0x80(%rbp)
cmpq   $0x4,-0x80(%rbp)
jae    1213 <main+0xb3>
mov    -0x78(%rbp),%rcx
lea    -0x70(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
mov    -0x80(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
lea    0xe7e(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
mov    -0x80(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x80(%rbp)
jmp    11d1 <main+0x71>
lea    0xe63(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
mov    -0x78(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x78(%rbp)
jmp    11b0 <main+0x50>
lea    0xe47(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    -0x70(%rbp),%rdi
mov    $0x3,%esi
mov    $0x4,%edx
call   1030 <gauss_jordan@plt>
movq   $0x0,-0x88(%rbp)
cmpq   $0x3,-0x88(%rbp)
03
jae    12fb <main+0x19b>
lea    0xe04(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
movq   $0x0,-0x90(%rbp)
cmpq   $0x4,-0x90(%rbp)
04
jae    12d6 <main+0x176>
mov    -0x88(%rbp),%rcx
lea    -0x70(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
mov    -0x90(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
lea    0xdc1(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
mov    -0x90(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x90(%rbp)
jmp    1285 <main+0x125>
lea    0xda0(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
mov    -0x88(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x88(%rbp)
jmp    125e <main+0xfe>
lea    0xd8e(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    -0xb0(%rbp),%rdi
xor    %esi,%esi
mov    $0x18,%edx
call   1060 <memset@plt>
lea    -0x70(%rbp),%rdi
lea    -0xb0(%rbp),%rsi
mov    $0x3,%edx
mov    $0x4,%ecx
call   1080 <back_substitution@plt>
movsd  -0xb0(%rbp),%xmm0
ff
movsd  -0xa8(%rbp),%xmm1
ff
movsd  -0xa0(%rbp),%xmm2
ff
lea    0xd4c(%rip),%rdi
mov    $0x3,%al
call   1040 <printf@plt>
xor    %eax,%eax
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %r15
push   %r14
push   %rbx
sub    $0x80,%rsp
movaps 0xedd(%rip),%xmm0
movaps %xmm0,0x50(%rsp)
movaps 0xec1(%rip),%xmm0
movaps %xmm0,0x40(%rsp)
movaps 0xea5(%rip),%xmm0
movaps %xmm0,0x30(%rsp)
movaps 0xe89(%rip),%xmm0
movaps %xmm0,0x20(%rsp)
movaps 0xe6d(%rip),%xmm0
movaps %xmm0,0x10(%rsp)
movaps 0xe51(%rip),%xmm0
movaps %xmm0,(%rsp)
mov    %rsp,%r14
mov    %r14,%rdi
mov    $0x3,%esi
mov    $0x4,%edx
call   1070 <gaussian_elimination@plt>
lea    0xea1(%rip),%rdi
call   1040 <puts@plt>
mov    $0x5b,%edi
call   1030 <putchar@plt>
movsd  (%rsp),%xmm0
lea    0xe76(%rip),%rbx
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x8(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x10(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x18(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
lea    0xe7a(%rip),%r15
mov    %r15,%rdi
call   1040 <puts@plt>
mov    $0x5b,%edi
call   1030 <putchar@plt>
movsd  0x20(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x28(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x30(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x38(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
mov    %r15,%rdi
call   1040 <puts@plt>
mov    $0x5b,%edi
call   1030 <putchar@plt>
movsd  0x40(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x48(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x50(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x58(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
mov    %r15,%rdi
call   1040 <puts@plt>
lea    0xda8(%rip),%rdi
call   1040 <puts@plt>
mov    %r14,%rdi
mov    $0x3,%esi
mov    $0x4,%edx
call   1050 <gauss_jordan@plt>
mov    $0x5b,%edi
call   1030 <putchar@plt>
movsd  (%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x8(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x10(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x18(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
mov    %r15,%rdi
call   1040 <puts@plt>
mov    $0x5b,%edi
call   1030 <putchar@plt>
movsd  0x20(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x28(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x30(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x38(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
mov    %r15,%rdi
call   1040 <puts@plt>
mov    $0x5b,%edi
call   1030 <putchar@plt>
movsd  0x40(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x48(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x50(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
movsd  0x58(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1060 <printf@plt>
mov    %r15,%rdi
call   1040 <puts@plt>
lea    0xca4(%rip),%rdi
call   1040 <puts@plt>
xorps  %xmm0,%xmm0
movaps %xmm0,0x60(%rsp)
movq   $0x0,0x70(%rsp)
lea    0x60(%rsp),%rsi
mov    %r14,%rdi
mov    $0x3,%edx
mov    $0x4,%ecx
call   1080 <back_substitution@plt>
movsd  0x60(%rsp),%xmm0
movsd  0x68(%rsp),%xmm1
movsd  0x70(%rsp),%xmm2
lea    0xc2d(%rip),%rdi
mov    $0x3,%al
call   1060 <printf@plt>
xor    %eax,%eax
add    $0x80,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
@__const.main.a = [3 x [4 x double]] [[4 x double] [double 3.000000e+00, double 2.000000e+00, double -4.000000e+00, double 3.000000e+00], [4 x double] [double 2.000000e+00, double 3.000000e+00, double 3.000000e+00, double 1.500000e+01], [4 x double] [double 5.000000e+00, double -3.000000e+00, double 1.000000e+00, double 1.400000e+01]]
@.str = [23 x i8] c"Gaussian elimination:\0A\00"
@.str.1 = [2 x i8] c"[\00"
@.str.2 = [4 x i8] c"%f \00"
@.str.3 = [3 x i8] c"]\0A\00"
@.str.4 = [16 x i8] c"\0AGauss-Jordan:\0A\00"
@.str.5 = [17 x i8] c"\0ASolutions are:\0A\00"
@.str.6 = [12 x i8] c"(%f,%f,%f)\0A\00"
define i32 @main() {
%1 = alloca i32
%2 = alloca [3 x [4 x double]]
%3 = alloca i64
%4 = alloca i64
%5 = alloca i64
%6 = alloca i64
%7 = alloca [3 x double]
store i32 0, i32* %1
%8 = bitcast [3 x [4 x double]]* %2 to i8*
%9 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %2, i64 0, i64 0
%10 = bitcast [4 x double]* %9 to double*
%11 = call i32 @gaussian_elimination(double* %10, i32 3, i32 4)
%12 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([23 x i8], [23 x i8]* @.str, i64 0, i64 0))
store i64 0, i64* %3
br label %13
13: ; preds = %33, %0
%14 = load i64, i64* %3
%15 = icmp ult i64 %14, 3
br i1 %15, label %16, label %36
16: ; preds = %13
%17 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.1, i64 0, i64 0))
store i64 0, i64* %4
br label %18
18: ; preds = %28, %16
%19 = load i64, i64* %4
%20 = icmp ult i64 %19, 4
br i1 %20, label %21, label %31
21: ; preds = %18
%22 = load i64, i64* %3
%23 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %2, i64 0, i64 %22
%24 = load i64, i64* %4
%25 = getelementptr inbounds [4 x double], [4 x double]* %23, i64 0, i64 %24
%26 = load double, double* %25
%27 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %26)
br label %28
28: ; preds = %21
%29 = load i64, i64* %4
%30 = add i64 %29, 1
store i64 %30, i64* %4
br label %18
31: ; preds = %18
%32 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.3, i64 0, i64 0))
br label %33
33: ; preds = %31
%34 = load i64, i64* %3
%35 = add i64 %34, 1
store i64 %35, i64* %3
br label %13
36: ; preds = %13
%37 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([16 x i8], [16 x i8]* @.str.4, i64 0, i64 0))
%38 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %2, i64 0, i64 0
%39 = bitcast [4 x double]* %38 to double*
%40 = call i32 @gauss_jordan(double* %39, i32 3, i32 4)
store i64 0, i64* %5
br label %41
41: ; preds = %61, %36
%42 = load i64, i64* %5
%43 = icmp ult i64 %42, 3
br i1 %43, label %44, label %64
44: ; preds = %41
%45 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.1, i64 0, i64 0))
store i64 0, i64* %6
br label %46
46: ; preds = %56, %44
%47 = load i64, i64* %6
%48 = icmp ult i64 %47, 4
br i1 %48, label %49, label %59
49: ; preds = %46
%50 = load i64, i64* %5
%51 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %2, i64 0, i64 %50
%52 = load i64, i64* %6
%53 = getelementptr inbounds [4 x double], [4 x double]* %51, i64 0, i64 %52
%54 = load double, double* %53
%55 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %54)
br label %56
56: ; preds = %49
%57 = load i64, i64* %6
%58 = add i64 %57, 1
store i64 %58, i64* %6
br label %46
59: ; preds = %46
%60 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.3, i64 0, i64 0))
br label %61
61: ; preds = %59
%62 = load i64, i64* %5
%63 = add i64 %62, 1
store i64 %63, i64* %5
br label %41
64: ; preds = %41
%65 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([17 x i8], [17 x i8]* @.str.5, i64 0, i64 0))
%66 = bitcast [3 x double]* %7 to i8*
%67 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %2, i64 0, i64 0
%68 = bitcast [4 x double]* %67 to double*
%69 = getelementptr inbounds [3 x double], [3 x double]* %7, i64 0, i64 0
%70 = call i32 @back_substitution(double* %68, double* %69, i32 3, i32 4)
%71 = getelementptr inbounds [3 x double], [3 x double]* %7, i64 0, i64 0
%72 = load double, double* %71
%73 = getelementptr inbounds [3 x double], [3 x double]* %7, i64 0, i64 1
%74 = load double, double* %73
%75 = getelementptr inbounds [3 x double], [3 x double]* %7, i64 0, i64 2
%76 = load double, double* %75
%77 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([12 x i8], [12 x i8]* @.str.6, i64 0, i64 0), double %72, double %74, double %76)
ret i32 0
}

LLVM IR (-O2):
@__const.main.a = [3 x [4 x double]] [[4 x double] [double 3.000000e+00, double 2.000000e+00, double -4.000000e+00, double 3.000000e+00], [4 x double] [double 2.000000e+00, double 3.000000e+00, double 3.000000e+00, double 1.500000e+01], [4 x double] [double 5.000000e+00, double -3.000000e+00, double 1.000000e+00, double 1.400000e+01]]
@.str.2 = [4 x i8] c"%f \00"
@.str.6 = [12 x i8] c"(%f,%f,%f)\0A\00"
@str = [22 x i8] c"Gaussian elimination:\00"
@str.7 = [15 x i8] c"\0AGauss-Jordan:\00"
@str.8 = [16 x i8] c"\0ASolutions are:\00"
@str.10 = [2 x i8] c"]\00"
define i32 @main() {
%1 = alloca [3 x [4 x double]]
%2 = alloca [3 x double]
%3 = bitcast [3 x [4 x double]]* %1 to i8*
%4 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %1, i64 0, i64 0, i64 0
%5 = call i32 @gaussian_elimination(double* %4, i32 3, i32 4)
%6 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([22 x i8], [22 x i8]* @str, i64 0, i64 0))
%7 = call i32 @putchar(i32 91)
%8 = load double, double* %4
%9 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %8)
%10 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %1, i64 0, i64 0, i64 1
%11 = load double, double* %10
%12 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %11)
%13 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %1, i64 0, i64 0, i64 2
%14 = load double, double* %13
%15 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %14)
%16 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %1, i64 0, i64 0, i64 3
%17 = load double, double* %16
%18 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %17)
%19 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([2 x i8], [2 x i8]* @str.10, i64 0, i64 0))
%20 = call i32 @putchar(i32 91)
%21 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %1, i64 0, i64 1, i64 0
%22 = load double, double* %21
%23 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %22)
%24 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %1, i64 0, i64 1, i64 1
%25 = load double, double* %24
%26 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %25)
%27 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %1, i64 0, i64 1, i64 2
%28 = load double, double* %27
%29 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %28)
%30 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %1, i64 0, i64 1, i64 3
%31 = load double, double* %30
%32 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %31)
%33 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([2 x i8], [2 x i8]* @str.10, i64 0, i64 0))
%34 = call i32 @putchar(i32 91)
%35 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %1, i64 0, i64 2, i64 0
%36 = load double, double* %35
%37 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %36)
%38 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %1, i64 0, i64 2, i64 1
%39 = load double, double* %38
%40 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %39)
%41 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %1, i64 0, i64 2, i64 2
%42 = load double, double* %41
%43 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %42)
%44 = getelementptr inbounds [3 x [4 x double]], [3 x [4 x double]]* %1, i64 0, i64 2, i64 3
%45 = load double, double* %44
%46 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %45)
%47 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([2 x i8], [2 x i8]* @str.10, i64 0, i64 0))
%48 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([15 x i8], [15 x i8]* @str.7, i64 0, i64 0))
%49 = call i32 @gauss_jordan(double* %4, i32 3, i32 4)
%50 = call i32 @putchar(i32 91)
%51 = load double, double* %4
%52 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %51)
%53 = load double, double* %10
%54 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %53)
%55 = load double, double* %13
%56 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %55)
%57 = load double, double* %16
%58 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %57)
%59 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([2 x i8], [2 x i8]* @str.10, i64 0, i64 0))
%60 = call i32 @putchar(i32 91)
%61 = load double, double* %21
%62 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %61)
%63 = load double, double* %24
%64 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %63)
%65 = load double, double* %27
%66 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %65)
%67 = load double, double* %30
%68 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %67)
%69 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([2 x i8], [2 x i8]* @str.10, i64 0, i64 0))
%70 = call i32 @putchar(i32 91)
%71 = load double, double* %35
%72 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %71)
%73 = load double, double* %38
%74 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %73)
%75 = load double, double* %41
%76 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %75)
%77 = load double, double* %44
%78 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %77)
%79 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([2 x i8], [2 x i8]* @str.10, i64 0, i64 0))
%80 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([16 x i8], [16 x i8]* @str.8, i64 0, i64 0))
%81 = bitcast [3 x double]* %2 to i8*
%82 = getelementptr inbounds [3 x double], [3 x double]* %2, i64 0, i64 0
%83 = call i32 @back_substitution(double* %4, double* %82, i32 3, i32 4)
%84 = load double, double* %82
%85 = getelementptr inbounds [3 x double], [3 x double]* %2, i64 0, i64 1
%86 = load double, double* %85
%87 = getelementptr inbounds [3 x double], [3 x double]* %2, i64 0, i64 2
%88 = load double, double* %87
%89 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([12 x i8], [12 x i8]* @.str.6, i64 0, i64 0), double %84, double %86, double %88)
ret i32 0
}