2021년 2월 23일 화요일

CalculusMadeEasy_CH1

1-6

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

b=np.array([[-1],[-3],[-4]])
la.norm(b)
5.0990195135927845

c=np.array([[1],[-2],[4]])
la.norm(c)
4.58257569495584

1-7

a=np.array([[1],[2],[-3]])
b=np.array([[-4],[-3],[6]])
ab_dot=np.dot(a.T, b)
ab_dot
array([[-28]])
cos_ab=ab_dot/(la.norm(a)*la.norm(b))
cos_ab
array([[-0.95814028]])
rad=np.arccos(cos_ab)
rad
array([[2.85123057]])
theta=np.rad2deg(rad)
theta
array([[163.36347829]])

1-8

t=symbols("t")
A=Matrix([[2,2],[-3,-3],[4,4]])
B=Matrix([[2,1],[4,-1],[5,-1]])
x=Matrix([[1],[t]])
A*x
$\left[\begin{matrix}2 t + 2\\- 3 t - 3\\4 t + 4\end{matrix}\right]$
B*x
$\left[\begin{matrix}t + 2\\4 - t\\5 - t\end{matrix}\right]$
eq=(A*x).T*(B*x)
eq
[(4−𝑡)(−3𝑡−3)+(5−𝑡)(4𝑡+4)+(𝑡+2)(2𝑡+2)]
eq=expand(eq)
eq
[𝑡2+13𝑡+12]
solve(eq)
[{t: -12}, {t: -1}]

선형대수코드_CH2

2-4

다음 행렬 A, B에 대해 AB=BA가 성립하도록 a, b, c, d를 계산하세요.

a,b,c,d=symbols("a b c d")
A=Matrix([[a, b],[c,d]])
B=Matrix([[1, 2], [3, 4]])
ab=A*B
ab
$\left[\begin{matrix}a + 3 b & 2 a + 4 b\\c + 3 d & 2 c + 4 d\end{matrix}\right]$
ba=B*A
ba
$\left[\begin{matrix}a + 2 c & b + 2 d\\3 a + 4 c & 3 b + 4 d\end{matrix}\right]$
Eq(ab,ba)
$\left[\begin{matrix}a + 3 b & 2 a + 4 b\\c + 3 d & 2 c + 4 d\end{matrix}\right] = \left[\begin{matrix}a + 2 c & b + 2 d\\3 a + 4 c & 3 b + 4 d\end{matrix}\right]$
solve(Eq(ab,ba))
{a: -c + d, b: 2*c/3}


2-11

2$\times$2 정방행렬이 가역행렬이기 위해서는 다른 동일차원의 행렬과의 곱으로 항등행렬이 존제해야 합니다. 즉, 다음이 성립해야 합니다.

$\left[\begin{matrix}a&b\\c&d\end{matrix}\right] \left[\begin{matrix}x&y\\z&w\end{matrix}\right]=\left[\begin{matrix}1&0\\0&1\end{matrix}\right]$

위 선형 결합의 선형독립이 성립하기 위해서는 확대행렬의 rref에서 해의 존재를 보여야 합니다.

a,b,c,d,x,y,z,w=symbols("a,b,c,d,x,y,z,w")
A=Matrix([[a, b],[c,d]])
B=Matrix([[x, y], [z, w]])
I2=eye(2)
A

$\left[\begin{matrix}a & b\\c & d\end{matrix}\right]$

I2
$\left[\begin{matrix}1 & 0\\0 & 1\end{matrix}\right]$

au=A.row_join(I2)
au
$\left[\begin{matrix}a & b & 1 & 0\\c & d & 0 & 1\end{matrix}\right]$

au.rref()
(Matrix([[1, 0, d/(a*d - b*c), -b/(a*d - b*c)],
        [0, 1, -c/(a*d - b*c), a/(a*d - b*c)]]),
        (0, 1))

위 결과는 다음과 같이 정리되므로 역행렬이 존재하기 위해서는 $ad-bc \neq 0$ 이어야 합니다.
$\begin{align}x&=\frac{d}{(ad - bc)}\\w&=\frac{a}{(ad - bc)} \end{align}$

2_12

다음 선형시스템의 해?

2x1 + 5x2 − 8x3 + 2x4 + 2x5  =  0
6x1 + 2x2 − 10x3 + 6x4 + 8x5  =  6
3x1 + 6x2 + 2x3 + 3x4 + 5x5  =  6
3x1 + x2 − 5x3 + 2x4 + 2x5  =  3
6x1 + 7x2 − 3x3 + 6x4 + 9x5  =  9
A=Matrix([[2,5,-8,2,2],[6,2,-10,6,8],[3,6,2,3,5],[3,1,-5,2,2],[6,7,-3,6,9]])
C=Matrix(5, 1, [0,6,6,3,9])
A
$\left[\begin{matrix}2 & 5 & -8 & 2 & 2\\6 & 2 & -10 & 6 & 8\\3 & 6 & 2 & 3 & 5\\3 & 1 & -5 & 2 & 2\\6 & 7 & -3 & 6 & 9\end{matrix}\right]$

C
$\left[\begin{matrix}0\\6\\6\\3\\9\end{matrix}\right]$

au=A.row_join(C)
au
$\left[\begin{matrix}2 & 5 & -8 & 2 & 2 & 0\\6 & 2 & -10 & 6 & 8 & 6\\3 & 6 & 2 & 3 & 5 & 6\\3 & 1 & -5 & 2 & 2 & 3\\6 & 7 & -3 & 6 & 9 & 9\end{matrix}\right]$

au.rref()
(Matrix([[1, 0, 0, 0, -3/7, 12/7],
        [0, 1, 0, 0, 0, 0],
        [0, 0, 1, 0, 1/7, 3/7],
        [0, 0, 0, 1, 2, 0],
        [0, 0, 0, 0, 0, 0]]),
        (0, 1, 2, 3))

위 결과를 정리하면 다음과 같습니다.

$\begin{align}x_1&=\frac{3}{7}x_5+\frac{12}{7}\\ x_2&=0\\x_3&=-\frac{1}{7}x_5+\frac{3}{7}\\x_4&=-2x_5\end{align}$

선형대수코드_CH3

3-1

(a)

A=Matrix([[1,2,4,4,-7], [0,1,1,2,-3],[1,0,2,0,-1]])
A.rref()
(Matrix([[1, 0, 2, 0, -1],
        [0, 1, 1, 2, -3],
        [0, 0, 0, 0, 0]]),
        (0, 1))

