matrix using numpy python

#This code is valid for only 3x3 matrix

import numpy as np

R1 = int(input("Enter the number of rows:"))
C1 = int(input("Enter the number of columns:"))

print("Enter the entries in a single line (separated by space): ")

# User input of entries in a
# single line separated by space
entries = list(map(int, input().split()))

# For printing the matrix
matrix1 = np.array(entries).reshape(R1, C1)
print(matrix1)

R2 = int(input("Enter the number of rows:"))
C2 = int(input("Enter the number of columns:"))
print("Enter the entries in a single line (separated by space): ")
entries = list(map(int, input().split()))
matrix2 = np.array(entries).reshape(R2, C2)
print(matrix2)
result = [[0,0,0],
         [0,0,0],
         [0,0,0]]

#adding
for i in range(len(matrix1)):
   # iterate through columns
   for j in range(len(matrix1[0])):
       result[i][j] = matrix1[i][j] + matrix2[i][j]

for r in result:
   print(r)

# Subtracting
for i in range(len(matrix1)):
   # iterate through columns
   for j in range(len(matrix1[0])):
       result[i][j] = matrix1[i][j] - matrix2[i][j]

for a in result:
   print(a)


# multiplication
for i in range(len(matrix1)):
   # iterate through columns
   for j in range(len(matrix1[0])):
       result[i][j] = matrix1[i][j] * matrix2[i][j]

for m in result:
   print(m)

#division
for i in range(len(matrix1)):
   # iterate through columns
   for j in range(len(matrix1[0])):
       result[i][j] = matrix1[i][j] / matrix2[i][j]

for d in result:
   print(d)



#inverse of matrix1

print(np.linalg.inv(matrix1))


#transpose of matrix1
matrix1_transpose = np.transpose(matrix1)
print(matrix1_transpose)