import math
import numpy as np
from scipy.linalg import solve

print("Math Library:")
print("Square root of 16:", math.sqrt(16))
print("Factorial of 5:", math.factorial(5))
print("Value of Pi:", math.pi)

print("\nNumPy Library:")
arr = np.array([1, 2, 3, 4, 5])
print("Array:", arr)
print("Mean of array:", np.mean(arr))
print("Standard deviation of array:", np.std(arr))

print("\nSciPy Library:")
# Solve the system of equations:
#  2x + y = 10
#  x + 3y = 13
A = np.array([[2, 1], [1, 3]])  # Coefficients
B = np.array([10, 13])  # Constants
solution = solve(A, B)  # Solve Ax = B
print("Solution for x and y:", solution)