(b) 이 rref로부터 선도변수가 2개, 3개의 자유변수가 존재함을 알 수 있습니다. 그러므로 위 결과를 정리하면 다음과 같습니다. 즉, x1, x2는 자유변수 x3, x4에 따라 다양한 해를 가집니다. 즉, 자명하지 않은 해를 가지므로 선형 종속입니다.

$\begin{align}x_1&=-2x_3+1\\x_2&=-x_3-2x_4+3\end{align}$

(c)위 시스템이 자명한 해를 갖는 선형독립이기 위해서는 자유변수로 이루어진 방정식이 선형독립인 경우입니다. 다음과 같이 이 식은 선형독립이며 $x_3=-\frac{1}{2}$와 $x_4=-\frac{4}{5}$인 경우 선형독립이 됩니다.

B=Matrix([[-2, 0, -1], [-1, -2, -3]])
B.rref()
(Matrix([[1, 0, 1/2],
        [0, 1, 5/4]]),
        (0, 1))

3-2

x1,x2,x3,x4=symbols('x:4')
A=Matrix([[1,1,2,-1],[1,3,-1,2],[1,1,1,3],[1,2,1,1]])
x=Matrix([[x1],[x2],[x3],[x4]])
A*x
$\left[\begin{matrix}x_{0} + x_{1} + 2 x_{2} - x_{3}\\x_{0} + 3 x_{1} - x_{2} + 2 x_{3}\\x_{0} + x_{1} + x_{2} + 3 x_{3}\\x_{0} + 2 x_{1} + x_{2} + x_{3}\end{matrix}\right]$
b=Matrix([[0],[0],[0],[0]])
Au_a=A.row_join(b)
Au_a
$ \left[\begin{matrix}1 & 1 & 2 & -1 & 0\\1 & 3 & -1 & 2 & 0\\1 & 1 & 1 & 3 & 0\\1 & 2 & 1 & 1 & 0\end{matrix}\right]$
Au_a.rref()
(Matrix([[1, 0, 0, 0, 0],
        [0, 1, 0, 0, 0],
        [0, 0, 1, 0, 0],
        [0, 0, 0, 1, 0]]),
        (0, 1, 2, 3))


B=Matrix([[1,2,3],[0,1,0],[1,0,3]])
x=Matrix([[x1],[x2],[x3]])
B*x
$\left[\begin{matrix}x_{0} + 2 x_{1} + 3 x_{2}\\x_{1}\\x_{0} + 3 x_{2}\end{matrix}\right]$
c=Matrix([[0],[0],[0]])
Au_b=B.row_join(c)
Au_b
$\left[\begin{matrix}1 & 2 & 3 & 0\\0 & 1 & 0 & 0\\1 & 0 & 3 & 0\end{matrix}\right]$
Au_b.rref()
(Matrix([[1, 0, 3, 0],
        [0, 1, 0, 0],
        [0, 0, 0, 0]]),
        (0, 1))


C=Matrix([[1,2,-1],[2,-1,3],[5,-4,3]])
x=Matrix([[x1],[x2],[x3]])
C*x
$\left[\begin{matrix}x_{0} + 2 x_{1} - x_{2}\\2 x_{0} - x_{1} + 3 x_{2}\\5 x_{0} - 4 x_{1} + 3 x_{2}\end{matrix}\right]$
c=Matrix([[0],[0],[0]])
Au_c=C.row_join(c)
Au_c
$\left[\begin{matrix}1 & 2 & -1 & 0\\2 & -1 & 3 & 0\\5 & -4 & 3 & 0\end{matrix}\right]$
Au_c.rref()
(Matrix([[1, 0, 0, 0],
        [0, 1, 0, 0],
        [0, 0, 1, 0]]),
        (0, 1, 2))

3-7

a=symbols('a')
A=Matrix([[1,1,1],[2,3,2],[2,3,a**2-5]])
b=Matrix([[2],[3],[a]])
solve(A.det())
[-sqrt(7), sqrt(7)]

aug=A.row_join(b)
print(aug)
Matrix([[1, 1, 1, 2],
        [2, 3, 2, 3],
        [2, 3, a**2 - 5, a]])


aug.rref()
(Matrix([[1, 0, 0, (3*a**2 - a - 18)/(a**2 - 7)],
        [0, 1, 0, -1],
        [0, 0, 1, (a - 3)/(a**2 - 7)]]),
        (0, 1, 2))


solve(3*a**2 - a - 18)
[1/6 - sqrt(217)/6, 1/6 + sqrt(217)/6]

3-9

다음 선형시스템의 해의 수를 결정
$\begin{align}x+2y-3z&=4\\4x+y+2z&=6\\x+2y+(a^2-19)z&=a\end{align}$

a=symbols("a", real=True)
A=Matrix([[1, 2, -3],[4, 1, 2],[1, 2, a**2-19]])
C=Matrix(3,1,[4, 6, a])
A
$\left[\begin{matrix}1 & 2 & -3\\4 & 1 & 2\\1 & 2 & a^{2} - 19\end{matrix}\right]$

au=A.row_join(C)
au
$\left[\begin{matrix}1 & 2 & -3\\4 & 1 & 2\\1 & 2 & a^{2} - 19\end{matrix}\right]$

au.rref()[0]
$\left[\begin{matrix}1 & 0 & 0 & \frac{k \left(- k \left(k + 2\right) - 4 k - 5\right) - \left(k + 2\right) \left(k \left(k + 2\right) - k \left(4 k + 4\right) + 1\right)}{- k \left(k + 2\right) - 4 k - 5}\\0 & 1 & 0 & \frac{- k \left(k + 2\right) + k \left(4 k + 4\right) + k \left(- k \left(k + 2\right) - 4 k - 5\right) - 1}{- k \left(k + 2\right) - 4 k - 5}\\0 & 0 & 1 & \frac{k \left(k + 2\right) - k \left(4 k + 4\right) + 1}{- k \left(k + 2\right) - 4 k - 5}\end{matrix}\right]$

solve(112 - 7*a**2, a)
[-4, 4]
sol=au.rref()[0][:,3]
sol
$\left[\begin{matrix}\frac{- 8 a^{2} + 7 a + 100}{112 - 7 a^{2}}\\\frac{- 10 a^{2} - 14 a + 216}{112 - 7 a^{2}}\\\frac{a - 4}{a^{2} - 16}\end{matrix}\right]$

