def PrintSourceMatrix(S,params):
    """Auxiliary function to print the source matrix (S)"""

    dim = params["dim"]	# Spatial dimensions
    print("The source term matrix is:\n")
    for i in range (0,dim+2):
        for j in range (0,dim+2):
            print("S[",i,",",j,"]=",S[i,j],"\n")

    return 0