import numpy as np
from scipy.linalg import lstsq
from scipy.linalg import norm
A = np.random.randint(0,50,size=(5,5))
b = np.random.randint(0,50,size=(5,1))
print("A:")
print(A)
print("\nb:")
print(b)
x, res, rnk, s = lstsq(A, b)
y = norm(b-np.dot(A,x),ord=2)
print("\nSolution x for Ax=b:")
print(x)
print("\ny:")
print(y)
import scipy.optimize as op
import numpy as np
import math
def fun(x):
return (-(math.sin(x-2)**2)*math.exp(-(x**2)))
Max = op.minimize_scalar(fun)
print(-Max.fun)
from scipy.spatial.distance import cdist
import numpy as np
import math
X = np.random.randint(0,50,(5,5))
print(X)
print("\nSolution:")
print(cdist(X,X,metric="euclidean"))