sol.subs(a, 4)
$\left[\begin{matrix}\text{NaN}\\\text{NaN}\\\text{NaN}\end{matrix}\right]$

sol.subs(a, -4)
$\left[\begin{matrix}\tilde{\infty}\\\tilde{\infty}\\\tilde{\infty}\end{matrix}\right]$
  • a = 4 인 경우 해가 없음
  • a = -4 인 경우 무한히 많은 해가 존재
  • $a \neq \pm 4$인 경우 단일한 해 존재

3_10

u=Matrix([[5],[10]])
v=Matrix([[3], [11]])
w=Matrix([[-1],[-12]])
au=(u.row_join(v)).row_join(w)
au
$\left[\begin{matrix}5 & 3 & -1\\10 & 11 & -12\end{matrix}\right]$
latex(au.rref())


$\left( \left[\begin{matrix}1 & 0 & 1\\0 & 1 & -2\end{matrix}\right], \ \left( 0, \ 1\right)\right)$

벡터 u, v, w의 선형결합은 선형독립입니다. 즉, u, v는 w의 기저.
a=False
b=True

선형대수코드_CH6

1

 1. 행렬 A에 대해
(a) 0벡터에 대응하는 모든 정의역의 벡터?

A=np.array([[3,-2, 1, -1], [-3, 5, 2, 4],[-9, 12, 3, 11],[3, -8, -5, -7]])
print(A)
[[ 3 -2 1 -1]
   [-3 5 2 4]
   [-9 12 3 11]
   [ 3 -8 -5 -7]]


b=np.array([[0],[0],[0], [0]])
print(b)
[[0]
   [0]
   [0]
   [0]]


au=np.hstack([A, b])
print(au)
[[ 3 -2 1 -1 0]
   [-3 5 2 4 0]
   [-9 12 3 11 0]
   [ 3 -8 -5 -7 0]]


Matrix(au).rref()
(Matrix([[1, 0, 1, 0, 0],
        [0, 1, 1, 0, 0],
        [0, 0, 0, 1, 0],
        [0, 0, 0, 0, 0]]),
        (0, 1, 3))

Matrix(A).nullspace()
[Matrix([[-1],
        [-1],
        [ 1],
        [ 0]])]

위 코드에서 나타낸 것과 같이 Ax=b의 확대행렬 au의 기약행사다리꼴로부터 선형종속임을 알 수 있습니다. 그러므로 이 시스템의 해집합의 기저에 대한 스칼라배한 결과가 모두 x벡터가 될 수 있습니다.

$scalar \cdot \left[\begin{matrix}-1\\-1\\1\\0\end{matrix}\right]$

(b) 다음 벡터 b가 치역에 존재합니까?
$b = \left[\begin{matrix}-2\\14\\36\\-27\end{matrix}\right]$

b=np.array([[-2],[14],[36],[-27]])
print(b)
[[ -2]
   [ 14]
   [ 36]
   [-27]]


Matrix(np.hstack([A, b])).rref()
(Matrix([[1, 0, 1, 0, 0],
        [0, 1, 1, 0, 0],
        [0, 0, 0, 1, 0],
        [0, 0, 0, 0, 1]]),
        (0, 1, 3, 4))

위 결과에 의하면 4행 즉, $0x_1+0x_2+0x_3+0x_4=1$은 성립하지 않습니다. 그러므로 T(x)=b는 존재할 수 없습니다.

2

T: ℝ2  → ℝ2에서 벡터 u, T(u), v, T(v)는 다음과 같습니다.
(a) T(2u)?
 (b)변환 T의 표준행렬인 A?

(a) 선형변환의 경우 다음 두 조건이 만족합니다.
T(a+b)=T(a)+T(b)
T(ca)=cT(a), c: 스칼라
그러므로 T(2u)=2T(u)

u=Matrix([[-2], [2]])
v=Matrix([[-1],[-3]])
Tu=Matrix(2,1,[6,4])
Tv=Matrix(2,1, [7, 6])
2*Tu
$\left[\begin{matrix}12\\8\end{matrix}\right]$

(b) 표준행렬 A의 경우 다음 코드와 같이 sympy 모듈의 symbols() 함수를 사용하여 생성할 수 있습니다. 이 행렬을 사용하여 Au=Tu, Av=Tv의 식들을 계산하여 미지수 4개의 해를 결정합니다. 이 과정에서 solve() 함수를 적용합니다.

x1,x2,x3,x4=symbols('x1 x2 x3 x4', real=True)
A=Matrix([[x1, x2],[x3, x4]])
A
$\left[\begin{matrix}x_{1} & x_{2}\\x_{3} & x_{4}\end{matrix}\right]$

eq1=Eq(A*u, Tu)
eq1
$\left[\begin{matrix}- 2 x_{1} + 2 x_{2}\\- 2 x_{3} + 2 x_{4}\end{matrix}\right] = \left[\begin{matrix}6\\4\end{matrix}\right]$

eq2=Eq(A*v, Tv)
eq2
$ \left[\begin{matrix}- x_{1} - 3 x_{2}\\- x_{3} - 3 x_{4}\end{matrix}\right] = \left[\begin{matrix}7\\6\end{matrix}\right]$

sol=solve((eq1, eq2), (x1, x2, x3, x4))
sol
{x1: -4, x2: -1, x3: -3, x4: -1}

A1=A.subs(sol)
A1
$\left[\begin{matrix}-4 & -1\\-3 & -1\end{matrix}\right]$

#위 (a)를 확인
A1*2*u
$\left[\begin{matrix}12\\8\end{matrix}\right]$

3

다음 행렬 A의 대각화가 가능합니까?

다음 코드에서 각각의 결과 객체는 latex형태로 나타내기 위해 sympy의 Matrix() 객체로 전환하였습니다.

A=np.array([[6,-2, 0],[0, 6, 2],[0,0,6]])
Matrix(A)
$\left[\begin{matrix}6 & -2 & 0\\0 & 6 & 2\\0 & 0 & 6\end{matrix}\right]$
val, vec=la.eig(A)
D=np.diag(val)
Matrix(D)
$\left[\begin{matrix}6.0 & 0.0 & 0.0\\0.0 & 6.0 & 0.0\\0.0 & 0.0 & 6.0\end{matrix}\right]$

