#!/usr/bin/env python3
# -*- coding: utf-8 -*-

import numpy as np

# ---------------------------------
# Lecture 11 - Gauss Elimination
# (C) Solve the matrix using 
#     Gauss elimination
# ---------------------------------
#

A = np.array([[ 5.,  5.,  5., -6.,  1.],
              [-2.,  2.,  4.,  4.,  5.],
              [-2.,  5., -4.,  0.,  0.],
              [ 1., -4.,  0.,  4., -4.],
              [-2., -2.,  1., -2., -2.]])

b = np.array([ 7., -8., 46., -49., 4.])

# Forward Elimination
       
# Back Substitution