C - Fakepathc - fakepathALL LECTURES LA and AG
C - Fakepathc - fakepathALL LECTURES LA and AG
C - Fakepathc - fakepathALL LECTURES LA and AG
Lecture 1
Complex Numbers and operations on complex numbers
Imaginary Unit
All of the normal operations of mathematics – addition, subtraction, multiplication and division –
can be applied to the number i .
Addition:
i +i=2 i , 2 i+i =3i , 5 i+ 3i=8i
Subtraction:
i−i=0 ,10 i−i=9 i , 3 i−5 i=−2
2 2
Multiplication:
i⋅i=i 2 =−1 , (3 i)⋅(5 i)=15i =−15 ,(4 i)⋅(6 i)=24 i =−24
i 5i 2i 2
=1 =5 =
Division:
i , i ,
7 i 7
Powers of the Imaginary unit
i 4 n =1 , i 4 n±2=−1
Complex Numbers
Definition. A complex number is an expression of the form
For the real and imaginary parts we use the notation x=Re z and y=Im z .
Definition. The complex conjugate of the complex number is the number with an equal real part and
an imaginary part equal in magnitude but opposite in sign.The complex conjugate of the complex
number z=x +iy is given by z̄=x −iy
For example:
Definition. Two complex numbers are equal if and only if their real parts are equal and their
imaginary parts are equal: if
z 1=x 1 + iy 1 and z 2 =x2 +iy 2 then the equalities
{
x 1 =x 2
z 1=z 2 ⇔ y 1= y 2 are equivalent.
⇒ x=3
Solution. 2x=6
( 4 +2 i ) x + ( 5−3 i ) y=13+i
Solution. Select the real and imaginary parts
Hence
{4 x +5 y=13
2 x−3 y=1 Solving this system, we find x=2 , y=1
Complex numbers may be represented in a plane. To do this we take a Cartesian coordinate system
x , y which enables us to represent any number of the form z=x +iy as a point.
Real numbers are a special case of complex numbers: if we put y=0 in formula (1) we may regard
a complex number x +i⋅0 as representing real number x ; real numbers are represented as points
lying on the real axis i.e. on the x -axis. Complex numbers that not real are called imaginary. Thus,
every complex number is either real or imaginary. A complex number without a real part ( i.e. with
the real part equal to zero) is called pure imaginary; such numbers are represented by points of the
y -axis which is called imaginary axis .
It is often convenient to introduce polar coordinates in a complex plane. The polar coordinates ρ
and θ of the point, which represents the complex number z=x +iy are denoted as ρ=|z| and
θ=arg z . ρ is called the modulus or the absolute value of z and θ is called the argument, or
amplitude ,or phase, of z.
As is known, ρ=√ x 2 + y 2 , x=ρ cos θ , and y= ρsin θ .This implies, by (1) , that
z=ρ( cos θ+i sinθ ) (2)
Hence, each complex number can be written in the so-called trigonometric form (2).
The modulus of a complex number is a certain uniquely defined non-negative real number whereas
the argument is defined within an integral multiple of 2 π . For instance, |i|=1 and
π
Argi= +2 kπ
2 ( k=0 ,±1 ,±2 ,. .. ) . The sign Argz denotes the totality of all the possible values
of the argument of a complex number z .Thus, Argz has infinitely many different values. There is,
however, one and only one value of Argz , denoted as arg z ,which satisfies the inequality
−π <arg z ≤π ; arg z is called the principal value of the argument.
{
y
arctan ,if x >0
x
y
π + arctan,if x <0 , y≥0
x
y
arg z = −π + arctan , if x <0 , y <0
x
π
,if x=0 , y >0
2
π
− ,if x=0 , y <0
2
ρ=|z|=√ 22 +2 2=√ 8
y 2 π
tanθ= = =1 θ=arctan 1=
x 2 ⇒ 4
Example 4
y −2
tanθ= = =−1 θ=arctan (−1 ) =− π
ρ=|z|=√ 2 + (−2 ) = √8=2 √2 ,
2 2
x 2 , 4
z=2 √ 2[cos − ( ) π
4 ( )
π
+i sin − ]
4
Example 5
tanθ= = √ = √3
y − 3
√
ρ=|z|= (−1 ) + (−√3 ) =2 ,
2 2
x −1 ,
π −2 π
θ=−π + arctan √ 3=−π + =
3 3
z z
1.Addition .Two complex numbers 1 and 2 are most easily added by separately adding their real
and imaginary parts of the summands. That is to say:
z 1 + z 2 =( x 1 +iy 1 ) + ( x2 +iy 2 ) =( x 1 + x 2 ) +i ( y 1 + y 2 )
z z̄=x 2 + y 2
If we use the trigonometric form z 1= ρ1 ( cos θ 1 +i sinθ 1 ) z 2 =ρ2 ( cos θ2 +i sin θ2 ) then the
,
product
z 1 z 2 can be written as
Therefore,
ρ=ρ1 ρ2 , θ=θ1 θ2 i.e. |z 1 z 2|=|z 1|⋅|z 2|, Arg ( z 1 z 2 ) = Argz1 + Argz 2
Hence, when complex numbers are multiplied their moduli are multiplied and their arguments are
added.
4.Reciprocal and division. Using the conjugation, the reciprocal of a nonzero complex number
z=x +iy can always be broken down to
1 z̄ z̄ z̄ x y
= = 2 = 2 2 = 2 2 −i 2 2
z z z̄ |z| x + y x + y x +y
z1 z 1 z̄ 2 z 1 z̄ 2 z1 x 1 x 2 + y 1 y 2 x 2 y 1 −x 1 y 2
= = = +i
z2 z 2 z̄ 2 |z 2|2 or
z2 x 22 + y 22 x 22 + y 22
z1 ρ1
z2 ρ2 [
= cos ( θ1 −θ2 ) +i sin ( θ1 −θ2 ) ]
In trigonometric form
( )
z |z | z
| 1|= 1 Arg 1 = Argz1 −Argz 2
z |z 2| ,
i.e. 2
z2
5. The real part Re z and the imaginary part Im z of the complex number z are expressed through
conjugate complex numbers as follows:
z̄ + z z̄ −z z− z̄
Re z= Im z=i =
2 , 2 2i
i 4 =1 , i 5 =i , i 6 =−1 , i 7 =−i ,
ix
and so on, and by considering the Taylor series expansions of e ,cos x andsin x :
( x2 x4 x6 x8
2 ! 4 ! 6 ! 8!
x 3 x 5 x7
= 1− + − + −.. . . +i x− + − +. ..
3! 5! 7! ) ( )
=cos x+ i sin x
The rearrangement of terms is justified because each series is absolutely convergent.
The rule of multiplication of complex numbers extends automatically to an arbitrary number of
factors. In particular, if we take equal factors then we have
[ ρ ( cos θ+i sin θ ) ]n =ρn ( cos nθ+i sin nθ ) ( n=2 ,3 , . .. )
,
In case ρ=1 , we obtain the formula ( cos θ+ i sin θ )n =( cos nθ+i sin nθ )
Which was named De Moivre`s formula.
60
Example 6. Calculate ( −1+i √ 3 ) .
θ+2 kπ
n ψ=
Since r and ρ are non-negative numbers, we have r =√ ρ and n .
Thus,
n n
√ z=√ ρ cos(θ+2 kπ
+i sin
n
θ+ 2 kπ
n ) ,where k=0,1,2,3,…,n-1
4
Example 7. Find all values of √ 1−i .
Solution. We represent complex number 1−i in a trigonometric form
[ ( ) ( )]
1−i= √ 2 cos −
π
4
+i sin −
π
4
[ ( ) ( )]
π π
− +2 kπ − +2 kπ
4 8 4 4
√ 1−i= √2 cos + i sin
4 4
Hence,
( k =0 )
4 8
(π
√ 1−i= √2 cos −isin
π
16 16 )
( k =1 )
4 8
(
7π
√ 1−i= √2 cos +isin
7π
16 16 )
( k =2 )
4 8
(
√ 1−i= √2 cos
15 π
+isin
16
15 π
16 )
( k =3 )
4 8
√ 1−i= √2 cos(23 π
16
+ isin
23 π
16 )
Examples:
z2
z =3+4 i , z 2 =1−2 i ,calculate: z1
Example 1. Given that 1
Solution.
z2 z 2 z̄ 1 z 2 z̄ 1
= =
z1 z 1 z̄ 1 |z 1|2
Example 2.
Find modulus and the argument of z :
π π
z=−sin −icos
8 8
Solution. We have
π π
x=−sin <0 y=−cos <0
8 8
,
In this case principal value of the argument is
y
arg z =−π +arctan
x
π
−cos
arg z =−π + arctan
−sin
8
π
π
=−π +arctan cot =
8 ( )
8
[ ( )]
π π
=−π + arctan tan − =−π +arctan tan
2 8
3π
8 (
=−π + =
8 8 )
3 π −5 π
5π
Argz=− +2 kπ
Hence, 8 ( k=0,±1,±2,... )
Absolute value of z is
π
√ π
ρ=|z|= sin2 + cos2 =1
8 8
Example 3.
a)-2 b) 2 i
Solution.
ρ=|z|=2
In this case principal value of the argument is
y
arg z =π + arctan =π
x
Hence, trigonometric form is z=ρ( cos θ+i sinθ )=2 ( cos π +i sin π )
ρ=|z|=2
In this case principal value of the argument is
π
arg z =
2
( )
40
1+i √ 3
Calculate: 1−i
Solution.
1+i √3
Let`s consider
¿ 240
[ 3 3 ]
cos (− )+ isin (− ) =2 (− −i √ )
2π 2π 1
2 2
3
2 20
49
√
ρ=|z|= (−1 )2 + (−√3 ) =2
2
,
In this case principal value of the argument is
Test
1.Given that
z 1=3+4 i , z 2 =1−2 i ,calculate:
z1 z2
z 1−z 2 b) 2 c)
|z 1| d)
z1
a)
form: a)
2.Perform the indicated operation and write the answers in standard
(−4+7 i ) + ( 5−10 i ) b) ( 4 +12 i )−( 3−15 i ) c) 5 i−(−9+i )
3.Multiplay each of the folowing and write the answers in standard form:
π π
z=−sin −icos
a) 8 8 b) z=4+3 i c) z=−2+2 √ 3 i
5.Write complex numbers in trigonometric form:
a)-2 b) 2 i c) − √ 2+i √ 2
6. Write complex numbers in exponential form:
a) i b) −1−i √3 c) -2
7. Calculate:
( )
40
1+i √ 3
6
a) ( 2−2i )7
b)
( √ 3−3 i ) c)
1−i
8. Find all values of root
a)
4
√−1 b) √i 3
c) √ −1+i
( )
z1 z
= 1
z 1 + z 2 =z 1 + z 2 z z =z 1 z 2 z2 z2
a) b) 1 2 c)
Lecture 2
Matrices
Definitions. A matrix is a rectangular array of real (or complex) numbers. We’ll only deal with real
number matrices. The individual values in the matrix are called entries. The size of the array is
written as m× n , where m− is number of rows, n is number of columns.
A matrix can be put down in the general form as
( )
a11 a12 . .. a 1n
A= a21 a 22 . .. a 2n
. . . . .. . .. ...
am 1 am 2 . .. amn
The numbers
a ij are called entries the i indexes the rows of matrix and the j indexes the
v=( v n1 )=( v n )
column of matrix. An n×1 matrix is called a column vector, written
()
v1
v= v 2
vn
( )
a11 a12 .. . a1 n
A= a 21 a22 .. . a2 n
... .. . .. . .. .
an 1 a n2 .. . ann
If m= n the matrix is called square
In this case we have:
a =0
1) A matrix is said to be diagonal if ij
i≠ j
The diagonal matrix diag ( 1,1,...,1 ) is called identity matrix and is usually denoted by
( )
1 0 0 .. 0
0 1 0 .. 0
I n= 0 0 1 .. 0
.. .. .. .. ..
0 0 0 .. 1
A matrix whose all elements are equal to zero is called a zero matrix. In case of square matrix we
introduce the notions of principal and secondary diagonals. Principal diagonal is diagonal
a 11 , a22 , . .. a nn connecting the left uppermost element with the right lowermost element.
I.
a =b ij for all i , j
Two matrices A and B are equal if they have the same size and ij
II. If A and B are matrices of the same size then the sum A and B is defined by
C= A+ B , where c ij =aij + bij for all i , j
( a11
a21
a12
a22
a13
a23
b
)(
+ 11
b 21
b12
b22
b13
b23 )(
a +b
= 11 11
a21+ b21
a 12+b 12
a 22+b 22
a13 +b13
a23 +b23 )
III. If A is any matrix and α - is a scalar (real number), then the scalar multiplication
B=αA is defined by b ij=αa ij for all i , j
α
( a11 a12 a13
a 21 a22 a23
αa
)(
= 11
αa 12 αa 13
αa21 αa 22 αa 23 )
IV. We can also compute the difference D= A−B by summing A and (−1 ) B ;
D= A−B= A+ (−1 ) B matrix subtraction
M ( F ) , where F is the underlying field (usually
The set of all m× n matrices is demoted by m , n
R or C ). In the case where m= n we write M n ( F ) to demote the matrices of size n×n .
Properties. If
A , B ,C ∈ M m , n ( F ) and α , β ∈ F then
n 1≤i≤m
c ij =∑ aik bkj
defined by k =1 1≤ j≤ p
so
( a11 a12
a21 a22 )( b11 b12 b 13
b21 b22 b 23 )(
a b +a b a b +a b a b +a b
= 11 11 12 21 11 12 12 22 11 13 12 23
a21 b11+a22 b 21 a21 b12+a22 b22 a21 b13 +a 22 b23 )
Properties of matrix multiplication
Assume A , B ,C are matrices for which all products below make sense. Then
(1) AB≠BA
(2) ( AB ) C= A ( BC )
(3) A ( B+C )=AB+ AC
A2 = A×A
(4) if A is square matrix then A3 = A 2× A
AI n = A
(5)
I n A= A where
I n − is identity matrix
Examples
(13 24 ) B=(−10 11 )
A=
,
AB=( ) BA=(
3 4)
−1 3 2 2
−3 7 , so AB≠BA .
The Transpose of a matrix
T
Given an m× n matrix A , the transpose of A is the n×m matrix, denoted by A whose
columns are formed from the corresponding rows of A .
Determinants
Every square matrix has associated with it scalar called its determinant. Given a matrix A , we use
|A| or det ( A ) to designate its determinant. The determinant of the first order of square matrix A
The determinant of the second order of square matrix A is the scalar defined by
a a
det ( A )=| 11 12|=a11 a22−a 21 a12
a 21 a22
Properties of determinants
I. Interchanging two rows or two columns of a determinant is equivalent to multiplying the
determinant by −1
For example,
III. If the elements in any row (or column) have common factor λ , then the determinant equals the
determinant of the corresponding matrix in which λ=1 , multiplied by λ .
i th row removed
M ij=det A| th
j column removed
i+ j
( i , j ) cofactor of A is the number
C ij given by
C ij = ( −1 ) M ij
Property IX. A determinant is equal to the sum of the products of the elements of any row
or column of the determinant by their cofactors
Property X . The sum of products of the elements of any row or column of the determinant
by cofactors of corresponding elements of other row or column equals zero.
For example
a 11C 21 +a12 C22 +a 13 C 23=0 .
TEST
Matrices and Determinants
A=(
3 2)
1 4
2
2. Find the matrix A if
A)
(19 164 ) B)
(155 2010 ) (1312 169 ) (139 1216 ) (26 84 )
C) D) E)
4 5
| |
3. Calculate the determinant −1 2
A) -3 B) 11 C) -13 D) 3 E)13
4. Find A−B if
A= ( 1 0
−2 3 ) and
B= (−1 2
3 −1 )
A)
( 0 −2
1 2 ) B)
( 2 −2
1 2 ) (
C)
0 2
−5 4 ) (
D)
2 −2
−5 4 ) E)
( 2 2
5 −4 )
5. The matrices A , B and C are square ones of the same order. Each of the following must be
true EXCEPT
A) A+ B=B+ A
B) C ( A+B )=CA+CB
C) A ( BC )=( AB ) C
D) AB=BA
E) ( A+C ) + B= A+ ( B+C )
5 7
| |
6. Calculate the determinant 1 3
A) 1 B)-8 C) 22 D) -1 E) 8
A)
(5317 2637 ) (
B) 18 24
27 10 ) C)
(2510 2637 ) D)
(5332 2614 )
E)
(1726 3753 )
8. Find the matrix 3 A−4 B if
A= (−52 34 ) B=(−15 72 )
,
A)
(25 3720 ) B) (10−35 −19
4 ) C)
(1035 2720 ) D)
(60 16 )
(
2 19
E) 5 4
)
9. Find the matrix A−B if
A= (−45 −36 ) B=(−44 −35 )
,
(
9 −6
A) −8 11
) B) ( )
1 6
8 1 ( )
1 0
C) 0 1
( )
9 0
D) 8 1
(
1 −6
E) 0 1 )
()
4
B= 0
10. Find the matrix AB if A=( 1, 2, 3 ) ,
7
()
4
0
A) 25 B)
21 C) ( 5 2 10 ) D) 17 E) ( 4 21 )
A)
(23 −12 ) B) ( 34 03 ) C)
(35 −1
−1 ) D)
(12 21 )
E)
( 34 13 )
12. Find the product MN if
M= (36 −15 ) N=(−32 12 )
,
A)
(63 52 ) B) (53 61 ) C)
(−9
12 −2 )
13
D)
(12 03 )
E)
(−915 134 )
5 3 2
|−1 2 4 |
13. Calculate the determinant 7 3 6
A) 50 B) 68 C) 0 D) 72 E)1
2 5 1
|0 3 4 |
14. Calculate the determinant 3 4 1
A) 77 B) 24 C) 35 D) 2 E) 25
1 1 1
|5 3 6 |
15. Calculate the determinant 7 4 2
A) 13 B) 14 C) 49 D) 50 E) 15
0 1 1
|1 0 1 |
16. Calculate the determinant 1 1 0
A) 0 B) 1 C) 2 D) 3 E) 4
( ) ( )
1 3 1 2 1 0
A= 2 0 4 B= 1 −1 2
17. Find the product AB if
1 2 3 , 3 2 1
( ) ( )( ) ( )
8 0 7 2 3 0 3 4 1 1 0 7
16 10 4 2 0 8 3 −1 6 10 1 4
A)
13 5 7 B)
1 4 3 C) 4 4 4 D)
2 3 7
( )
8 0 7
1 6 1
E)
1 3 5
18. Find A if
A= 3 2
3 1 4 ( )
A)
(19 648 ) B) ( 4739 7886 ) C)
(39 136 ) D)
( 4130 1881 )
(
20 13
E) 17 16
)
2
19. Find f ( A ) if f ( x )=x −5 x+ 3 ,
A= ( 2 −1
−3 3 )
(
−3 9
A) 3 3
) B) ( ) 1 1
2 2 ( )
5 8
C) 1 0
( )
3 0
D) 0 3
( )
0 0
E) 0 0
n
20. Find A if
(10 11 )
A=
A)
(10 11 ) ( )
B) n 1
0 n C)
(10 n1 ) D)
(21 10 )
E)
(11 0n )
Lecture 3
Inverse Matrix
Let
A ∈ M n ( F ) . The matrix A is said to be invertible if there is a matrix B∈ M n ( F )
such that AB=BA=I .
−1
In this case B is called the inverse of A and the notation for the inverse is A .
A matrix that is not invertible is sometimes called a singular matrix, and invertible matrix is
called a nonsingular matrix.
A matrix A is invertible if and only if det A≠0 . The inverse of 2×2 matrix
A= ( ) a b
c d is
A−1 = (
1 d −b
|A| −c a ) |A|=ad −bc≠0 .
Let A be an invertible 3×3 matrix. Then
( )
C 11 C 21 C 31
1
−1
A = C C 22 C 32
det A 12
C 13 C 23 C 33 det A≠0
1
A−1 = adjA
or det .
The matrix of cofactors on the right side is called the adjugate.
For example:
( )
2 5 7
A= 6 3 4
5 −2 −3 −1
Find A .
2 5 7
Δ=|6 3 4 |=−18+100−84−105+16+90=−1
5 −2 −3
The nine cofactors are
3 4
C 11=(−1 )1+1| |=−9+8=−1
−2 −3
( ) (
−1 1 −1
)
1
−1
1 −1 1
A = 38 −41 34 −1
A = −38 41 −34
−1
−27 29 −24 or
27 −29 24
3.Continue until you form the identity matrix. Keep repeating linear row reduction operations
until the left side of your augmented matrix displays the identity matrix (diagonal of 1s, with other
terms 0). When you have reached this point, the right side of your vertical divider will be the inverse
of your original matrix
4. Write out the inverse matrix. Copy the elements now appearing on the right side of the vertical
divider as the inverse matrix.
( )
3 0 2
A= 2 0 −2
−1
Example. Find inverse matrix A , if
0 1 1
( ) ( )
3 0 2 1 0 0 5 0 0 1 1 0
' '
2 0 −2 | 0 1 0 →R 1=R 1 +R2 → 2 0 −2 | 0 1 0 →R 1=R 1÷5 →
0 1 1 0 0 1 0 1 1 0 0 1
( ) ( )
1 0 0 0. 2 0 . 2 0 1 0 0 0. 2 0.2 0
' '
2 0 −2 | 0 1 0 → R2 =R3 → 0 1 1| 0 0 1 → R3 =R1⋅(−2 ) +R3 →
0 1 1 0 0 1 2 0 −2 0 1 0
( ) ( )
1 0 0 0.2 0.2 0 1 0 0 0 .2 0. 2 0
' '
0 1 1| 0 0 1 →R 3 =R3 ÷(−2 ) → 0 1 1 | 0 0 1 →R 2=R 3⋅(−1 )+R 2 →
0 0 −2 −0. 4 0 . 6 0 0 0 1 0 .2 −0 . 3 0
( )
1 0 0 0 .2 0. 2 0
'
0 1 1| 0 0 1 →R 2=R 3⋅(−1 )+R 2 →
0 0 1 0 .2 −0 . 3 0
( )
1 0 0 0.2 0.2 0
0 1 0 |−0.2 0.3 1
0 0 1 0.2 −0.3 0
( )
0 .2 0 .2 0
−1
A = −0 .2 0. 3 1
So,
0 .2 −0 .3 0
Rank of a Matrix
Given the rectangular matrix
( )
a11 a12 . .. a 1n
A= a21 a 22 . .. a 2n
. . . . .. . .. ...
am 1 am 2 . .. amn
( )
2 4 3 5
4 8 6 10
6 12 9 15
Solution. All the second - and third-order minors of the given matrix are equal to zero since
the elements of the rows of these minors are proportional, while the first order minors (the elements
of the matrix themselves) are nonzero. Consequently, the rank of the matrix is equal to unity:
r ( A )=1 .
Example №2 Determine the rank of the matrix
( )
1 0 0 0 5
0 0 0 0 0
2 0 0 0 11
Solution. Deleting the 2nd row and then 2nd, 3rd and 4th columns from the matrix, we
1 5
| |=1≠0
Since 2 11 the rank of the given matrix is equal to 2.
( )
3 5 7
1 2 3
Example №3 Determine the rank of the matrix
1 3 5 .
Solution. Add up the elements of the 1st and the 3rd row respectively
( )( )
3 5 7 4 8 12
1 2 3 ~ 1 2 3
1 3 5 1 3 5
Then divide by the 4 the elements of the first row
( )( )
3 5 7 1 2 3
1 2 3 ~1 2 3
1 3 5 1 3 5
From elements of the 1st row subtract the corresponding elements of the 2nd row
( )( )(
1 5 7 0 0 0
1 2 3 ~1 2 3 ~
1 3 5 1 3 5
1 2 3
1 3 5 )
We deleted the first row.
1 2
| |≠0
The rank of the last matrix is equal to 2 since, for instance 1 3 .
Consequently, r ( A )=2 .
( )
4 3 2 2
0 2 1 1
Example №4 Determine the rank of the matrix
0 0 3 3 .
Solution. Subtract the elements of the 3rd column from the elements of the 4th column
( )( )( )
4 3 2 2 4 3 2 0 4 3 2
0 2 1 1 ~ 0 2 1 0 ~ 0 2 1
0 0 3 3 0 0 3 0 0 0 3
4 3 2
|0 2 1|=24≠0
Since 0 0 3 the rank of the matrix is equal to 3 r ( A )=3 .
TEST
Inverse matrix and rank of matrix
A)
( )
5 2
2 1 (0 1 )
1 0
B) C)
(
−5 2
2 −1 ) D)
(5 −2
−2 1 )
E)
( )
1 5
2 2
( )
1 2 −3
A= 0 1 2
2) Given the matrix
0 0 1 . Find its inverse
( ) ( ) ( )
1 −2 7 1 0 0 1 0 3
0 1 −2 2 1 0 0 1 2
A)
0 0 1 B)
−3 2 1 C)
2 0 1
( ) ( )
1 2 7 1 −2 7
0 1 2 2 1 −2
D)
0 0 1 E)
2 0 1
( )
3 2 2
A= 1 3 1
3) Given the matrix
5 3 4 . Find its inverse
( ) ( ) ( )
9 −2 −4 9 2 4 1 2 7
1 1 1
1 2 −1 1 2 1 0 1 2
5 5 4
A)
−12 1 7 B)
12 1 7 C)
5 1 7
( ) ( )
9 −2 −4 0 0 1
1
1 4 −1 1 0 0
4
D)
15 1 7 E)
0 1 0
( 35 −2
4) Solve the matrix equation−4 −5 6 )
X) =(
−1 2
A)
(19 25 ) (35 −4
B)
−2
) (12 −13 ) (00 05 ) (35 24 )
C) D) E)
( )
1 2 3 4
2 4 6 8
6) Determine the rank of the matrix
3 6 9 12
A) 1 B) 2 C) 3 D)4 E) 0
( )
5 0 0 0 4
0 0 0 0 0
7) Determine the rank of the matrix
2 0 0 0 3
A) 0 B) 1 C) 2 D)3 E) 7
( )
1 0 2 1
0 2 4 2
8) Find the rank of the matrix
0 2 2 1
A) 1 B) 2 C) 3 D)4 E) 0
( )
1 2 1 −1
9 5 2 2
9) Find the rank of the matrix
7 1 0 4
A) 0 B) 1 C) 3 D)2 E) 5
( )
3 1 0
2 0 5
10) Find the rank of the matrix
1 2 1
A) 3 B) 0 C) 2 D)5 E) 4
Lecture 4
System of linear equations. General case.
a 1 x 1 + a2 x 2 +. . .+ an x n =b , (1)
{
a11 x 1 + a12 x 2 +. . .+ a1 n x n=b1
a21 x 1 + a22 x 2 + .. .+a2 n x n =b2
⋮
am 1 x 1 + am 2 x 2 +. ..+ amn x n=bm
(2)
The set of all solutions of a linear system is called the solution set of the system.
Theorem 1. Any system of linear equations has one of the following exclusive conclusions
(a) No solution
(b) Unique solution
(c) Infinitely many solutions
A linear system is said to be consistent (or compatible) if it has at least one solution; and it
said to be inconsistent (or incompatible) if it has no solution
Definition. The augmented matrix of the general linear system (2) is the table
( )
a11 a12 . .. a 1n b1
Ā= a21 a 22 . .. a 2n b2
. . . . .. . .. ... ...
am 1 am 2 . .. amn bn
( )
a11 a12 . .. a 1n
a a 22 . .. a 2n
A= 21
. . . . .. . .. ...
am 1 am 2 . .. amn
and the coefficient matrix of (2) is .
Kronecker-Capelli Theorem. A system of linear equations (2) is consistent if and only if
the rank of the coefficient matrix A is equal to that of the augmented matrix Ā , obtained from A
by adding the column of free terms
b i rank ( A )=rank ( Ā ) .
For use this theorem, we compute first the rank of the matrix A , suppose rank A=r . The
non-zero minor which gives us the rank is called the principal minor. Then we compute the
characteristic minors.
Let system (2) is consistent. We suppose that rank A=rank Ā=r and the principal
minor is of the form:
{
a 11 x 1 +a12 x 2 +.. .+a1 n x n =b1
a21 x1 +a22 x2 +.. .+a 2n x n =b 2
.. .. . .. .. . .. .. . .. .. . .. .. . .. .. .. . .. .. . .. .
ar 1 x 1 +a r2 x 2 +.. .+arn x n=br
(3)
The matrix of the coefficients for the system (3) has a non-zero minor of order r (formed by
first r columns), then has rank equal with r , where r ≤n .
2) If r < n , supposing that the minor formed with the coefficients of first
r unknowns is principal, we called the unknowns 1 2
x , x ,. .. , x
r principal unknowns and the
remaining unknowns are called secondary unknowns. We put in the right side all terms which has
secondary unknowns
x r+1 , .. . , x n , we give then arbitrary values λ r+1 , .. . , λ n and we find a system
with r equations and r unknowns 1 2
x , x ,. .. , x r :
{
a11 x 1 +a12 x 2 +. . .+a1 r x r =b1 −a1r +1 λ r+1 −.. .−a1n λn
a21 x1 +a22 x2 +.. .+a 2n x n =b2 −a 2r +1 λr+1 −.. .−a2 n λn
.. .. . .. .. . .. .. . .. .. . .. .. . .. .. .. . .. .. . .. .
ar 1 x 1 +ar 2 x 2 +. ..+a rn x n =br −arr +1 λ r+1 −.. .−arn λ n
(4)
{
a11 x 1 +a12 x 2 +...+a1 n x n=0
a21 x 1 +a22 x 2 +...+a2 n x n =0
........................................
a m1 x 1 +am2 x 2 +...+amn x n =0
(5)
we have the following properties for the homogeneous systems:
1) A homogeneous system is always consistent (we can see this if we replace
x 1 , .. . , x n with zero)
2) Supposing that rank A is r
i) if r =n , then the zero solution is the unique solution of system (5)
ii) if r < n , then the system (5) has non-zero solutions.
A homogeneous system with the same number of equations and unknowns has non-
zero solutions if and only if its determinant is zero.
For a homogeneous system if the number of equations is smalls than the number of
unknowns then the system has non-zero solutions.
Matrix Representation of a Linear System
Matrices are helpful in rewriting a linear system in a very simple form. The algebraic
properties of matrices may then be used to solve systems. First, consider the linear system
{
a 11 x 1 + a12 x 2 + a13 x3 =b 1
a21 x1 + a22 x2 +a23 x 3=b2
a31 x 1 + a32 x 2 + a33 x 3 =b 3
(1)
Set the matrices
( ) () ()
a11 a12 a 13 x1 b1
A= a21 a22 a 23 X = x2 B= b 2
a31 a32 a 33 x3 b3
, and
Using the matrix multiplications, we can rewrite the system (1) as the matrix equation
A⋅X= B
−1
Now, if the above A has an inverse, then both sides can be left-multiplied by A , to get
I 3 X =A−1 B
X =A−1 B
For example
{ ( ) () ( )
x +2 y=10 1 2 0 x 10
3 x+2 y +z=23 A= 3 2 1 X= y B= 23
y+2 z=13 0 1 2 , z , 13
( )
c 11 c 21 c 31
−11
A = c 12 c 22 c 32
A
c 13 c 23 c 33
First, find inverse
1 2 0
Δ=|3 2 1 |=4−1−12=−9
0 1 2
2 1 2 0 2 0
c 11 =| |=3 c 21=−| |=−4 c 31=| |=2
1 2 1 2 2 1
3 1 1 0 1 0
c 12=−| |=−6 c 22=| |=2 c 32=−| |=−1
0 2 0 2 3 1
3 2 1 2 1 2
c 13=| |=3 c 23=| |=−1 c 33=| |=2−6=−4
0 1 0 1 3 2
( )
1 4 2
− −
3 9 9
( )
3 −4 2
1 2 2 1
A−1 =− −6 2 −1 = − −
9 3 9 9
3 −1 −4
1 1 4
−
3 9 9
( )
1 4 2
− −
3 9 9
( )()
10 4
2 2 1
X=
3
−
9 9
23 = 3 x=4
13 5 y=3
1 1 4
−
, i.e. z=5
3 9 9
Lecture 5
Cramer’s rule
Cramer’s rule is an explicit formula for the solution of a system of linear equations
with as many equations as unknowns, valid whenever the system has a unique solution. It
expresses the solution in terms of determinants of the coefficient matrix and of matrices
obtained from it by replacing one column by the vector of right hand sides of the equations.
{
a 11 x + a12 y +a 13 z=b 1
a21 x+ a22 y + a23 z=b2
a31 x +a 32 y +a 33 z=b 3
(1)
a11 a12 a 13
Δ=|a21 a22 a 23|
a31 a32 a 33 is determinant of the coefficient matrix.
Let multiply the first of equations (1) by the cofactor C 11 of entry a 11 , the second one by C 21
and the third equation by
C 31 and them sum together the results. Doing this we derive
But the coefficient of x will become Δ , while the coefficients of all other unknowns
become zero by property 10 in lecture 1. The lefthand side becomes simply Δ⋅x . The right hand
b1 a12 a13
|b 2 a22 a23|
C b +C b +C b
side is 11 1 21 2 31 3 which is equal to b 3 a32 a33 .
b1 a12 a13
Δ x =|b2 a22 a23 |
Thus b 3 a32 a33 is determinant obtained by substituting the constant column for the
coefficient of x .
Δx
So
Δ⋅x=Δ x . Let suppose that Δ≠0 , we find x= Δ .
The formulas
Δx Δy Δ
x= y= z= z
Δ , Δ , Δ are called Cramer’s rule.
For example
{
x − y+ z=5 1 −1 1
2 x + y +z=6 Δ=|2 1 1|=5≠0
x + y +2 z=4 1 1 2
5 −1 1 1 5 1
Δx =|6 1 1 |=15 Δ y =|2 6 1 |=−5
4 1 2 1 4 2
1 −1 5
Δz =|2 1 6 |=5 15 −5 5
x= =3 y= =−1 z= =1
1 1 4 , so 5 , 5 , 5 .
Gauss’ method (the method of elimination)
Consider a linear system
{
a 11 x 1 + a12 x 2 + a13 x3 =b 1
a21 x1 + a22 x2 +a23 x 3=b2
a31 x 1 + a32 x 2 + a33 x 3 =b 3
(1)
{
a12 a b
x 1+ x 2 + 13 x 3 = 1
a11 a 11 a11
a21 x1 +a22 x2 +a23 x 3=b2
a31 x 1 +a32 x 2 +a33 x 3 =b 3
(2)
−a 21 (−a31 )
Multiply the first equation of (2) by and add it to the second (third) equation of system
(2). This eliminates x 1 and yields the system
{
x 1 +a'12 x 2 +a'13 x 3=b'1
a'22 x 2 +a '23 x 3 =b'2
a'32 x2 +a'33 x 3=b'3
(3)
(the primes designate the new coefficients here)
Now we divide the second equation of (3) by a 22 , then multiply by (−a32 ) and add with
'
{
x 1 +b12 x2 +b13 x 3=B1
x 2 + b23 x 3=B2
x 3 =B3
Example
{ {
2 x 1 + x 2 −x 3=1 x1 −x 2 +2 x 3 =5
3 x1 + 2 x 2 −2 x3 =1 2 x 1 + x 2 −x 3 =1
x 1 −x 2 +2 x 3 =5 ⇒ 3 x 2 +2 x 2 −2 x3 =1 ⇒
{ {
x −x +2 x =5
x1 −x 2 +2 x 3 =5 1 2 3
5
3 x 2−5 x 3 =−9 x 2 − x 3 =−3
3
⇒ 5 x 2 −8 x 3 =−14 ⇒ 5 x 2 −8 x 3 =−14 ⇒
{
x1 −x 2 +2 x 3 =5
5 5
x 2 − x 3 =−3 x 2 − ⋅3=−3
3 3
1 x 2 =2
x 3 =1
⇒ 3 ⇒ x 3 =3 , x 1=1
TEST
A) 0
{x−3y+z=−6¿{5x+y−2z=4¿ ¿
B) 6 C) 7 D) 8 E) 9
A) 5
{x+3y−z=4¿{2x−4y−z=−7¿ ¿
)5 B) 6 C) 7 D) 8 E) 9
{5x+y−3z=3¿{2x−3y+7z=19¿ ¿
B) 5 C) 6 D) 7 E) 8
A) 4
Lecture 6 (LA)
Linear Spaces
Basic notions. Let us consider a set R of elements x , y , z, . .. such that the sum x+ y∈ R is defined
for any two elements x ∈ R and y ∈R and the product λx ∈ R is defined for any element x ∈ R
and any real number λ .
If summation of the members of the set R and multiplication of a member of that set by a real
number satisfy the conditions:
1. x + y= y + x
2. ( x+ y )+ z=x + ( y + z )
4. there exists an element y ∈R for every element x ∈ R such x + y=0 ( in what follows we`ll
write y=− x i.e. x + (−x )=0 )
5. 1⋅x=x
6. λ ( μx )=( λμ ) x
7. ( λ+ μ ) x= λx+ μx
8. λ ( x + y )= λx+ λy
then the set R is said to be linear(or vector) space, and the elements x , y , z, . .. of that space are
called vectors.
The difference of two vectors ⃗
x and ⃗y of a linear space is a vector ⃗v of that space such that
⃗y +⃗v =⃗x . The difference of two vectors ⃗x and ⃗y is designated as ⃗x −⃗y ,i.e ⃗x −⃗y =⃗v
If equality (1) can be also satisfied in the case when not the numbers α , β , γ , . .. , λ are equal to zero,
the vectors ⃗x ,⃗y ,⃗z ,...⃗
u are said to be linearly dependent.
Dimensionality and the basis of a linear space. If there are n linearly independent vectors in the
linear space R but any n+1 vectors of that space are linearly dependent, the space R is said to be n-
dimensional. It is also customary to say that the dimensionality of space R is equal to n and write
d ( R )=n . The space containing infinitely many linearly independent vectors is called infinite-
dimensional. If R is infinite-dimensional space, thend ( R )=∞ .
A collection of n linearly independent vectors of an n-dimensional linear space is known as the
basis. The following theorem hold true: every vector of a linear n-dimensional space can be
uniquely represented as a linear combination of the vectors of the basis. Thus, if 1 2
⃗e , ⃗e ,.. .,⃗e
n is the
basis of an n-dimensional linear space R, then any vector x ∈ R can be uniquely represented in the
form
⃗x =ξ 1 ⃗e 1 + ξ2 ⃗e 2 +. . .+ ξ n ⃗e n
⃗x =ξ 1 ⃗e 1 + ξ2 ⃗e 2 +. . .+ ξ n ⃗e n ⃗y =η1 ⃗e 1 +η2 ⃗e 2 +. . .+ ηn ⃗e n
If , , then
x + y=( ξ1 +η1 ) e1 +( ξ 2 +η 2 ) e 2 +. ..+( ξ n +η n ) e n
{
α 1 ξ11+ α 2 ξ 12+. . .+ α n ξ1 n = 0
α 1 ξ 21 + α 2 ξ 22+. . .+ α n ξ3 n = 0
.. . ... . .. .. .
α 1 ξ n1 + α 2 ξ n 2 +. ..+ α n ξ nn= 0
ξ 11 ξ 12 . .. ξ1 n
|ξ 21 ξ 22 . .. ξ2 n |≠0
. .. . .. . .. ...
ξn 1 ξn 2 . .. ξ nn
x 1=( ξ 11 ; ξ 21 ) x 2 =( ξ 12 ; ξ 22 )
In particular , are linearly independent if and only if
,
ξ 11 ξ 22−ξ 12 ξ 21≠0
x 1=( ξ 11 ; ξ 21)
Example 2. Given a linear space of various pairs of the ordered real numbers ,
x 2 =( ξ 12 ; ξ 22 ) x 3 =( ξ 13 ; ξ23 )
, …,and addition of vectors and multiplication of a vector by a real
x i + x k =( ξ 1i + ξ 1k ; ξ 2 i +ξ 2k ) λx i= ( λξ 1i ; λξ 2 i )
number are determined by the equalities :
It is easy to see that there exists a unique pair of values ( λ , μ ) for which this equality is satisfied.
This follows from the fact that the system of equations
{ λ+ 3 μ=
2 λ+4 μ=
η1
η2
is determinate.
e e
Thus, the vectors 1 and 2 form a basis. Determine the coordinates of the vector x=( 7 ;10 ) in this
basis. The problem reduces to determining λ and μ from the system of equations
{ λ+3 μ= 7
2 λ+4 μ= 10
λ=1 , μ=2 i.e x=e 1 +2 e 2
This yields
Linear Transformations
Basic notions. A transformation A is said to be given in the linear space R if in accordance with
some rule each vector x ∈ R is put into correspondence with the vector Ax ∈ R . A transformation
A is called linear if the equalities
A ( x + y )= Ax+ Ay A ( λx )=λ Ax
,
are satisfied for any two vectors x and y and for any real number λ .
itself. An
A linear transformation is called identity transformation if it maps any vector x into
identity linear transformation is denoted by I. Thus, Ix=x .
( )
a11 a12 .. . a1n
A= a 21 a22 .. . a 2 n
... .. . .. . .. .
an 1 a n2 .. . a nn
Since
Ax ∈ R , the vector Ax can be also resolved with respect to the vectors of the basis
1. A ( 0 ) =0
2. A (−x )=− A ( x )
3. A ( y−x )= A ( y )− A ( x )
4. If
x=α 1 x 1 +α 2 x 2 +. ..+α n x n then
A ( x )= A ( α 1 x 1 +α 2 x 2 +. ..+ α n x n )=α 1 A ( x 1 ) +α 2 A ( x 2 ) +. ..+ α n A ( x n )
Proof
To prove the first property, note that 0 x=0 . Then it follows that A ( 0 ) =A ( 0 x ) =0 A ( x )=0 .
The third property follows from y−x = y+ (−x ) , which implies that
A ( x + y )= Ax+ Ay A ( λx )=λ Ax
The fourth property follows from , ,
A ( x )= A ( α 1 x 1 +α 2 x 2 +. ..+ α n x n )= A ( α 1 x 1 ) + A ( α 2 x 2 ) +. ..+ A ( α n x n ) =
α 1 A ( x 1 ) + α 2 A ( x 2 ) +.. .+α n A ( x n )
which implies that
3 3
Example 4. Let A : R → R be a linear transformation such that
Suppose R is a given n-dimensional linear space. The nonzero vector x ∈ R is called an eigenvector
of the linear transformation A if there is a number λ such that the equality Ax=λx .The number
λ is called a characteristic number ( eigenvalue) of the linear transformation A corresponding to
the vector x .
Eigenvalues and eigenvectors can be found as follows. Since x=Ix we can rewrite equation
Ax=λx in the form ( A−λI )x =0 .
( )
a11 a12 .. . a1n
a a22 .. . a 2 n
A= 21
... .. . .. . .. .
an 1 a n2 .. . a nn
then the characteristic numbers of the linear transformation A are the real roots
λ 1 , λ2 , .. . , λn of the
nth degree equation which can be written in the form
a11− λ a 12 . .. a1 n
a a 22−λ . .. a2 n
| 21 |=0
. .. . .. . .. . . .
an 1 an 2 . .. ann−λ
characteristic equation and its left-hand side, as a characteristic polynomial of the
It is known as
linear transformation A. The eigenvector
⃗x k corresponding to characteristic number λ k is any vector
ξ 1 ⃗e 1 +ξ 2 ē 2 +. . .+ ξ n ⃗e n whose coordinates
satisfy the system of homogeneous equations
{
( a 11−λ k ) ξ1 + a12 ξ 2 +.. .+ a1 n ξ n = 0
a21 ξ1 + ( a 22−λk ) ξ 2 +. ..+ a2 n ξ n = 0
.. . .. . . .. ...
a n1 ξ 1 + a n2 ξ2 + .. .+ ( ann− λk ) ξ n= 0
If the matrix A of the linear transformation A is symmetrical, then all the roots of the
characteristic equation|A−λI|=0 are real numbers.
Example 5. Find characteristic numbers and eigenvectors of the linear transformation A specified
' '
by the equations x =5 x+ 4 y , y =8 x +9 y
Solution. The transformation matrix is written as follows:
( )
A= 5 4
8 9
The characteristic equation has the form
5−λ 4
| |=0 2
8 9−λ or λ −14 λ+13=0
{( 5− λ1 ) ξ1 +
8 ξ1 +
4 ξ 2=0
( 9−λ1 ) ξ2 =0 { ( 5− λ2 ) ξ1 +
8 ξ1 +
4 ξ 2 =0
( 9−λ2 ) ξ 2=0
Since
λ 1=1 , the first system can be written as follows:
{ 4 ξ1 +
8 ξ 1+
4 ξ2 =0
8 ξ 2 =0
The value
λ 2=13 leads to a system of equations
{−8 ξ 1 +
8 ξ1 −
4 ξ 2=0
4 ξ 2=0
Euclidean Space
The linear space is said to be Euclidean if there is rule which makes it possible to construct, for
x and ⃗
⃗ y belonging to R, a real number called a scalar product of the vectors⃗x
every two vectors
and ⃗y and designated as ⃗x⋅⃗y , the rule complying with
the following conditions:
⃗x⋅⃗y =⃗y⋅⃗x
x⋅( ⃗
⃗ y +⃗z )=⃗x⋅⃗ y +⃗ x⋅⃗z
( λ ⃗x )⋅⃗y=λ ( ⃗x⋅⃗y ) for any real number λ
⃗
x⋅⃗
x >0 if ⃗
x ≠0
It follows from these conditions that
(a) ( ⃗y +⃗z )⋅⃗x =⃗y⋅⃗x +⃗z⋅⃗x
(b) ⃗
x⋅( λ ⃗y )= λ(⃗x⋅⃗y )
x =0 for any vector ⃗
(c) 0⋅⃗ x
a ij are some numbers, among which at least one is nonzero. This sum is called quadratic form
where
of the real variables
x 1 , x2 ,. .. , x n is denoted F ( x1 , x 2 , .. . , x n )
n n
F ( x 1 , x 2 , .. . , x n )=∑ ∑ a ij x i x j
Thus, i=1 j−1 (1)
a
The numbers ij are called the coefficients of quadratic form. The quadratic form can be written in
xx x x x x =x j x i ,
such a way that the coefficients at i j and j i are equal to each other. Indeed, since i j
then
1 1
a ij x i x j +a ji x j x i=( aij + a ji ) x i x j= ( aij + a ji ) x i x j + ( aij +a ji ) x j xi
2 2
In what follows we will assume that in quadratic form (1)
a ij=a ji (2)
number, then
F ( λx1 , λx 2 , . .. , λx n )=λ 2 F ( x 1 , x2 ,. .. , x n )
( )
a11 a12 .. . a1n
a a22 .. . a 2 n
A= 21
... .. . .. . .. .
an 1 a n2 .. . a nn
By virtue of equality (2) the matrix A has equal elements located symmetrically to the principal
diagonal. Such matrix is called symmetric matrix. The symmetric matrix
( )
a11 a12 .. . a1n
a a .. . a 2 n
A= 12 22
. . . .. . .. . .. .
a1n a 2n .. . a nn
Quadratic form
n n
∑ ∑ aij x i x j a il=0 for i≠ j .Hence ,canonical quadratic form is
i=1 j−1 is called canonical if
n
a 11 x 21 + a22 x22 +.. .+a nn x 2n =∑ aii x 2i
i=1 and its matrix is diagonal.
If n=2 then quadratic form is
F ( x 1 , x 2 ) =a 11 x 21 + 2 a12 x 1 x 2 +a 22 x 22
( )
a11 a12 a13
A= a12 a22 a23
a13 a23 a33
()
y1
¿ y 1 x 1 + y 2 x 2 + y 3 x 3 =( x 1 x2 x3 ) y2
y3
where
( )( )( )( ) ( )
y1 a11 x1 +a12 x2 +a13 x 3 a11 a 12 a13 x 1 x1
y 2 = a 12 x 1 +a 22 x 2 +a 23 x 3 = a12 a 22 a23 ⋅ x 2 =A x 2
y3 a13 x1 +a23 x 2 +a 33 x 3 a13 a 23 a33 x3 x3
()
x1
X = x2
x3
Thus, if we introduce the number vector we obtain
F=X T AX (4)
Conversely, if a form is represented as (4) and if the matrix A is symmetric then A is the matrix of
this quadratic form.
Let us now perform an arbitrary linear transformation of the variables of form
'
This transformation can be put down in the matrix form as X =B X (5)
T T T
Then by formula ( AB ) =B A
T T ¿
(6) , we have X = X B which implies
F=X BT AB { X = X ( BT AB ) X ¿
¿ ' ¿ '
diagonal elements are the eigenvalues of the matrix A . We obtain the quadratic form
' ' ' ¿ ¿ ¿
F ( x1 , x 2 , . x 3 )= λ1 x 1 + λ 2 x2 + λ3 x 3 which does not contain terms with products
x '1 x '2 , x '1 x '3 , x '2 x'3 .
It is customary to say that quadratic form F ( x1 , x 2 , . x 3 ) has been reduced to the canonical form by
means of the orthogonal transformation B. The argument was carried out under the assumption that
eigenvalues
λ 1 , λ2 , λ3 were distinct
F ( x 1 , x 2 , x 3 ) =( 2 x 1 +4 x 2 + x 3 ) x 1 + ( 4 x 1−3 x 2 + 0 ) x 2 + ( x 1 + 0+ x 3 ) x 3
Solution .
Matrix of this quadratic form is
( )
2 4 1
A= 4 −3 0
1 0 1
( )( )
2 4 1 x1
F ( x 1 , x 2 , x 3 ) =( x 1 x2 x3 ) 4 −3 0 x 2
1 0 1 x3
Hence,
Quadratic form
F ( x 1 , x 2 ) =2 x21 +4 x 1 x 2 −3 x22 is given
Example 2. .
Find quadratic form obtained from given form by transforming
{ x1 =2 y 1 −3 y 2
x 2= y 1+ y 2
A= 2 2(
2 −3 )
and matrix of transformation is
B= ( 2 −3
1 1 )
Hense, matrix of the required quadratic form is
F=BT AB= ( )(
2 1 2 2 2 −3 13 −17
−3 1 2 −3 1 1
=
−17 3 )( )( )
F ( y 1 , y 2 )=13 y 21 −34 y1 y 2 +3 x 22
And quadratic form have form
Definiteness of Quadratic Forms
The quadratic formQ ( x , y )=x + y is positive for all nonzero (that is( x , y ) ≠( 0 , 0 ) ) arguments
2 2
Let
(
Q ( x 1 , x2 )=ax 21 +2 bx 1 x2 +cx 22 =( x 1 , x 2 )⋅
a
b
b x1
⋅
c x2)( )
be a 2 variable quadratic form.
Here
A= ( )
a b
b c is the symmetric matrix of the quadratic form.
a b
A=| |=ac−b2
The determinant b c is called discriminant of Q.
Easy to see that
ax 12 +2 bx 1 x 2 + cx 22 =a ( )
b 2 ac−b2 2
x 1+ x 2 +
a a
x2
.
D '1 =c
(of degree 1) .
Then
b 2 D
(
Q ( x 1 , x2 )=D 1 x 1 + x2 + 2 x 22
a D1 )
From this expression we obtain:
1. If
D 1 > 0 and D 2 > 0 then the form is of x 2 + y 2 type, so it is positive definite;
2. If
D 1 < 0 and D 2 > 0 then the form is of −x 2 − y 2 type, so it is negative definite;
3. If
D 1 > 0 and D 2 < 0 then the form is of x 2 − y 2 type, so it is indefinite;
If
D 1 < 0 and D 2 < 0 then the form is of−x 2 + y 2 type, so it is also indefinite;
Thus if
D 2 < 0 then the form is indefinite.
Q ( x 1 , x2 , x 3 )=x 21 +2 x 22 −7 x23 −4 x 1 x 2 +8 x 1 x 3
( )
1 −2 4
A= −2 2 0
The symmetric matrix of this quadratic form is
4 0 −7
The leading principal minors of this matrix are
1 −2 4
1 −2 |D |=|−2 2 0 |=−18
|D2|=| |=−2 3
|D1|=|1|=1 , −2 2 4 0 −7
Now look:
Q ( x 1 , x2 , x 3 )=x 21 +2 x 22 −7 x32−4 x 1 x 2 +8 x 1 x3 =
x 21 −4 x 1 x 2 +8 x 1 x 3 +2 x 22−7 x 23 =x 21 −4 x 1 ( x 2 −2 x 3 ) +2 x 22−7 x 23 =
[ x 21−4 x1 ( x 2−2 x 3) + 4 ( x 2−2 x 3)2−4 ( x 2−2 x 3 )2] + 2 x 22−723=
2
[ x 1 −2 x 2 + 4 x 3 ] −2 x 22 +16 x 2 x 3−23 x 23 =
[ x 1 −2 x 2 + 4 x 3 ]2−2 ( x 22 −8 x 2 x 3 )−23 x 23=
[ x 1 −2 x 2 + 4 x 3 ] −2 [ x 22 −8 x 2 x 3 +16 x 23−16 x 23 ]−23 x 23
2
()( )( )
l1 1 −2 4 x1
l2 = 0 1 −4 ⋅ x2
l3 0 0 1 x3
( )
1 −2 4
P= 0 1 −4 ⋅¿ ¿
where
0 0 1
is a nonsingular matrix (changing variables).
Now turn to general 3 variable quadratic form
( )( )
a11 a 12 a13 x 1
Q ( x 1 , x2 , x 3 )=( x 1 , x 2 , x 3 ) a21 a 22 a23 ⋅ x 2
a31 a 32 a33 x 3
a11 a 12 a13
a a 12 |D3|=|a21 a 22 a23|
|D2|=| 11 |
|D1|=|a1|, a21 a22 , a31 a 32 a33
2 |D 2| 2 |D 3| 2
Q ( x 1 , x2 , x 3 )=|D 1|l 1 + l+ l
|D1| 2 |D2| 3
a11− λ a 12 a13
| a 21 a 22−λ a23 |=0
a 31 a 32 a33−λ are real numbers.
Assume that
e '1 = b 11 e1 + b 21 e 2 + b 31 e 3
'
e 2 = b12 e1 + b 22 e 2 + b 32 e 3
e '3 = b13 e 1 + b23 e 2 + b 33 e 3
( )
b 11 b 12 b13
B= b21 a 22 a23
a31 a 32 a33
Solution . Here
a 11=27 ,a 12=−5 ,a 22=3 . Derive a characteristic equation
27−λ ¿5
| |=0
−5 3−λ
0r λ 2−30 λ +56=0
That is, the eigenvalues
λ 1=2 , λ 1=28 .
{25 ξ1 −5 ξ 2 =
−5 ξ 1 +ξ 2 =
0
0
Thus we have
ξ 2 =5 ξ1 .Putting ξ 1 =c ,we obtain ξ 2 =5 c ,that is,
⃗ =c ( ⃗
u e 2)
e 1 +5 ⃗
.
If
λ 1=28 ,we arrive at a system
{ −ξ 1−5 ξ 2 = 0
−5 ξ 1−25 ξ 2 = 0
v =c (−5 ⃗
⃗ e2 )
e 1 +⃗
In this case we obtain an eigenvector
1 1
c= =
To normalize the vectors ⃗
u ,⃗v , we must assume √ 12 +52 √ 26 . Thus we have found the
' 1 ' 1
e 1= ( ⃗e 1 +5 ⃗e 2 ) e 2 = (−5 ⃗e 1 +⃗e 2 )
normalized eigenvectors √ 26 , √ 26 .
( )
1 −5
B= √ √ 26
26
5 1
√26 √ 26
1 ' 5 '
x 1= x 1− x
√ 26 √26 2
5 ' 1 '
x 2= x 1+ x
√ 26 √ 26 2
It follows that
'
(√26 )
2
5 1
+3 x'1 + x'2 =2 x ' 2 +28 { x
√ 26 1
22
' '
F= λ 1 x 2+ λ 2 x 2
This result could be obtained directly since 1 2
Lecture 8(LA)
Coordinate systems
The rectangular coordinate system consists of two real number lines that intersect at a right angle.
The horizontal number line is called the x-axis, and the vertical number line is called the y-axis.
These two number lines define a flat surface called a plane, and each point on this plane is
associated with an ordered pair of real numbers (x, y). The first number is called the x-coordinate,
and the second number is called the y-coordinate. The intersection of the two axes is known as
the origin, which corresponds to the point (0, 0).
An ordered pair (x, y) represents the position of a point relative to the origin. The x-coordinate
represents a position to the right of the origin if it is positive and to the left of the origin if it is
negative. The y-coordinate represents a position above the origin if it is positive and below the
origin if it is negative. Using this system, every position (point) in the plane is uniquely identified.
For example, the pair (2, 3) denotes the position relative to the origin as shown:
This system is often called the Cartesian coordinate system, named after the French mathematician
René Descartes (1596–1650).
The x- and y-axes break the plane into four regions called quadrants, named using roman numerals I,
II, III, and IV, as pictured. In quadrant I, both coordinates are positive. In quadrant II, the x-
coordinate is negative and the y-coordinate is positive. In quadrant III, both coordinates are negative.
In quadrant IV, the x-coordinate is positive and the y-coordinate is negative.
M1 P M 1 M x−x 1
= =λ =λ
MQ MM 2 x
i.e 2
−x and
x−x 1 =λx 2− λx .From the last relations we
M2
deduce
x 1 + λx 2 y 1 + λy 2 M
x= y= Q
1+ λ , 1+ λ
M1
P
where (x, y) are old coordinates [i.e. coordinates relative to xy system], (x',y') are new coordinates
[relative to x'y' system] and (x0, y0) are the coordinates of the new origin 0' relative to the
old xy coordinate system.
where the origins of the old [xy] and new [x'y'] coordinate systems are the same but the x' axis
makes an angle α with the positive x axis.
where the new origin O' of x'y' coordinate system has coordinates (x0, y0) relative to the
old xy coordinate system and the x' axis makes an angle α with the positive x axis
Area of a triangle
x y 1
1 1 1 1 x −x y − y
A=(± ) | x2 y 2 1 |=(± ) | 2 1 2 1 |
2 2 x 3 −x 1 y 3 − y 1
x3 y 3 1
where ( 1 1 ) , ( 2 2 ) and ( 3 3 ) are the vertices of the triangle
A x ,y B x ,y C x ,y
ABC, and
the sign in the right side is chosen so that the area of the triangle is nonnegative.
Area of a quadrilateral
1
2[ 1 2 1 2
A=(± ) ( x −x )( y + y ) + ( x 2 −x 3 )( y 2 + y 3 ) + ( x 3 −x 0 )( y 3 + y 0 ) + ( x 0 −x 1 )( y 0 + y 1 ) ]
where ( 1 1 ) , ( 2 2 ) , ( 3 3 ) , ( 0 0 ) are the vertices of the
A x ,y B x ,y C x ,y D x ,y
quadrilateral ABCD. The sign in the right side is chosen so that the area of the
quadrilateral is nonnegative.
Polar coordinates
In two dimensions, the Cartesian coordinates (x,y) specify the location of a point P in the plane.
Another two-dimensional coordinate system is polar coordinates. Instead of using the signed
distances along the two coordinate axes, polar coordinates specifies the location of a point P in the
plane by its distance r from the origin and the angle θ made between the line segment from the
origin to P and the positive x-axis. The polar coordinates (r,θ)of a point PP are illustrated in the
below figure.
As r ranges from 0 to infinity and θ ranges from 0 to 2π, the point P specified by the polar
coordinates (r,θ) covers every point in the plane. Adding 2π to θ brings us back to the same point, so
if we allowed θ to range over an interval larger than 2π, each point would have multiple polar
coordinates. Hence, we typically restrict θ to be in the interval 0≤θ<2π. However, even with that
restriction, there still is some non-uniqueness of polar coordinates: when r=0, the point P is at the
origin independent of the value of θ
Conversion formulas
We can calculate the Cartesian coordinates of a point with polar coordinates (r,θ) by forming the
right triangle illustrated in the below figure. The hypotenuse is the line segment from the origin to
the point, and its length is r. The projection of this line segment on the x-axis is the leg of the
triangle adjacent to the angle θ, so x=rcosθ. The y-component is determined by the other leg,
so y=rsinθ .Our conversion formula
is
x=r cos θ , y=r sin θ (1)
.
To go the other direction, one can use the same right triangle. Since r is the distance from the origin
√ 2 2
to (x,y), it is the magnitude r = x + y . Alternatively, from the equation (1), one can calculate
directly that
x 2 + y 2 =r 2 cos 2 θ+r 2 sin2 θ=r 2 ( cos 2 θ+sin 2 θ )=r 2
Taking the ratio of y and x from equation (1), one can obtain a formula for θ
y r sin θ
= =tan θ
x r cosθ
,
One can also see this relations from the above right triangle. We can set
y
θ=arctan
x
but a problem is that arctan gives a value between −π/2and π/2. One might need to add π or 2π to
get the correct angle. With this caveat (and also mapping points where x=0to θ=π/2 or −π/2), one
obtains the following formula to convert from Cartesian to polar coordinates
r =√ x + y
2 2
(2)
y
θ=arctan
x
the formula
d= √( x −x ) +( y − y
2 1
2
2 1
2
) + ( z 2 −z 1 )
2
3.Dividing a line segment in the ratio λ
Suppose that the point ( 0 0 0 ) divides the segment AB in the ratio λ. The coordinates of the
C x , y ,z
point C are given by the expressions
x 1 + λx2 y 1 + λy 2 z 1 + λz 2 AC
x 0= y 0= z 0= λ= , λ≠−1
1+ λ , 1+ λ , 1+ λ , where CB
where (
x1 , y1 , z1)
are the coordinates of the point A, and (
x2 , y2 , z2)
are the coordinates of the
point B.
The coordinates of the midpoint of the line segment are obtained from the previous formulas
4.
at λ=1 and are written as
x 1+ x 2 y 1+ y 2 z 1+ z 2 AC
x 0= y 0= z 0= λ= =1
2 , 2 , 2 , where CB
5.Area of a triangle
A ( x 1 , y 1 , z1 ) B ( x 2 , y 2 , z 2 ) C ( x3 , y 3 , z 3 )
The area of a triangle with the vertices and is found by
, ,
the formula
[ ]
y z1 1 2 z1 x1 1 2 x1 y1 1 2
1 1
A= | y 2 z 2 1| +|z 2 x 2 1| +|x 2 y2 1|
2
y3 z3 1 z 3 x3 1 x3 y3 1
6. Volume of a pyramid
A ( x 1 , y 1 , z1 ) B ( x2 , y2 , z2 )
The volume of a pyramid whose vertices have the coordinates and
, ,
C ( x3 , y 3 , z 3 ) D ( x 4 , y 4 , z 4 )
is determined by the expression
,
x1 y1 z1 1
1 x y2 z2 1 x −x y 1− y 4 z 1−z 4
V =± | 2 | 1 1 4
6 x3 y3 z3 1 V =± |x 2 −x 4 y 2− y 4 z 2−z 4 |
6
x4 y4 z4 1 x 3 −x 4 y 3 − y 4 z 3−z 4
or
The sign in the right side of the formulas above is chosen so that to get a positive value for the area
or volume.
Lecture 9
Vectors and operations on the vectors
A vector is a directed line segment drawn from a point P (called its initial point) to a point
Q (called its terminal point), with P and Q being distinct points. The vector is denoted by ⃗
PQ or
⃗v . Its magnitude is the length of the line segment, denoted by |⃗
PQ|. The
⃗v Q zero vector is just a point, and it is denoted by 0 .
P
Two nonzero vectors
with zero magnitude is equal to the zero vector.
Vector Arithmetic
A scalar is a quantity can be represented by single number
1) For a scalar k and a nonzero vector ⃗v , the scalar multiple of ⃗v by k denoted by k ⃗v , is the
⃗v
1) For a vector
⃗PQ in R2 with initial point P ( x 1 , y 1 ) and terminal point Q ( x 2, y 2 ) , the
⃗
magnitude of PQ is
⃗
|PQ|= √ 2
( x2−x 1) + ( y 2 − y 1 )
2
v is |⃗v|= a +b
2) For a vector ⃗v ( a , b ) ,the magnitude of ⃗ √ 2 2
d= √( x −x ) +( y − y ) +( z −z )
2 1
2
2 1
2
2 1
2
v is |⃗v|=√ a +b +c
2 2 2
4) For a vector ⃗v ( a , b , c ) in R the magnitude of ⃗
3
(a) ⃗
v +w
⃗ =⃗
w +⃗
v commutative law
⃗
There are specific unit vectors which we will often use, called the basis vectors: i =( 1, 0, 0 ) ,
⃗j=( 0, 1, 0 ) , ⃗k=( 0, 0, 1 ) in R3 .
They are mutually perpendicular, since they lie on distinct coordinate axes
|⃗i |=|⃗j|=|⃗k|=1
Every vector can be written as a unique scalar combination of the basis vectors:
v=( a , b )=a ⃗i +b ⃗j in R2
⃗v ( a , b ) ⃗v ( a , b , c )
2
k̄
ī
j̄ y
1 j̄
ī x x
1
When a vector ⃗v =( a ,b , c ) is written as ⃗v =a ⃗i +b ⃗j+c ⃗k , we say that ⃗v is in component
⃗
form, and that a,b,c are ⃗i ; ⃗j ; k components, respectively of ⃗
v.
Let a vector ⃗
AB and an axis l be given. The projection l proj ⃗
of the vector ⃗
AB
AB on the
axis l is the length of the segment A B , connecting the feet of the perpendiculars drawn from the
' '
О А/ В/ u
{
' '
|A B | if the direction of the segment
projl ⃗
AB= ' '
A B coincides with positive direction of the axis
' '
−|A B | in the direction opposite to it
Properties of projections
1)
projl ⃗
AB=|AB|cos ϕ , ϕ is the angle between the vector and the axis
А
В
В
А
О А/ В/ B/ A/
⃗v⋅⃗w =|v||w|cos θ ,
where θ is the angle between ⃗ ⃗.
v and w
Notice that the dot product of two vectors is a scalar, not a vector.
Properties of the dot product
(a) ⃗
v⋅⃗
w=w
⃗⋅⃗v commutative law
(b) ( k ⃗v )⋅w
⃗ =⃗v ( k w
⃗ )=k ( ⃗v⋅⃗
w ) associative law
u ( ⃗v + w
(c) ⃗ ⃗ )=⃗u ⃗v +⃗u w
⃗ distributive law
(d)
⃗v⋅⃗w =|v|proj v ω=|w|proj w v
⃗ ⃗ is
v and w
⃗v⋅⃗
w =v 1 w1 +v 2 w2 +v 3 w 3 .
v 1 w1 +v 2 w 2 +v 3 w 3=0 .
Corollary 2. Let ⃗ ⃗ be nonzero vectors and let θ be the angle between them.
v and w
⃗v⋅⃗w
cosθ=
Then |v|⋅|w| or
v 1 w 1 + v 2 w2 +v 3 w3
cos θ=
√ v 21 + v 22 + v 23⋅√w 21 + w22 + w23 .
The cross product v×w is defined as a vector c̄ that is perpendicular to both v and w
with a direction given by the right-hand rule and a magnitude equal to the area of the parallelogram
that the vectors span.
If ϕ is the angle between nonzero vectors v̄ and w̄ in R3 , then
c̄
|v̄×w̄|=|v||w|sin ϕ .
Properties.
ū , v̄ , w̄ in R3 , we have
1
A= |v̄ ×w̄|
2
(b) The area A of a parallelogram with adjacent sides v̄ and w̄ is
A=|v̄×w̄|
Computing the cross product
Coordinate notation. The standard basis vectors ī , j̄ , k̄ satisfy the following equalities in a right
hand coordinate system:
z ī= j̄×k̄
j̄= k̄× ī
k̄= ī× j̄
k̄
which imply, by the anticommutativity of the cross product, that
ī
j̄ y
x
k̄× j̄=− ī , ī× k̄=− j̄ , j̄× ī =−k̄
The definition of the cross product also implies ī× ī = j̄× j̄=k̄×k̄=0 .
These equalities are sufficient to determine the cross product of any two vectors ū and v̄ .
Let
ū=u1 ī +u2 j̄+u3 k̄ , v̄=v 1 ī +v 2 j̄+v3 k̄ .
+ ( u 1 v 2 −u2 v 1 ) k̄
or
u u u u u u
ū×v̄=| 2 3 | ī−| 1 3 | j̄+| 1 2| k̄
v2 v3 v1 v3 v1 v2
it is cofactor expansion by first raw of determinant
i j k
ū×v̄=|u1 u2 u3 |
v1 v2 v3
Scalar triple product (Mixed product)
The scalar triple product is defined as the dot product of one of the vectors with the cross
product of the other two
ū ( u1 ,u2 ,u 3 ) v̄ ( v 1 , v 2 , v 3 ) w̄=( w1 , w 2 , w3 ) R3
Theorem. For any vectors , , in :
u 1 u2 u3
ū⋅( v̄× w̄ )=| v 1 v 2 v 3 |
w1 w 2 w 3
ī j̄ k̄
ū×v̄=|1 −5 −1 |= ī|−5 −1|− j̄|1 −1 |+k̄|1 −5|=−18 { ī −5 j̄+7 k̄ ¿
−3 −3 2 3 2 −3
2 −3 3
Example 2. Find area of the parallelogram with adjacent sides ū=2 ī +3 j̄+5 k̄ and
v= ī +2 j̄ + k̄
ī j̄ k̄
ū×v̄=|2 3 5 |= ī|3 5 |− j̄|2 5 |+ k̄|2 3|=−7 ī +3 j̄+ k̄
2 1 1 1 1 2
1 2 1
A=|w̄× v̄|=√ 49+9+1=√ 59
Example 3. Find the volume of the parallelepiped with adjacent sides ū ( 2, 1, 3 ) , v̄ (−1, 3, 2 ) ,
w̄ (1, 1, −2 )
2 1 3
u⋅( v×w )=|−1 3 2 |=−28
1 1 −2 .
So V =|−28|=28 .
TEST
3) Find the angle between the vectors ā (−1 ,−2 ) and b̄ (−2,−4 )
A) 00 B) 600 C) 300 D) 450 E) 900
A) 4 B) 5 C) 7 D) 8 E) 9
Consider a line l on plane XOY . Suppose that we know a point M 1 ( x 1 , y 1 ) and vector
⃗n =( A , B ) normal to the line. Let M ( x , y ) be any point on line. Vector M 1 M=( x−x 1 ; y − y 1 ) that
joins the two points and thus lies on the line is perpendicular to ⃗ n . Therefore scalar product
⃗n⋅M 1 M =0
That is
A( x x1 ) B( y y1 ) 0 (1)
Equation (1) is the straight line equation passing through the given point with given vector n ,
normal to the line.
Changing, we get
l y
Ax Ax1 By By1 0
Ax By ( Ax1 By1 ) 0 M n
M
Let’s designate
M1
Ax1 By1 C
Ax By C 0 (2)
Intercept form
Changing the general form, we get
Ax By C Divide each side by C
A B
x y 1
C C
x y
1
C C (0,b)
or A B . Let’s designate
b
(а,0)
а
C C
a, b
A B
We get
x y
1
a b (3)
(3) is the intercept form for equation of a line. A straight line intersects x-axis at ( a,0 ) and y-
axis at ( 0, b )
The slope – intercept form
We can change the general form into slope intercept form.
Ax By C 0
A C
m b
Designate B and B
y mx b
y 2 y1 (mx2 b) (mx1 b) m( x2 x1 )
Point- tan m
Slope form x2 x1 x2 x1 ( x2 x1 ) of a straight line
We use this form when we need to find the equation of a line passing through a point ( x1 , y1 )
with slope m : y y1 m( x x1 ) .
Symmetric (Canonical) equation of a line
l
Consider a line on plane XOY .Suppose that we know a point
M 0 ( x0 , y 0 )
and vector (a, b)
that is parallel to the line. This vector is called direction vector of a line. Note, in all likelihood,
will not be on the line itself. We only need to be parallel to the line. Finally, let M ( x, y ) be any
M M ( x x0 ; y y 0 )
point on the line. Since vector 0 and (a, b) are parallel, then their coordinates
are proportional. Therefore
x x0 y y 0 y
a b
This is called the symmetric (canonical) equations of the
line. M
M0
x x1 y y1
x 2 x1 y 2 y1
Normal form
Let length of perpendicular from origin to straight line is ‘P’ and let this perpendicular make
an angle with positive x- axis ‘α’, then equation of a line can be
x y x cos y sin
1 1
p sec p cos ec p p x cos y sin p 0
М0
d
p
Angle between two lines
When two lines intersect, the angle between then is defined as the acute angle between lines
L1 and L2 .
L2
The equation of line L1 is y m1 x b1 L1
φ
The equation of line L2 is y m2 x b2 .
α2
α1
tan 2 - tan 1
tan tan( 2 1 )
1 tan 1 tan 2
Recalling that the tangent of the angle of inclination is the slope of the line, we have
tan 1 m1 , tan 2 m2 , such that m2 m1
m2 m1
tan
1 m1 m2
If and only if 1 2 , the lines are parallel. Then tan 1 tan 2 m1 m2 . The slopes of
parallel lines are the same.
2 1
If and only if 90 , the line are perpendicular. Then
0
2 and
1
tan 2 tan 1 cot 1 m2 m1 1
2 tan 1 .
R(x, y)
If point lies on the bisector, then length of perpendicular from the point R to both the lines
should be equal.
|A 1 x +B1 y +C 1| |A2 x +B2 y +C 2|
=
√ A 2 +B 2
1 1 √
A 2 +B 2
2 2
Generalizing for any point (x, y), the equation of the angle bisector is obtained as
A 1 x+ B1 y+C 1 A2 x + B2 y +C 2
=±
√ A12+ B12 √ A 22 +B 22
This equation gives two bisectors: one-acute angle bisector and the other obtuse bisector.
Example. Find the angle bisector of the angle between the straight lines
x+y-5=0 and 7x−y−19=0
Solution.
x+ y−5 7 x− y−19
=±
√12+12 √7 2+(−1 )2
x+ y−5 7 x− y−19
=±
√2 √50
x+ y−5 7 x − y−19
=
The equation of one bisector is √ 2 5 √2
LECTURE 11
General equation of a plane
P0 ( x0 , y 0 , z 0 )
Suppose that we know a point that is on the plane, . Let’s also suppose that we
have a vector that is orthogonal (perpendicular) to the plane, n ( A, B, C ) .
P ( x, y , z )
This vector is called the normal vector. Now, assume that is any point the plane.
Let’s consider the vector z
M 0 M {x x 0 ; y y 0 ; z z 0 }
.
n
Now, because n is perpendicular to the
plane, it also perpendicular to any vector that lies
in the p plane.
M0 M
The dot product of two perpendicular
M 0M n 0
vectors is equal to zero. So
or y
A( x x0 ) B ( y y 0 ) C ( z z 0 ) 0
(1)
This is called the scalar equation of plane. x
D Ax0 By0 Cz 0
normal form, all the terms of the equation must be multiplied by the normalizing factor
1
A 2 B 2 C 2 where the sign before the radical is opposite to the sign of the constant
term D in general equation of a plane.
Special cases of the position of the plane specified by the general equation
Ax By Cz D 0 :
If in the general equation of a plane the coefficient D 0 ,then, the division of all the term by D
can reduce the equation to the form
A B C
x y z 1 0
D D D
D D D
a b c
Let’s designate A , B , C
We get the intercept form for equation of a plane
x y z
1
a b c
a,b.c are ,respectively , the abscissa, the ordinate and the z- coordinate of the point of intersection
of the plane with the
Ox, Oy, Oz axes.
A1 x B1 y C1 z D1 0
A2 x B2 y C 2 z D2 0
Their normal vectors are
n1 A1i B1 j C1 k and n2 A2 i B2 j C 2 k
1) Two planes are parallel if and only if their normal vectors are parallel i.e. the corresponding
coordinates were proportional
A1 B1 C1
A2 B2 C 2 .
2) Two planes are perpendicular if and only if their normal vectors are perpendicular i.e.
n1 n2 0 or A1 A2 B1 B2 C1C 2 0
3) The angle between two planes is
A1 A2 B1 B2 C1C 2
cos
A B12 C12 A22 B22 C 22
1
2
( x1 , y1 , z1 ), ( x 2 , y 2 , z 2 )( x3 , y 3 , z 3 )
I. Plane through 3 points :
x x1 y y1 z z1
x2 x1 y 2 y1 z 2 z1 0
x3 x1 y3 y1 z 3 z1
and parallel to the vectors 1 (a1 , b1 , c1 ) and 2 (a 2 , b2 , c 2 )
( x0 , y 0 , z 0 )
II. Plane through
x x0 y y0 z z0
a1 b1 c1 0
a2 b2 c2
x x1 y y1 z z1
x 2 x1 y 2 y1 z 2 z1 0
a b c
Solution. We find the normalizing factor ( whose sign s negative since D=21>0)
1 1
2 2 32 6 2 7
Thus, normal form of equation of the given plane is
2 3 6
x y z 30
7 7 7
Example 2. Set up equation of a plane which passes through the point
M (2,3,5) and is
A( x x0 ) B ( y y 0 ) C ( z z 0 ) 0
Solution. It suffices to use equation of a plane which passes
though a given point and is perpendicular to a given vector
4( x 2) 3( y 3) 2( z 5) 0
i.e
4 x 3 y 2 z 27 0
TEST
1) Write down the equation of the plane, that passes through the two points M 1 (1,2,0) and
M 2 (2,1,1) and parallel the direction vectors a (3,0,1) .
A) x 2 y 0 B) x 2 y 3 z 3 0 C) x 5 y 5 z 0
D) x y E) x 5 y 3 z 5 0
M 0 (5,1. 1)
2) Find the distance from the point to the plane x 2 y 2 z 4 0
A) 5 B) 3 C) 2 D) 4 E) 1
A) 14 B) 4,8 C) -1,5 D) 4 E) 6
x x0 y y 0 z z 0
a b c
x x0 y y 0 z z 0 z
t
a b c
P L
We get the parametric representation of L with the
parameter t : P0
x x0 at ; υ
y y 0 bt ;
z z ct y
0
x
If P1 ( x1 , y1 , z1 ) and P2 ( x 2 , y 2 , z 2 ) are two
distinct points in R3, then the equation of line through two points is
x x1 y y1 z z1
x 2 x1 y 2 y1 z 2 z1
For this you need to find an arbitrary point of the line and numbers
( a , b, c )
In this case, the directing vector of the straight line can be found as the vector product of the
normal vectors to the given planes.
i j k
B C1 A1 C1 A1 B1
v n1 n2 A1 B1 C1 i 1 j k i a jb kc
B2 C2 A2 C2 A2 B2
A2 B2 C2
Solution. We find the vector (a, b, c) ,which is parallel to the sought-for straight line. Since it
n 2i j 3k and n2 5i 4 j k of the given
must be perpendicular to the normal vectors 1
planes, we can take as the vector product of the vectors 1 and 2 .
n n
i j k
v n1 n2 2 1 3 11i 17 j 13k
5 4 1
P ( x0 , y 0 , z 0 )
We can take as the point , through which the desired line passes, its point of
intersection with any of the coordinate planes, say the
yOz plane. Since in this case x0 0 , we can
find the coordinates
y0 and z 0 of that point from the system of equations of the given planes if we
x 0
put in them:
y 3 z 1 0
4y z 7 0
x y 2 z 1
11 17 13
1) Two lines are parallel if and only if their direction vectors are parallel and
a1 b1 c1
a 2 b2 c 2
2) Two lines are perpendicular if and only if 1 and 2 are perpendicular and
a1 a 2 b1b2 c1c 2 0
a1 a 2 b1b2 c1c 2
cos
a12 b12 c12 a 22 b22 c 22
x x1 y y1 z z1 x x2 y y 2 z z 2
a1 b1 c1 and a 2 b2 c2
x2 x1 y 2 y1 z 2 z1
a1 b1 c1 0
a2 b2 c2
TEST
2) Write down the equation of the line that passes through the points (2,1,3) and (1,4,3)
x 2 y 1 z 3 x 1 y 4 y 3 x y z 3
A) 1 5 6 B) 5 7 8 C) 5 7 6
x y 4 z 3
D) 1 5 8 E) x y 5
LECTURE 13
Ellipses
Definition of an ellipse. Let P and Q be points in the plane and k a number greater than
the distance from P to Q . The ellipse with foci P and Q is the set of all points X such that the
sum of the distance from X to P and the distance from X to
Q is k written algebraically with X representing a point x, y
on the ellipse
XP XQ k .
The midpoint of the segment PQ joining the foci is the center of the ellipse. The points where the
line through the foci intercept the ellipse are its vertices. The segment connecting the vertices is the
major axis and the segment through the center of the ellipse perpendicular to the major axis is the
minor axis.
k
If the point P and Q coincide, the ellipse generated is a circle with radius 2 . Thus a circle
is a special case of an ellipse.
Equation of an ellipse
The simples case is an ellipse centered at the origin with its foci on the x or y axis.
k
a
Suppose that the foci are on the x axis at the points P c;0 and Q c,0 where c 0 . Let 2
so that k 2a . Then x, y is on the ellipse exactly when XP XQ 2a .
Written algebraically
x c 2 y 02 x c 2 y 02 2a
x c 2 y 2 2a x c 2 y 2 .
Square both sides and simplify the result
x c 2 y 2 4a 2 4a x c 2 y 2 x c 2 y 2
or
x 2 2 xc c 2 y 2 4a 2 4a x c 2 y 2 x 2 2 xc c 2 y 2
a x c 2 y 2 a 2 cx
a 2 x c a 2 y 2 a 4 2a 2 cx c 2 x 2
2
a 2 x 2 2a 2 xc a 2 c 2 a 2 y 2 a 4 2a 2 cx c 2 x 2
a 2
c2 x2 a2 y2 a2 a2 c2 (1)
b 2 x 2 a 2 y 2 a 2b 2 .
2 2
Dividing both sides by a b shows that the coordinates of every point on the ellipse satisfy the
equation
x2 y2
1
a2 b2 standard form of an ellipse.
Characteristics of ellipses
1) x - intercepts: x a
y 2) y -intercepts: y b
dir dir
ect 3) major axis is on the x - axis
ect
b
4) vertices a, 0 and a, 0
rix rix
a -a -c c a
x a
e x
e
-b
x 5) foci c, 0 and c, 0 , where
between foci is 2c
c
e
6) The eccentricity of an ellipse is a , e 1
a a
x x
7) The directrices of an ellipse are the line e, e , that perpendicular x -axis
xx 0 yy 0
2
2 1 x ,y
8) The tangent of ellipse a b in the point 0 0
The standard form of the equation of an ellipse with center( 0 , 0 ) and major axis parallel to the y-axis
x2 y 2
2
+ 2 =1
is b a
where
a>b
the length of the major axis is 2a
the coordinates of the vertices are (0,±a)(0,±a)
the coordinates of the co-vertices are (±b,0)
2 2 2
the coordinates of the foci are (0,±c), where c =a −b
Hyperbolas
Definition of a hyperbola
The hyperbola with foci P and Q is the set of all points X such that the absolute value of
the difference of the distance from X to P and the distance from X to Q is k.
Equation of Hyperbola
The simplest case is a hyperbola centered at the origin with its foci on the x and y - axis.
Suppose that the foci are on the x - axis at the points P c,0 and Qc,0 , where c 0 .
k
a
Let 2 , so that k 2a .
XP XQ 2a
Then
Written algebraically
x c 2 y 2 x c 2 y 2 2a
x c 2 y 2 2a x c 2 y 2
Square both sides
x 2 2 xc c 2 y 2 4a 2 4a x c 2 y 2 x 2 2 xc c 2 y 2
c 2 x 2 2 xca 2 a 4 a 2 x 2 2 xca 2 a 2 c 2 a 2 y
a 2
c2 x2 a2 y2 a2 a2 c2 (2)
b 2 x 2 a 2 y 2 a 2b 2 .
2 2
Dividing both sides by a b shows that the coordinates of every point on the ellipse satisfy
the equation
x2 y2
1
a2 b2 standard from of an hyperbola.
Characteristics of Hyperbola
y b M
y x 1) x - intercepts: x a
a
2) y -intercepts: none
b 3) focal axis is on the axis x - axis
4) vertices a, 0 and a, 0
Parabolas
Definition of a parabola. Let L be a line in the plane and P be a point not on L . If X is any
point not on L , the distance from X to L is defined to be the length of the perpendicular line
segment from X to L . The parabola with focus P and directrix L is the set of all points X such
that distance from X to P = distance from X to L .
The line through P perpendicular to L is called
y the axis.
X
L The intersection of the axis with the parabola
which is midpoint of the segment of the axis
from P to L is the vertex of the parabola.
P x
Equation of a parabola
p
P ,0
Suppose that the focus is on the x - axis at the point 2 , where P is a nonzero
p
x
constant and that the directrix is line 2.
If M x, y is any point on the parabola, then distance from M x, y to the vertical
p
x
line 2 is the length MN . By the definition of a parabola MN MP
p
y MN MQ QN x
directri 2
x N Q
M x , y
p p x
P ,0
2 2
2
p
MP x y 2
2
So,
2
p p
x x y
2
2 2 .
Square both sides of the equation and simplify
2 2
p p
x x y
2
2 2
p2 p2
x xp
2
x xp
2
y2
4 4
y 2 2 xp standard form of parabola
Characteristics of parabola
p
,0
1) focus on the -axis at 2
x
p
x
2) directrix 2
3) axis of symmetry is the x -axis
p
0,
If focus is at 2 on the y axis then
x 2
2 py
4) The tangent of parabola at the point x0 , y 0
yy 0 px x0 .
TEST
2. find the standart form of the ellipse if b 4 and distance between foci 2c 6
x2 y2 x2 y2 x2 y2 x2 y2
1 1 1 1
A) 16 9 B) 16 4 C) 25 16 D) 9 4 E)
x2 y2
1
25 9
x2 y2
1
3. Find the distance between foci 2 c of the ellipse 25 16
A) 9 B) 3 C)8 D) 10 E) 6
x2 y2
1
4. Find the equation of the tangent to the ellipse 18 8
at the point M 0 3,2
x2 y2
1
5. Find the eccentricity of ellipse 25 9
3 3 4 2 9
A) 5 B) 4 C) 5 D) 3 E) 25
x2
y2 1
6. Find the equation of the tangent to the hyperbola 8
at the point M 0 4,1
A)
2x y 1 0 B)
x 2y 1 0 C)
x 2y 1 0
D)
x 2y 2 0 E)
4x 5 y 1 0
x2 y2
1
8. Find the eccentricity of the hyperbola 16 9
3 4 5 4 3
A) 4 B) 5 C) 4 D) 3 E) 2
9. Find the standard form of parabola if the point 3,3 is on the parabola
y 2 px 2
A)
y 2 3x B)
y 2 9x C)
y 2 6x D)
y 2 12 x
E)
y 2 4x
D)
2x y 4 0 E)
2x y 4 0
Lecture 14(LA)
General Equation of a Curve of the Second Order
Let us put down the general form of an equation of an algebraic curve of the second order
Ax 2 2 Bxy Cy 2 Dx Ey F 0 (1)
Now our aim is to transform the Cartesian coordinates in such a way that equation (1) should take
the simplest form; this will enable us to find out what curve is determined by the equation.
The canonical equation containing no terms with the product of the coordinates ,we first of all try to
turn the coordinate axes in such a manner that this product should be eliminated. According to
formulas
x x cos y sin
y x sin y cos
,the equation in new coordinates will have the form
After the axes are turned through an angle
(2)
Removing the parentheses
Ax cos ) 2 2 Axy sin cos A( y sin 2 Bx 2 cos sin 2 Bxy cos 2
2
2 Bxy sin 2 2 By 2 sin cos C x sin 2Cxy sin cos C y cos
2 2
Or
2B
tan 2
2 B cos 2 C Asin 2 0 i.e AC (3)
Now we find the angle from equation (3) and thus determine through what angle the coordinate
axes should be turned.
After the axes are turned through the angle the equation takes the form
Ax 2 C y 2 Dx E y F 0
(4)
Where A, C , D and E are some coefficients which can be found by collecting similar terms in
(2):
But a rotation of the coordinate axes does not change the quantity AC B although the
2
coefficients A, B, C in terms of the second degree may vary , that is the expression is an invariant.
There is no term with x y in equation (4) and therefore, since B 0 , we obtain
AC B 2 AC AC B 2
Consequently, if the expression AC B written for original equation (1) is positive then the
2
coefficients A, C in equation (4) have the same signs since their product is positive. Such a case is
called elliptic.
Finally, if AC B 0 then one of the coefficients AorC is equal to zero( the so-called
2
parabolic case)
Let us turn to the elliptic case. Completing the square in equation (4) we arrive at equation of the
form
Ax a C y b F 0
2 2
x x a
and y y b . Then the equation turns into
Ax 2 C y 2 F
x 2 y 2
1
F F
A C
That is
For the sake of definiteness let us suppose that both A and C are positive. Then if F 0 we
obtain the canonical equation of ellipse. Consequently, the original curve is an ellipse but displaced
and turned with respect to the axes x and y .
3 x 2 4 xy 2 y 2 5 x 4 y 1 0
A 3, B 2, C 2
AC B 2 6 4 0 (ellips)
Example 2.
2 x 2 6 xy y 2 4 x y 3 0
A 2, B 3, C 1
AC B 2 2 9 7 0 (hiperbola)
Example 3.
4 x 2 12 y 9 y 2 5 x 4 y 6 0
A 4, B 6, C 9
AC B 2 36 36 0 (parabola)
4x 2 2 x 9y 2 4 y 4
4x 2 2 x 1 1 9y 2 4 y 4 4 4
4x 1 9 y 2 4 4 36
2 2
4x 1 9 y 2 36
2 2
We carry out a parallel translation of the coordinate axes, taking the point O 1,2 as the new origin
of coordinates and pass to the new coordinates x x 1
and y y 2 . Then the equation
turns into
x 2 y 2
1
4 x 2 9 y 2 36 or 9 4
Thus, we obtain the canonical equation of ellipse.
x 2 2 x 1 1 9y 2 4 y 4 4 44
x 12
9 y 2 44 1 36
2
x 12
9 y 2 9
2
We carry out a parallel translation of the coordinate axes, taking the point O 1,2 as the new
origin of coordinates and pass to the new coordinates x x 1
and y y 2 . Then the
equation turns into
x 2
y2 1
x 2 9 y 2 9 or 9
Thus, we obtain the canonical equation of hyperbola.
x x cos y sin
y x sin y cos
( x cos y sin )( x sin y cos ) 1 0
x 2 cos sin y 2 cos sin xy (cos 2 sin 2 ) 1 0
cos 2 sin 2 0 cos 2 0 2
2 4
1 2 1 2
x y 1 0
2 2
x 2 y 2
1
2 2 .
Lecture 15(LA)
A Quadratic surfaces
In this lecture we will study several families of so-called quadratic surfaces, namely surfaces z =
f(x, y) which are defined by equations of the type
with A, B, C, D, E, F, H, I, J and K being fixed real constants and x, y, z being variables. These
surfaces are said to be quadratic because all possible products of two of the variables x, y, z appear
in (1).
In fact, by suitable translations and rotations of the x, y and z coordinate axes it is possible to
simplify the equation (1) and hence classify all the possible surfaces into the following ten types:
1. Spheres
2. Ellipsoids
3. Hyperboloids of one sheet
4. Hyperboloids of two sheets
5. Cones
6. Elliptic paraboloids
7. Hyperbolic paraboloids
8. Parabolic cylinders
9. Elliptic cylinders
10. Hyperbolic cylinders
You should be able to recognize, classify and sketch at least some of these surfaces The best way to
do that is to look for identifying signs which tell you what kind of surface you are dealing with.
Those signs are:
• The intercepts: the points at which the surface intersects the x, y and z axes
. • The traces: the intersections with the coordinate planes (xy-, yz- and xz-plane).
• The sections: the intersections with general planes.
• The centre: (some have it, some not).
• If they are bounded or not.
• If they are symmetric about any axes or plane
Spheres
x x 0
2
y y 0 z z 0 r 2
2 2
(2)
the point x0 , y 0 , z 0 in the space is the centre of the sphere. The points (x, y, z) in the sphere are all
points whose distance to the centre is given by r. Therefore:
(a) The intercepts of the sphere with the x, y, z-axes are the points x0 r ,0,0 , 0, y 0 r ,0 and
0,0, z 0
r
(b) The traces of the sphere are circles or radius r.
Comparing this equation with the general formula (1) we see it has the form
x 2 y 2 z 2 Hx Iy Jz K 0 , (4)
with H, I, J and K being some real constants. Therefore, the rule to recognize a sphere is the
following: any quadratic equation such that the coefficients of the
x 2 , y 2 and z 2 terms are equal
and that no other quadratic terms exist corresponds to a sphere.
x2 y2 z2
1
a2 b2 c2 (5)
You can see that for a = b = c = 1 we recover the equation of a sphere centered at the origin.
Therefore an ellipsoid is a “deformation” of the sphere such that the sphere gets either stretched or
squeezed (depending on the values of a, b, c) in the x, y, z-directions.
(a) The intercepts of the ellipsoid with the x, y, z-axes are the points
(±a, 0, 0), (0, ±b, 0) and (0, 0, ±c).
(b) The traces of the ellipsoid are ellipses which satisfy the equations:
x2 y2
1
a2 b2 for z = 0 (xy-plane) (6)
x2 z2
1
a2 c2 for y = 0 (xz-plane), (7)
y2 z2
1
b2 c2 for x = 0 (yz-plane). (8)
(c) The sections of the ellipsoid are also ellipses.
(d) The ellipsoid is bounded (all points (x, y, z) in the ellipsoid correspond to finite values of x, y
and z.
(e) The centre of the ellipsoid in the picture is the origin of coordinates. It can be changed by
shifting x, y, z by constant amounts.
(f) The ellipsoid is symmetric about all coordinate planes.
x2 y2 z2
1
a2 b2 c2 (9)
(a) The intercepts of the hyperboloid of one sheet with the x, y, z-axes are the points (±a, 0, 0), (0,
±b, 0). Notice that there is no intersection with the z-axis. The reason is that if we set x = y = 0 in
the equation (9) we obtain the condition z c which admits no real solution for real c.
2 2
(b) The traces of the hyperboloid of one sheet are ellipses in the xy-plane
x2 y2
1
a2 b2 for z = 0 (xy-plane), (10)
and hyperbolas in the xz- and yz-planes:
x2 z2
1
a2 c2 for y = 0 (xz-plane), (11)
y2 z2
1
b2 c2 for x = 0 (yz-plane). (12)
(c) The sections of the hyperboloid of one sheet are ellipses for planes parallel to the xy plane and
hyperbolas for planes parallel to the yz- and xz-planes.
(d) The hyperboloid of one sheet is not bounded.
(e) The centre of the hyperboloid of one sheet in the picture is the origin of coordinates. It can be
changed by shifting x, y, z by constant amounts.
(f) The hyperboloid of one sheet is symmetric about all coordinate planes
x2 y2 z2
1
a2 b2 c2 (13)
(a) The intercepts of the hyperboloid of two sheets with the x, y, z-axes are the points (0, 0, ±c).
There are no intersections with the x, y-axes. The reason is that if we set
z = x = 0 in the equation (9) we obtain the condition
y 2 b 2 which can not be fulfilled for
any real values of y and b. Analogously if we set z = y = 0 we obtain x a which also has
2 2
no real solutions.
(b) In this case we have two sheets, in contrast to all examples we have seen so far. The reason is
that the equation (13) implies
x2 y2 z2
1
a2 b2 c2 (14)
2
z
2
1 0 z c
This equation can only be solved if c which implies that
(c) The traces of the hyperboloid of two sheets are hyperbolas in the xz- and yz-planes:
x2 z2
1
a2 c2 for y = 0 (xz-plane), (15)
2 2
y z
2
2 1
b c for x = 0 (yz-plane). (16)
(d) The sections of the hyperboloid of two sheets are hyperbolas for any planes parallel to the xz-
z c
or yz-planes and ellipses for planes parallel to the xy-plane with .
(e) The hyperboloid of two sheets is also not bounded.
(f) The centre of the hyperboloid of two sheets in the picture is the origin of coordinates. It can be
changed by shifting x, y, z by constant amounts.
(g) The hyperboloid of two sheets is again symmetric about all coordinate planes. General rule:
Any quadratic surface such that: the coefficients of
x2 , y2 2
and z are different, one of the
coefficients is negative and two of the coefficients are positive, and no other quadratic terms
appear describes a hyperboloid. In addition, if the constant term has negative sign the
hyperboloid has two sheets whereas if the constant term has positive sign the hyperboloid has
only one sheet.
Cones
A cone is a quadratic surface whose points fulfil the equation
x2 y2
2
2 z2 0
a b (17)
Comparing (17) with the equations for the hyperboloids of one and two sheet we see that the
cone is some kind of limiting case when instead of having a negative or a positive number on
the l.h.s. of the quadratic equation we have exactly 0.
(a) The only intercept of the cone with the x, y, z-axes is the origin of coordinates (0, 0, 0).
(b) The traces of the cone are lines in the xz- and yz-planes
x
z
a for y = 0 (xz-plane), (18)
y
z
b for x = 0 (yz-plane), (19)
and the origin (0, 0) in the xy-plane:
x2 y2
0
a2 b2 for z = 0 (xy-plane). (20)
(c) The sections of the cone are lines for any planes parallel to the xz- or yz-planes and ellipses
(or circles if a = b) for planes parallel to the xy-plane.
(d) The cone is not bounded.
(e) The centre of the cone in the picture is the origin of coordinates. It can be changed by
shifting x, y, z by constant amounts.
(f) The cone is symmetric about all coordinate planes.
Hyperbolic paraboloids
A hyperbolic paraboloid is defined by the equation
x2 y2
z
a2 b2 (25)
(a) The only intercept of the hyperbolic paraboloid with the x, y, z-axes is the origin of
coordinates (0, 0, 0).
(b) The traces of the paraboloid are parabolas in the xz- and yz-planes
x2
z 2
a for y = 0 (xz-plane), (26)
2
y
z 2
b for x = 0 (yz-plane), (27)
and two lines in the xy-plane corresponding to the equation:
x2 y2
0
a2 b2 for z = 0 (xy-plane), (28)
(c) The sections of the hyperbolic paraboloid with planes parallel to the xz- or yz-planes are
parabolas. The sections with planes parallel to the xy-plane are hyperbolas.
(d) The paraboloid is not bounded.
(e) The centre of the paraboloid in the picture is the origin of coordinates which in this case is
called a saddle point. It can be changed by shifting x, y, z by constant amounts.
(f) The hyperbolic paraboloid is symmetric about the xz- and yz-planes.
Rule: A hyperbolic paraboloid has the same features as the elliptic paraboloid with the only
difference that the coefficients of the quadratic terms (
x 2
, y 2
in our example) have opposite
signs.
Parabolic cylinders
Elliptic cylinders