P=vec
Matrix(np.around(P, 10))
$\left[\begin{matrix}1.0 & 1.0 & -1.0\\0.0 & 0.0 & 0.0\\0.0 & 0.0 & 0.0\end{matrix}\right]$

Matrix(np.dot(np.dot(P, D), la.inv(P)))
$\left[\begin{matrix}6.0 & 0.0 & 0.0\\0.0 & 6.0 & 0.5\\0.0 & 0.0 & 6.0\end{matrix}\right]$

위 결과는 행렬 A와 같지 않습니다. 즉, 대각화는 가능하지 않습니다. 유사변환에 의한 분해에서 P는 가역행렬이어야 하지만 생성된 행렬 P의 경우 비가역행렬입니다.

np.around(la.det(P), 10)
0.0

4

다음 벡터들의 집합 {u1, u2, u3}는 직교집합입니까? 그렇다면 이 벡터들의 배수들로 이루어진 직교행렬을 나타내십시요?
$u_1=\left[\begin{matrix}2\\1\\-2\end{matrix}\right]\quad u_2=\left[\begin{matrix}1\\2\\2\end{matrix}\right]\quad u_3=\left[\begin{matrix}2\\-1\\1\end{matrix}\right]$
각 벡터가 직교 관계이면 내적이 0이므로 이들의 계산을 위해 각 벡터를 결합하여 행렬로 전환한 후 계산합니다. 또한 직교관계인 벡터들의 선형결합은 선형독립입니다.

x=np.array([[2,1,2],[1,2,-2],[-2,2,1]])
for i in range(2):
for j in range(i+1, 3):
print(f'{i} and {j}: {np.dot(x[:,i], x[:,j].T)}')
0 and 1: 0
   0 and 2: 0
   1 and 2: 0

Matrix(x).rref()
(Matrix([[1, 0, 0],
      [0, 1, 0],
      [0, 0, 1]]),
      (0, 1, 2))

위 결과에 의하면 위 벡터 집합은 직교집합입니다. 이들 벡터의 배수로 형성된 직교행렬은 각 벡터를 단위 행렬로 전환한 후 결합한 행렬로 나타낼 수 있습니다.

x=np.array([[2,1,2],[1,2,-2],[-2,2,1]], dtype=np.float64)
print(x)
[[ 2. 1. 2.]
   [ 1. 2. -2.]
   [-2. 2. 1.]]


x_n=la.norm(x, axis=0)
x_n
array([3., 3., 3.])
for i in range(3):
 x[:,i]=x[:, i]/x_n[i]
print(np.around(x, 3))
[[ 0.222 0.111 0.222]
   [ 0.111 0.222 -0.222]
   [-0.222 0.222 0.111]]

5

다음 행렬 A는 직교행렬입니까?

U=np.array([[2,5],[-5, 2]])
print(U)
[[ 2 5]
   [-5 2]]


print(U.T)
[[ 2 -5]
   [ 5 2]]


print(la.inv(U))
[[ 0.06896552 -0.17241379]
   [ 0.17241379 0.06896552]]

6

다음 부분공간 W의 직교 기저를 발견하세요?
$W=\left\{ \left[\begin{matrix}3\\-1\\2\\-1\end{matrix}\right], \quad \left[\begin{matrix}-5\\9\\-9\\3\end{matrix}\right] \right\}$

x1=np.array([[3, -1, 2, -1]]).reshape(4,1)
x2=np.array([[-5, 9, -9, 3]]).reshape(4,1)
print(x1)
[[ 3]
   [-1]
   [ 2]
   [-1]]


print(x2)
[[-5]
   [ 9]
   [-9]
   [ 3]]


v1=x1
v2=x2-np.dot(v1.T, x2)/np.dot(v1.T, v1)*v1
print(v2)
[[ 4.]
   [ 6.]
   [-3.]
   [ 0.]]

7

다음 행렬은 정규직교입니까?

A=Matrix([[Rational('4/5'), 0, Rational('3/5')], [Rational('3/5'), 0, -Rational('4/5')],[0,1,0]])
A
$ \left[\begin{matrix}\frac{4}{5} & 0 & \frac{3}{5}\\\frac{3}{5} & 0 & - \frac{4}{5}\\0 & 1 & 0\end{matrix}\right]$

A.det()
1
A.T
$\left[\begin{matrix}\frac{4}{5} & \frac{3}{5} & 0\\0 & 0 & 1\\\frac{3}{5} & - \frac{4}{5} & 0\end{matrix}\right]$

A.inv()
$ \left[\begin{matrix}\frac{4}{5} & \frac{3}{5} & 0\\0 & 0 & 1\\\frac{3}{5} & - \frac{4}{5} & 0\end{matrix}\right]$

위 결과와 같이 행렬은 선형독립이며 전치행렬과 역행렬이 같으므로 정규직교 행렬입니다.

2021년 1월 2일 토요일

미적분이야기_CH2

.  액체 1.5L를 담을 수 있는 원통형 캔을 만듭니다. 최소의 재료를 사용할 캔의 차원을 결정?

캔의 높이: h, 반지름:r
부피 v = πr2h = 1.5 → h = 1.5/(πr2)
최소의 재료를 사용하여야 한다는 조건은 캔의 표면적을 최소로 하는
것과 같습니다.
표면적 A = 2πrh + 2πr2
r=symbols('r', real=True)
h=1.5/(pi*r**2)
A=2*pi*r*h+2*pi*r**2
A
$2\pi r^2+\frac{3.0}{r}$
dAdr=A.diff(r)
dAdr
$4\pi r - \frac{3.0}{r^2}$
sol=solve(dAdr, r)
sol
[0.620350490899400]
dA2dr=A.diff(r, 2)
dA2dr
$4\pi+\frac{6.0}{r^3}$
dAForm=FunIncDecS(dAdr, r, sol, rev=0.5, n=3)[0]
dA2Form=FunIncDecS(dA2dr, r, sol, rev=0.5, n=3)[0]
re=pd.concat([dAForm, dA2Form.loc[:,1]], axis=1)
re.columns=["r", "dA부호","dA2부호"]
re
r dA부호 dA2부호
0 0.120 - +
1 0.620 0 +
2 1.12 + +

임계점에서 A'(r)의 부호는 '-'에서 '+'로 전환되며 A''(r)>0 이므로 오목형입니다. 그러므로 이 점에 대응하는 표면적이 최소면적이 됩니다. 위의 결과 r을 대입하여 높이와 면적을 결정합니다.

