Implement a python program that uses the scipy library to solve a system of linear equations using Gaussian elimination. The input should be a matrix of coefficients and a vector of constants, and the output should be the solution vector. PLAINFORMAT
import numpy as np
from scipy.linalg import solve

def gaussian_elimination(coeff_matrix, const_vector):
    solution = solve(coeff_matrix, const_vector)
    return solution