python numpy作业

代码：

import numpy
import time
from scipy import linalg

n = 200
m = 500
A = numpy.random.normal(size = (n, m))
f = numpy.random.randint(1, 1000, size = m)
B = linalg.toeplitz(f)

#9.1
#A+A
A_A = A + A
'''print(A_A)'''
#AAT
AAT = numpy.dot(A, A.T)
'''print(AAT)'''
#ATA
ATA = numpy.dot(A.T, A)
'''print(ATA)'''
#AB
AB = numpy.dot(A, B)
'''print(AB)'''
#计算A(B-sI)
def func1(A, B, s):
    I = numpy.eye(B.s