N(h.subs(r, sol[0]), 3)
1.24

N(A.subs(r, sol[0]), 3)
7.25


.  다음 그림과 같은 창을 만들려고 합니다. 창틀 재료의 길이는 총 12 m입니다. 가장 크게 만들기 위한 창의 차원을 결정하세요?

창의 면적: $A=\frac{1}{2}πr^2+2hr$

조건: 12 = 2h + 2r + πr
h,r=symbols('h r', real=True)
eq=Eq(12, 2*h+2*r+pi*r)
h1=solve(eq, h)
h1
[-pi*r/2 - r + 6]

A=1/2*pi*r**2+2*h1[0]*r
A
$0.5\pi r^2+r(−\pi r−2r+12)$

dAdr=A.diff(r)
dAdr
$r(−\pi−2)−2r+12$

sol=solve(dAdr, r)
sol
[12/(π + 4)]

dA2dr=A.diff(r, 2)
dA2dr
−(π+4)

dAForm=FunIncDecS(dAdr, r, sol, rev=0.5, n=3)[0]
dA2Form=FunIncDecS(dA2dr, r, sol, rev=0.5, n=3)[0]
re=pd.concat([dAForm, dA2Form.loc[:,1]], axis=1)
re.columns=["r", "dA부호","dA2부호"]
re
w dA부호 dA2부호
0 1.18 + -
1 1.68 0 -
2 2.18 - -

임계점을 기준으로 A'(r)의 부호 전환, A"(r)은 음의 부호를 가지므로 볼록형이므로 임계점에 대응하는 값은 극대값이 됩니다. 즉, 그 지점에서는 최대 높이와 최대면적이 형성됩니다.

N(h1[0].subs(r, sol[0]), 3)
1.68

N(A.subs(r, sol[0]), 3)
10.1


.  y = x2 + 1의 곡선 위의 점 중에서 (0, 2)와 가장 가까운 점을 결정하세요?

곡선 위의 점을 (x, y)라고 하면 점 (0, 2)와의 거리는 다음과 같이 계산할 수 있습니다.
$d=\sqrt{(x-0)^2+(y-2)^2}=\sqrt{x^2+(y-2)^2}$

x, y=symbols('x y', real=True)
L=sqrt(x**2+(y-2)**2)
L
$\sqrt{x^2+(y-2)^2}$

y1=x**2+1
L=L.subs(y, y1)
L
$\sqrt{x^2+(x^2-1)^2}$

dLdx=L.diff(x)
dLdx
$\frac{2x(x^2-1)+x}{\sqrt{x^2+(x^2-1)^2}}$

#x > 0
sol=solve(dLdx, x)
sol
[0, -sqrt(2)/2, sqrt(2)/2]

dL2dx=L.diff(x, 2)
simplify(dL2dx)
$\frac{2x^6-3x^4+6x^2-1}{\sqrt{(x^4-x^2+1)^3}}$

dLForm=FunIncDecS(dLdx, x, sol[2], rev=0.2, n=3)[0]
dL2Form=FunIncDecS(dL2dx, x, sol[2], rev=0.2, n=3)[0]
re=pd.concat([dLForm, dL2Form.loc[:,1]], axis=1)
re.columns=["x", "dL부호","dL2부호"]
re
x dL부호 dL2부호
0 0.507 - +
1 0.707 0 +
2 0.907 + +

N(y1.subs(x, sol[2]), 3)
1.5
N(L.subs(x, sol[2]), 3)
0.866

2020년 12월 9일 수요일

미적분 이야기_CH1

1

1.  x=1에서 f(x)=x2의 접선의 기울기를 추정하기 위해
(a)(1,1), (2,4)를 통과하는 외선의 기울기?
(b)(1,1), (3/2, 9/4)를 통과하는 외선의 기울기?
(c)x=1에서의 접선의 기울기는 어떤 외선의 기울기에 근접할까요.

x=symbols('x')
f=x**2
f
x2
#(a)
a_slope=(1-4)/(1-2)
a_slope
3.0
#(b)
b_slope=(1-9/4)/(1-3/2)
b_slope
2.5
#(c)
limit(f, x, 1)
1

2

시간(t)에 따른 거리에 대한 식 s(t)=-16t2+64로부터 t=0.5초에서의 순간속도?

x, h=symbols('x, h')
f=-16*x**2+64
f2=(f.subs(x, x+h)-f)/h
f2
$\frac{16 x^{2} - 16 \left(h + x\right)^{2}}{h}$

v_avg=limit(f2, h, 0)
v_avg
−32x

v_avg.subs(x, 0.5)
−16.0

6

다음 극한의 연속 또는 불연속을 결정하세요?
(1)$\lim\limits_{x \to 2} \large \frac{x-3}{x(x-2)}$

x, h=symbols('x, h')
f=(x-3)/(x*(x-1))
f
$\frac{x - 3}{x \left(x - 1\right)}$
limit(f, x, 1, "-")

limit(f, x, 1, "+")
−∞

(2)$\lim\limits_{x \to 1} \large \frac{x+3}{(x-1)^2}$

x, h=symbols('x, h')
f=(x+2)/(x-1)**2
f
$\frac{x + 2}{\left(x - 1\right)^{2}}$
limit(f, x, 1, "-")

limit(f, x, 1, "+")

(3)$\frac{x-1}{x^2+2x}=\frac{x-1}{x(x+2)}$
x=-2, 0에서 분모가 0이 됩니다. 그러므로 $(-\infty, -2)\; \cup \;(2, 0)\; \cup \;(0, \infty)$구간에서 연속

7

다음 극한을 계산하세요?
극한 계산에서 함수들은 가장 간단한 형태로 전환후에 실시합니다.

(1)$\lim\limits_{x \to 3} \large \frac{x^2-3x}{x^2-5x-3}$

x=symbols('x')
f=(x**2-3*x)/(2*x**2-5*x-3)
f
$\frac{x^{2} - 3 x}{2 x^{2} - 5 x - 3}$
f1=simplify(f)
f1
$\frac{x}{2 x + 1}$
limit(f1, x, 3)
$\frac{3}{7}$
limit(f, x, 3)
$\frac{3}{7}$

(2)$\lim\limits_{ x \to -1} \large \frac{\sqrt{x + 2} - 1}{x + 1}$

