Write the python code for the following

import numpy as np from copy import deepcopy mat1 np . array ( [ [ 1 , 2, 4 ] , [ 4, 5, 6] , [2,2,2],[3,4,5]],dtype-np . float64) = def QRrot (A): #Replace this comment with code -write a function to compute the QR decomposition of an nxm matrix A using rotators -assume A is full rank and n>m -The function should work for any size matrix -The function should work for matrices with real entries - use a tolerance of 1e-15 when checking to see if a value is o return(Q,R