# -*- coding: utf-8 -*-
"""
Created on Tue Oct 11 19:16:37 2022

@author: Roman
"""

import numpy as np

a = np.array([1, 2, 3])

b = np.array([[1, 2], [3, 3]])

z = np.zeros((2, 3))

o = np.ones((2, 2))


print(a)
print(b)
print(z)
print(o)

print(o + b)
print(o * b)
print(o @ b)
print(np.dot(o, b))
print('transp', b.T)

print(np.linalg.det(b))
print(np.linalg.norm(b))
print(np.linalg.inv(b))
print(np.trace(b))


a = np.array([[1, 2], [3, 5]])
b = np.array([1, 2])
print(np.linalg.solve(a, b))

x = np.array([1, 2, 3, 4, 5])
y = np.array([1, 3, 7, 10, 13])

p = np.polyfit(x, y, 1)
print(p)

x1 = np.polyval(p, 2.5)
print(x1)

x2 = np.polyval(p, x)
print(x2)


b = np.arange(0, 1 + 0.1, 0.25)
print("b =\n", b)

у = np.eye(3)
print("e =\n", y)