$ \begin{matrix}\frac{\sqrt{x + 2} - 1}{x + 1}&=\frac{(\sqrt{x + 2} -1)(\sqrt{x + 2} +1)}{(x + 1)(\sqrt{x + 2} +1)}\\&=\frac{x+2-1}{(x + 1)(\sqrt{x + 2} +1)}\\&=\frac{1}{\sqrt{x + 2} +1}\end{matrix}$

$\lim\limits_{x \to -1} \large \frac{1}{\sqrt{x + 2} +1}=\frac{1}{2}$

x=symbols('x')
f=(sqrt(x+2)-1)/(x+1)
limit(f, x, -1)
$\frac{1}{2}$

(3)$\lim\limits_{ x \to 5} \large \frac{\sqrt{x - 1} - 2}{x - 5}$

$ \begin{matrix}\frac{\sqrt{x - 1} - 2}{x - 5}& =\frac{(\sqrt{x - 1} -2)(\sqrt{x - 1} +2)}{(x-5)(\sqrt{x- 1}+2)}\\ &=\frac{x-1-4}{(x -5)(\sqrt{x-1} +2)}\\ &=\frac{1}{\sqrt{x-1} +2}\end{matrix}$

$\lim\limits_{x \to 5} \large \frac{1}{\sqrt{x-1} +2}=\frac{1}{4}$

x=symbols('x')
f=(sqrt(x-1)-2)/(x-5) limit(f, x, 5)
$\frac{1}{4}$

(4)$\lim\limits_{ x \to 1} \large \frac{\frac{1}{x + 1}-\frac{1}{2}}{x - 1}$

$ \begin{matrix}\frac{\frac{1}{x + 1}-\frac{1}{2}}{x - 1}&=\frac{\frac{-(x-1)}{2(x+1)}}{x-1}\\&=\frac{-(x-1)}{2(x+1)(x-1)}\\&=\frac{-1}{2(x+1)}\end{matrix}$

$\lim\limits_{x \to 1} \large -\frac{1}{2(x+1)}=-\frac{1}{4}$

x=symbols('x')
f=(1/(x+1) -1/2)/(x-1)
limit(f, x, 1)
$−\frac{1}{4}$

8

다음 함수의 점근선을 결정하세요?

(1)$\frac{3x-1}{2x+5}$

$\begin{align}\lim\limits_{x \to \infty}\frac{3-\frac{1}{x}}{2+\frac{5}{x}}=\frac{3}{2}\\ \lim\limits_{x \to -\infty}\frac{3-\frac{1}{x}}{2+\frac{5}{x}}=\frac{3}{2}\end{align}$

(2)$\frac{3x^2-2x}{4x^3-5x+7}$

$\begin{align}\lim\limits_{x \to \infty}\frac{\frac{3}{x}-\frac{2}{x^2}}{4-\frac{5}{x^2}+\frac{7}{x^3}}=0\\\lim\limits_{x \to -\infty}\frac{\frac{3}{x}-\frac{2}{x^2}}{4-\frac{5}{x^2}+\frac{7}{x^3}}=0\end{align}$

9

다음을 계산하세요?

$(1) \lim\limits_{x \to 0} x^2 cos\left(\frac{1}{x}\right)$
$t=\frac{1}{x},\quad x \rightarrow 0 \Rightarrow t \rightarrow \infty$
$\begin{align}\lim\limits_{x \to 0} x^2 cos\left(\frac{1}{x}\right)&=\lim\limits_{t \to \infty}\frac{1}{t^2}cos(t)\\&=0\end{align}$

x=symbols('x')
f=x**2*cos(1/x)
limit(f, x, 0)
0

$(2)\quad\begin{align}\lim\limits_{x \to 0} \frac{tan(x)}{x}&=\lim\limits_{x \to 0}\frac{1}{x} \frac{sin(x)}{cos(x)}\\&=\lim\limits_{x \to 0}\frac{sin(x)}{x}\lim\limits_{x \to 0}\frac{1}{cos(x)}\\&=1\end{align}$

x=symbols('x')
f=tan(x)/x
limit(f, x, 0)
1

$(3)\quad\lim\limits_{x \to \infty} xsin\left(\frac{1}{x}\right)$
$\frac{1}{x}=t, \quad x \to \infty \Rightarrow t \to 0$
$\begin{align}\lim\limits_{x \to \infty} xsin\left(\frac{1}{x}\right)&=\lim\limits_{t \to 0}\frac{sin(t)}{t}\\&=1\end{align}$

x=symbols('x')
f=x*sin(1/x)
limit(f, x, oo)
1

$(4)\quad\begin{align}\lim\limits_{x \to 0}\frac{x+2sin(x)}{sin(x)cos(x)}&=\lim\limits_{x \to 0}\frac{x}{sin(x)cos(x)}+\lim\limits_{x \to 0}\frac{2sin(x)}{sin(x)cos(x)}\\ &=\lim\limits_{x \to 0}\frac{1}{\frac{sin(x)}{x}}\lim\limits_{x \to 0}\frac{1}{cos(x)}+\lim\limits_{x \to 0}\frac{2}{cos(x)}\\ &=1 \cdot 1+ 2\\&=3\end{align}$

x=symbols('x')
f=(x+2*sin(x))/(sin(x)*cos(x))
limit(f, x, 0)
3

9

다음을 계산하세요?

$(1) \lim\limits_{x \to 0} x^2 cos\left(\frac{1}{x}\right)$
$t=\frac{1}{x},\quad x \rightarrow 0 \Rightarrow t \rightarrow \infty$
$\begin{align}\lim\limits_{x \to 0} x^2 cos\left(\frac{1}{x}\right)&=\lim\limits_{t \to \infty}\frac{1}{t^2}cos(t)\\&=0\end{align}$

x=symbols('x')
f=x**2*cos(1/x)
limit(f, x, 0)
0

$(2)\quad\begin{align}\lim\limits_{x \to 0} \frac{tan(x)}{x}&=\lim\limits_{x \to 0}\frac{1}{x} \frac{sin(x)}{cos(x)}\\&=\lim\limits_{x \to 0}\frac{sin(x)}{x}\lim\limits_{x \to 0}\frac{1}{cos(x)}\\&=1\end{align}$

x=symbols('x')
f=tan(x)/x
limit(f, x, 0)
1

$(3)\quad\lim\limits_{x \to \infty} xsin\left(\frac{1}{x}\right)$
$\frac{1}{x}=t, \quad x \to \infty \Rightarrow t \to 0$
$\begin{align}\lim\limits_{x \to \infty} xsin\left(\frac{1}{x}\right)&=\lim\limits_{t \to 0}\frac{sin(t)}{t}\\&=1\end{align}$

