import numpy as np
import time


def main():
    
    n=input("Ordine della matrice :\n")
    n=int(n)

    ma=np.zeros((n,n),dtype=np.float64)


    for i in range(n):
        for j in range (n):
            ma[i,j]=np.sin(j)


    mb=np.zeros((n,n),dtype=np.float64)

    for i in range(n):
        for j in range (n):
            mb[i,j]=np.cos(j)
       
    mc=np.zeros((n,n),dtype=np.float64)

    print("Ora moltiplico")

    start_time = time.perf_counter()
    for i in range(n):
        for j in range (n):
            for k in range (n):
                mc[i,j]+=ma[i,k]*mb[k,j]

    end_time = time.perf_counter()

    elapsed_time = end_time - start_time

    print(str(mc[0,0])+ " "+str(mc[1,0])+"  "+ str(mc[2,0]))
    print("Elapsed time: ", elapsed_time)   

if __name__ == "__main__":
    main()