x=symbols('x')
f=x*sin(1/x)
limit(f, x, oo)
1

$(4)\quad\begin{align}\lim\limits_{x \to 0}\frac{x+2sin(x)}{sin(x)cos(x)}&=\lim\limits_{x \to 0}\frac{x}{sin(x)cos(x)}+\lim\limits_{x \to 0}\frac{2sin(x)}{sin(x)cos(x)}\\ &=\lim\limits_{x \to 0}\frac{1}{\frac{sin(x)}{x}}\lim\limits_{x \to 0}\frac{1}{cos(x)}+\lim\limits_{x \to 0}\frac{2}{cos(x)}\\ &=1 \cdot 1+ 2\\&=3\end{align}$

x=symbols('x')
f=(x+2*sin(x))/(sin(x)*cos(x))
limit(f, x, 0)
3

10

다음을 계산하세요?

$(1)\lim\limits_{n \to \infty} \large n \cdot ln\left(1+\frac{x}{n}\right)$

$u=\frac{x}{n}, \quad n \to \infty \Rightarrow u \to 0$
$\begin{align}\lim\limits_{n \to \infty} \large n \cdot \left(1+\frac{x}{n}\right)&=\lim\limits_{u \to 0} \frac{x}{u} \cdot ln(1+u)\\&=x\lim\limits_{u \to 0} ln(1+u)^{\frac{1}{u}}\\&=x \end{align}$

n, x=symbols('n x')
f=n*ln(1+x/n)
limit(f, n, oo)
x

$(2)\begin{align}\lim\limits_{n \to \infty} \frac{2x-7}{\sqrt{9x^2+7x+2}}&=\lim\limits_{n \to \infty}\frac{2-\frac{7}{x}}{\sqrt{9+\frac{7}{x}+\frac{2}{x^2}}}\\&=\frac{2}{3} \end{align}$

x=symbols('x')
f=(2*x-7)/(sqrt(9*x**2+7*x+2))
limit(f, x, oo)
$\frac{2}{3}$

(3) $\begin{align}\lim\limits_{x \to \infty} \sqrt{x^{2} + 5 x + 9}- \sqrt{x^{2} - 4 x + 2} &=\lim\limits_{x \to \infty}\frac{(x^{2} + 5 x + 9)-(x^{2} - 4 x + 2)}{\sqrt{x^{2} + 5 x + 9}+\sqrt{x^{2} - 4 x + 2}}\\&=\lim\limits_{x \to \infty}\frac{ 9x + 7}{\sqrt{x^{2} + 5 x + 9}+\sqrt{x^{2} - 4 x + 2}}\\&=\lim\limits_{x \to \infty}\frac{ 9 + \frac{7}{x}}{\sqrt{1 + \frac{5}{x} +\frac{9}{x^2}}+\sqrt{1 - \frac{4}{x}+\frac{2}{x^2}}}\\&=\frac{9}{2}\end{align}$

x=symbols('x')
f=sqrt(x**2+5*x+9)-sqrt(x**2-4*x+2)
limit(f, x, oo)
$\frac{9}{2}$

2020년 12월 7일 월요일

선형대수코드_CH 7

1

행렬 A의 QR분해를 실행하세요?

A=np.array([[-3, -19, 3],[9, -7, 11], [-3, 9, 11], [0, 4, 32]])
print(A)
[[ -3 -19 3]
   [ 9 -7 11]
   [ -3 9 11]
   [ 0 4 32]]


Q,R=la.qr(A)
print(np.around(Q, 3))
[[-0.302 -0.898 -0.265]
   [ 0.905 -0.18 -0.193]
   [-0.302 0.359 -0.313]
   [ 0. 0.18 -0.891]]


print(np.around(R, 3))
[[ 9.95 -3.317 5.729]
   [ 0. 22.271 5.029]
   [ 0. 0. -34.884]]


print(np.dot(Q, R))
[[ -3. -19. 3.]
   [ 9. -7. 11.]
    [ -3. 9. 11.]
   [ 0. 4. 32.]]

2

다음 각 이차식에 대응하는 대칭행렬을 나타내십시요.

a. ℝ2의 이차식, $3x_1^2-4x_1x_2+5x_2^2$
$\left[\begin{matrix}3&-2\\-2&5\end{matrix}\right]$

b.ℝ2의 이차식, $3x_1^2+2x_1x_2$
$\left[\begin{matrix}3&1\\1&0\end{matrix}\right]$

c.ℝ3의 이차식, $3x_1^2+2x_1x_2$
$\left[\begin{matrix}3&1&0\\1&0&0\\0&0&0\end{matrix}\right]$

3

$7x_1^2+8x_1x_2+4x_1x_3+10x_2^2+5x_2x_3+8x_3^2$에 대응하는 대칭행렬을 나타내십시요. 이 식은 양의 정부호 입니까?

A=np.array([[7, 4, 2],[4, 10, 2.5],[2, 2.5, 8]])
print(A)
[[ 7. 4. 2. ]
   [ 4. 10. 2.5]
   [ 2. 2.5 8. ]]


la.eig(A)
(array([14.37376873, 4.20939986, 6.41683141]),
   array([[ 0.51711708, 0.83649609, -0.1812849 ],
      [ 0.72914595, -0.54145794, -0.41853253],
      [ 0.44825897, -0.08424717, 0.88992489]]))

위 식에 대응하는 대칭행렬의 고유값은 모두 양수입니다. 그러므로 위 식은 양의 정부호 입니다.

4

다음은 xTAx가 yTDy의 변환입니다. D를 발견하고 성립함을 보이세요.
$5 x_{1}^{2} + 4 x_{1} x_{2} + 6 x_{2}^{2} - 4 x_{2} x_{3} + 7 x_{3}^{2}=3.0 y_{1}^{2} + 6.0 y_{2}^{2} + 9.0 y_{3}^{2}$

x1,x2,x3=symbols("x1 x2 x3")
x=Matrix(3,1, [x1,x2,x3])
x
$\left[\begin{matrix}x_{1}\\x_{2}\\x_{3}\end{matrix}\right]$
A=Matrix([[5,2,0],[2,6,-2],[0,-2,7]])

A
$\left[\begin{matrix}5 & 2 & 0\\2 & 6 & -2\\0 & -2 & 7\end{matrix}\right]$
simplify(x.T*A*x)

$\left[\begin{matrix}5 x_{1}^{2} + 4 x_{1} x_{2} + 6 x_{2}^{2} - 4 x_{2} x_{3} + 7 x_{3}^{2}\end{matrix}\right]$

d, p=la.eig(np.array(A, dtype=np.float64))
d
array([3., 6., 9.])
D=np.diag(d)
print(D)
[[3. 0. 0.]
   [0. 6. 0.]
   [0. 0. 9.]]


y1,y2,y3=symbols("y1 y2 y3")
y=Matrix(3,1,[y1,y2,y3])
y.T*Matrix(D)*y
$\left[\begin{matrix}3.0 y_{1}^{2} + 6.0 y_{2}^{2} + 9.0 y_{3}^{2}\end{matrix}\right]$

5

Q(x)에 대해 (a)xTx=1의 제한 조건에서 최대값을 계산하고 이 값들이 각각 Q(x)의 이차식 행렬 A의 최대 고유값과 같다는 것을 보이세요. (b)xTx=1와 xTu1=0의 제한 조건에서 최대값(u1은 최대 고유값에 대응하는 고유단위벡터)을 계산하세요.
$Q(x)=3 x_{1}^{2} + 2 x_{1} x_{2} + 2 x_{2}^{2} + 4 x_{2} x_{3} + 2 x_{3}^{2}$

x1,x2,x3=symbols("x1 x2 x3")
x=Matrix(3,1, [x1,x2,x3])
x
$\left[\begin{matrix}x_{1}\\x_{2}\\x_{3}\end{matrix}\right]$

A=Matrix([[3,1,0],[1,2,2],[0,2,2]])
A
$\left[\begin{matrix}3 & 1 & 0\\1 & 2 & 2\\0 & 2 & 2\end{matrix}\right]$

xTAx=x.T*A*x
expand(xTAx[0])
$3 x_{1}^{2} + 2 x_{1} x_{2} + 2 x_{2}^{2} + 4 x_{2} x_{3} + 2 x_{3}^{2}$

d, p=la.eig(np.array(A, dtype=np.float64))
d
array([-0.16424794, 2.77286556, 4.39138238])

print(p)
[[ 0.2260912 0.84604119 0.48280128]
   [-0.7154086 -0.19216509 0.6717612 ]
   [ 0.6611152 -0.49727948 0.56181831]]


u1=p[:,2]
la.norm(u1)
1.0

np.dot(np.dot(u1.T,np.array(A, dtype=np.float64)),u1)
4.391382380630901

u2=p[:,1]
u2
array([ 0.84604119, -0.19216509, -0.49727948])

np.dot(np.dot(u2.T,np.array(A, dtype=np.float64)),u2)
2.77286555782931

6

다음 행렬 A 의 특이값?
$A=\left[\begin{matrix}2&3\\0&2\end{matrix}\right]$

A=np.array([[2,3],[0,2]])
print(A)
[[2 3]    [0 2]]

AA=np.dot(A.T, A)
print(AA)
[[ 4 6]    [ 6 13]]

d, P=la.eig(AA)
print(d)
[ 1. 16.]

a=np.sqrt(d)
a
array([1., 4.])

특이값은 1, 4

7

다음 행렬 A의 특이값 분해에서 좌특이행렬, 특이값, Σ행렬, 우특이행렬을 계산하세요?
$A=\left[\begin{matrix}3 & -3\\0 & 0\\1 & 1\end{matrix}\right]$

A=np.array([[3, -3],[0,0],[1,1]])
AA=np.dot(A.T,A)
print(AA)
[[10, -8],
[-8, 10]]


d, V=la.eig(AA)
d
array([18., 2.])

특이값, 특이행렬, 우려특이행렬

s=np.sqrt(d)
s
array([4.24264069, 1.41421356])
print(V)
[[ 0.70710678 0.70710678]
   [-0.70710678 0.70710678]]


sMat=np.vstack([np.diag(s), np.array([0,0])])
print(sMat)
[[4.24264069 0. ]
   [0. 1.41421356]
   [0. 0. ]]

좌특이행렬

AV=np.dot(A, V)
print(AV)
[[ 4.24264069e+00 -8.88178420e-16]
   [ 0.00000000e+00 0.00000000e+00]
   [ 2.22044605e-16 1.41421356e+00]]


U1=AV/s
print(np.around(U1, 3))
[[ 1. -0.]    [ 0. 0.]    [ 0. 1.]]

U1은 3×3차원이 되어야 합니다. 그러나 위의 V로 부터의 산출한 U는 3×2차원이므로 각 열벡터와의 선형 결합 즉, U1Tx=0가 선형독립이고 정규직교기저인 벡터를 첨가하여 U를 생성합니다.
$\left[\begin{matrix}1&0&0\\0&0&0\\0&0&1\end{matrix}\right]\left[\begin{matrix}x_{1}\\x_{2}\\x_{3}\end{matrix}\right]=\left[\begin{matrix}0\\0\\0\end{matrix}\right]$
위 식의 해집합의 기저를 발견하기 위해 sympy의 nullspace()함수를 사용합니다.

u2=Matrix(U1.T).nullspace()
u2=np.array(u2, dtype=np.float64).reshape(3,1)
print(u2)
[[0.]
   [1.]
   [0.]]

결과를 U1에 첨가하면 3×3차원의 U가 생성됩니다.

U=np.hstack([U1, u2])
print(np.around(U, 3))
[[ 1. -0. 0.] [ 0. 0. 1.] [ 0. 1. 0.]]

확인

print(np.around(np.dot(np.dot(U,sMat),V.T),3))
[[ 3. 3.]
   [ 0. 0.]
   [-1. 1.]]

numpy.linalg.syd()함수를 사용하여 각 요소들을 계산

U,s,V=la.svd(A)
print(np.around(U, 3))
[[-1. -0. 0.]
   [ 0. 0. -1.]
   [ 0. -1. 0.]]

print(s)
[4.24264069 1.41421356]
print(V)
[[-0.70710678 0.70710678]
   [-0.70710678 -0.70710678]]

CalculusMadeEasy_CH1

1-6 a=np.array([[1],[2],[3]]) la.norm(a) 3.7416573867739413 b=np.array([[-1],[-3],[-4]]) la.norm(b) 5.0990195135927845 c=np.array([[...