Tải bản đầy đủ (.doc) (13 trang)

HƯỚNG DẪN THỰC HÀNH CƠ BẢN MATHEMATICA ppt

Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (141.21 KB, 13 trang )

HƯỚNG DẪN THỰC HÀNH CƠ BẢN MATHEMATICA
I/ Các phép toán số học
+, -, *, /, ^
II/Cách khai báo các hàm số cơ bản (có sẵn)
Hàm số cơ bản Được khai báo bằng lệnh
( )f x x=
f[x_]:=Abs[x]
( )f x x=
f[x_]:=Sqrt[x] hoặc f[x_]:=x^(1/2)
( )
n
m
f x x=
f[x_]:=x^(m/n)
( )f x sinx=
f[x_]:=Sin[x]
( )f x cosx=
f[x_]:=Cos[x]
( )f x tgx=
f[x_]:=Tan[x]
( )f x cotgx=
f[x_]:=Cot[x]
1
osx
( ) sec
c
f x x= =
f[x_]:=Sec[x]
1
sin
( ) csc


x
f x x
= =
f[x_]:=Csc[x]
( )f x arcsinx=
f[x_]:=ArcSin[x]
( )f x arccosx=
f[x_]:=ArcCos[x]
( )f x arctgx=
f[x_]:=ArcTan[x]
( )f x arccotgx=
f[x_]:=ArcCot[x]
( )
a
f x log x=
f[x_]:=Log[a,x]
( )f x lgx=
f[x_]:=Log[10,x]
( )f x lnx=
f[x_]:=Log[E,x] hoặc f[x_]:=Log[x]
( )
x
f x a=
f[x_]:=a^x
( )
x
f x e=
f[x_]:=E^x
2
( ) ( )

x x
e e
f x ch x

+
= =
f[x_]:=Cosh[x]
2
( ) ( )
x x
e e
f x sh x


= =
f[x_]:=Sinh[x]
( )
( )
( ) ( )
x x
x x
sh x
e e
ch x
e e
f x th x



+

= = =
f[x_]:=Tanh[x]
( )
( )
( ) ( )
x x
x x
ch x
e e
sh x
e e
f x cth x


+

= = =
f[x_]:=Coth[x]
1
III/Cánh khai báo một hàm số mới
1/ Khai báo hàm giá trị thực, biến thực
Thí dụ1(hàm một biến): hàm
3
( ) . . .
x
f x x sinx ln x e cosx= +
được khai báo bằng lệnh
f[x_]:=x*Sin[x]+(Log[x]^3)*(E^x)*Cos[x]
Thí dụ 2 (hàm nhiều biến):
Các hàm

2 2
( , ) . .f x y x y y sin x= +
,
( , , ) . . .
xy
g x y z x cotgy y cotgz z e= + +
được khai báo như
sau.
f[x_,y_]:= x*y^2+y*(Sin[x])^2 ]
g[x_,y_,z_]:=x*Cot[y]+y*Cot[z]+z*E^(x*y)
2/Khai báo hàm thực biến véc tơ
Thí dụ 3: Khai báo hàm chuẩn một biến véc
( )f x x

=
,
1 2
( , , , ) R
T n
n
x x x x= ∈
như
sau “giả thiết đã nhập n trước đó”(theo lý thuyết
{1, ,n}i
i
x Max x


=
)

f[x_]:=Max[ Table[ Abs[ x[[i]] ] , {i,1,n}] ]
Thí dụ 4: Khai báo hàm chuẩn hai biến véc tơ tơ (áp dụng tính khoảng cách giữa
hai điểm
x

y
trong không gian định chuẩn
R
n
)
( , )f x y x y

= −
,
1 2
( , , , ) R
T n
n
x x x x= ∈
,
1 2
( , , , ) R
T n
n
y y y y= ∈
“giả thiết đã nhập n trước đó” (theo lý
thuyết
{1, n}i
i i
x y Max x y



− = −
L
)
f[x_,y_]:=Max[ Table[ Abs[ x[[i]]-y[[i]] ] ,{i,1,n}] ]
Chú ý Giá trị của n có thể lấy bằng lệnh tính số phần tử của ma trận cột x
Length[x]
Thí dụ 5: Khai báo hàm chuẩn biến ma trận.
Cho ma trận
( )
ij
m n
A a
×
=
, theo lý thuyết
{1, ,m}
1
( )
i
n
ij
j
A Max a


=
=


. Khi đó hàm chuẩn
của ma trận này được được khai báo như sau (giả thiết đã nhập trước đó các giá
trị của m, n)
2
f[A_]:=Max[ Table[ Sum[Abs[A[[i,j]]],{j,1,n}] , {i,1,m}] ]
3/ Khai báo hàm giá trị véc tơ
Thí dụ 6: Khai báo hàm
1
.
2
3
.
F( ,y,z)= .
. .
y z
x y x z
f
x f x e
x siny y cosz
f
+ +
 
 
 ÷
 ÷
=
 ÷  ÷
+
 
 

bằng lệnh
F[x_,y_,z_]:={ x+y+z , x*E^(y*z) , x*Sin[y]+y*Cos[z] }
IV/Giải toán bằng Mathematica
1/ Giải toán đại số và giải tích
1.1/ Vẽ đồ thị hàm số trong mặt phẳng( trong không gian hai chiều)
-Vẽ đồ thị hàm một biến
( ), [a,b]y f x x= ∈

Plot[ f[x] , {x,a,b} ]
-Vẽ trên cùng một hệ trục tọa độ đồ thị của hai hàm số
( ), ( )y f x y g x= =
,
[a,b]x ∈
.
Plot[ {f[x],g[x]} , {x,a,b} ]
-Vẽ đồ thị của hàm cho bởi phương trình tham số ( hoành độ và tung độ là hàm
biến t )
{
( )
, [a,b]
( )
x x t
t
y y y
=

=
ParametricPlot[ {x[t],y[t]} , {t,a,b} ]
1.2/ Vẽ đồ thị hàm số trong không gian ba chiều
- Vẽ đồ thị của một hàm hai biến(đồ thị là một mặt trong không gian ba chiều)

( , )z f x y=
,
[a,b], [c,d]x y∈ ∈
.
Plot3D[ f[x,y] , {x,a,b} , {y,c,d} ]
- Vẽ đồ thị của một mặt cho bởi phương trình tham số(cả ba tọa độ của điểm
thuộc mặt đó đều là hàm của hai biến t,s. (t và s là hai tham số).
( , )
( , )
( , )
x x t s
y y t s
z z t s
=


=

=


,
[a,b]
s [c,d]
t ∈

3
ParametricPlot3D[ {x[t,s],y[t,s],z[t,s]} , {t,a,b} , {s,c,d}]
1.3/ Các giới hạn
( )

x a
lim f x
−>
,
( )
x a
lim f x

−>
,
( )
x a
lim f x
+
−>
,
( )
x
lim f x
−>+∞
,
( )
x
lim f x
−>−∞
lần lượt được tính bằng các lệnh.
Limit[f[x],x->a]
Limit[f[x],x->a, Direction->-1]
Limit[f[x],x->a, Direction->1]
Limit[f[x],x->Infinity]

Limit[f[x],x-> -Infinity]
1.4/ Tính đạo hàm cấp n của hàm
f
(
f
có thể có nhiều biến ) theo biến
x
bằng
lệnh
D[ f , {x,n} ]
Chú ý : Nếu tính đạo hàm cấp 1 có thể dùng lệnh D[ f ,x]
1.5/ Tính đạo hàm của hàm véc tơ
Thí dụ 7: Cho hàm số
1
.
2
3
.
F( ,y,z)= .
. .
y z
x y x z
f
x f x e
x siny y cosz
f
+ +
 
 
 ÷

 ÷
=
 ÷  ÷
+
 
 
Sau khi khai báo hàm bằng lệnh
F[x_,y_,z_]:={x+y+z,x*E^(y*z),x*Sin[y]+y*Cos[z]}
Ta tính được các ma trận cột
1
2
3
f
x
f
x
f
x






 
 ÷
 ÷
 ÷
 ÷
 

,
1
2
3
f
y
f
y
f
y






 
 ÷
 ÷
 ÷
 ÷
 

1
2
3
f
z
f
z

f
z






 
 ÷
 ÷
 ÷
 ÷
 
lần lượt bằng lệnh
D[F[x,y,z],x]
D[F[x,y,z],y]
D[F[x,y,z],z]
Từ đó muốn tính ma trận đạo hàm
4
1 1 1
2 2 2
3 3 3
F ( , , )
'
f f f
x y z
f f f
x y z
f f f

x y z
x y z
∂ ∂ ∂
∂ ∂ ∂
∂ ∂ ∂
∂ ∂ ∂
∂ ∂ ∂
∂ ∂ ∂
 
 ÷
 ÷
=
 ÷
 ÷
 ÷
 ÷
 
thì dùng lệnh
Transpose[ { D[F[x,y,z],x] , D[F[x,y,z],y] , D[F[x,y,z],z] } ]
Tuy nhiên, nếu muốn tính
F (1,3,7)
'
ta phải có mẹo nhỏ như sau
dh[x_,y_,z_]:=Transpose[{D[F[u,v,w],u] ,D[F[u,v,w],v] ,D[F[u,v,w],w]}]/.{u->x,v->y,w->z}
dh[1,3,7]
Chú ý: Người đọc tự kiểm tra nếu ta tính
F (1,3,7)
'
bằng một trong hai cách sau
thì mắc phải sai lầm gì ?

Cách 1: Transpose[{D[F[1,3,7],x] ,D[F[1,3,7],y] ,D[F[1,3,7],z]}]
Cách 2: dh[x_,y_,z_]:=Transpose[ {D[F[x,y,z],x] , D[F[x,y,z],y] , D[F[x,y,z],z]} ]
dh[1,3,7]
1.5/ Tính nguyên hàm của hàm
( )f x
theo biến
x
bằng lệnh
Integrate[ f[x] , x]
1.6/ Tính tích phân của hàm
( )f x
, trên đoạn
[a,b]
(kết quả là số thập phân)
bằng lệnh
NIntegrate[ f[x] , {x,a,b} ]
2/ Giải toán đại số tuyến tính
2.1/ Khai báo các ma trận biết trước các phần tử
5
Thí dụ 8: các ma trận
[ ]
X= 2 1 4 3
,
3
4
5
Y=
12
10
21

 
 
 
 
 
 
 
,
1 2 4
5 2 4
2 1 7
A
 
 
=
 
 
 
được khai báo lần lượt
bằng lệnh
X={2,1,4,3};
Y={{3},{4},{5},{12},{10},{21}};
A={{1,2,4},{5,2,4}, {2,1,7}};
Muốn lấy phần tử thứ k của X ta dùng lệnh
X[[k]]
Muốn lấy phần tử thứ k của Y ta dùng lệnh
Y[[k,1]]
Muốn lấy phần tử hàng i cột j của ma trận A ta dùng lệnh
A[[i,j]]
Chú ý : Trong khai báo trên X được hiểu là ma trận 1 cột 4 hàng, Y được hiểu là

ma trận 1 cột 6 hàng, A là ma trận 3 hàng 3 cột.
2.2/ Khai báo ma trận chỉ biết trước cỡ của ma trận, còn giá trị của phần tử trên
mỗi hàng, mỗi cột chưa biết.
- Sau khi khai báo giá trị của m và n thì khai báo ma trận A có m hàng n cột
bằng lệnh.
A=Table[a[[i,j]],{i,1,m},{j,1,n}]
2.3/ Khai báo ma trận đặc biệt
- ma trận đơn vị cấp n
IdentityMatrix[n]
6
- ma trận vuông cấp 5 mà các phần tử nằm trên đường chéo lần lượt là
a,b,c,d,e. Các phần tử nằm ngoài đường chéo bằng không.
DiagonalMatrix[a,b,c,d,e]
2.4/ Các phép toán ma trận
- Phép công, trừ và nhân hai ma trận A với B được thực hiện bởi lệnh
A+B
A-B
A.B
Chú ý : Nhân hai số thực là dấu sao còn nhân hai ma trận là dấu chấm.
- Chuyển vị của ma trận A
Transpose[A]
- Nếu A là ma trận khả nghịch thì ma trận nghịch đảo được tính
Inverse[A]
- Nếu A là ma trận vuông thì định thức được tính
Det[A]
Chú ý:
- Có bốn loại móc được sử dụng trong chương trình Mathematica
{ } cặp móc nhọn
[ ] cặp móc vuông đơn
[[ ]] cặp móc vuông kép

( ) cặp móc tròn
- Dấu móc vuông đơn [ ] bắt buộc phải dùng khi gọi một hàm có sẵn hoặc
khai báo một hàm số mới như trong các thí dụ trên.
- Dấu móc vuông kép [[ ]] bắt buộc được sử dụng khi gọi đến một phần tử
của một ma trận.
- Dấu móc nhọn { } được sử dụng theo đúng mẫu lệnh của phần mềm
7
Thí dụ 9: Lệnh tính tích phân NIntegrate[ f[x] , {x,a,b} ]
Lệnh nhập ma trận ở trên X={2,1,4,3}
- Dấu móc tròn ( ) được sử dụng khi trình bày các biểu thức toán
Thí dụ 10: Biểu thức
1 3
2
.( lg )
x
e sin x
lnx tgx x
+
+
+
được viết trong Mathematica như sau
((E^(x+1))+Sin[x]^3)/(Log[x](Tan[x]+Log[x]^2))

Sau khi bấm tổ hợp phím Shift +Enter ta có kết quả
Lỗi thường gặp khi gõ công thức trên
Gõ sai Gõ đúng
e^(x+1) E^(x+1)
E^x+1 E^(x+1)
E^[x+1] E^(x+1) hoặc Exp[x+1]
Sin^3[x] (Sin[x])^3 hoặc Sin[x]^3

Sin(x)^3 Sin[x]^3
sin[x]^3 Sin[x]^3
Ln(x) Log[x] hoặc Log[E,x]
ln[x] Log[x]
Log^2[x] Log[x]^2
lg[x] Log[10,x]
2.5/ Lệnh giải hệ phương trình A.X=B sau khi đã nhập hai ma trận A và B
LinearSolve[A,B]
8
V/ Lập trình đơn giản hỗ trợ môn Phương pháp tính
Vì hai lý do sau
- Thời gian của môn Phương pháp tính chỉ có hạn
- Phần mềm Mathematica chỉ là một trong những công cụ hỗ trợ cho môn học
cho nên ở đây chúng ta chỉ xét một số lệnh cơ bản đủ để giải quyết các bài toán
đặt ra của môn học. Những ai có nhu cầu cao hơn có thể tự học tập thêm qua
phần hướng dẫn Help của phần mềm Mathematica

1/ Muốn lặp lại các công việc “việc 1, việc 2, …, việc k” n lần ta dùng lệnh Do
như sau
Do[việc 1;việc 2;…; việc k, {n}]
Chú ý – các công việc cách nhau bởi dấu chấm phẩy
– Công việc cuối cùng được kết thúc bởi dấu phẩy
Các dạng khác của lệnh Do
- Thực hiện công thức theo một chỉ số A[i] trong đó
1,i n=
bằng lệnh
Do[A[i] ,{i,1,n}]
- Thực hiện công thức theo hai chỉ số A[i,j] trong đó
1,i n=
,

1,j m=
bằng lệnh
Do[A[i,j] ,{i,1,n},{j,1,m}]
2/ Chừng nào biểu thức lôgic “bt” còn có giá trị đúng thì ta còn thực hiện lặp lại
các công việc “việc 1, việc 2, …, việc k” . Khi đó ta sẽ dùng lệnh While để lập
trình như sau
While[bt ,việc 1;việc 2;…; việc k]
Chú ý:
- Sau bt là dấu phẩy.
- Giữa hai việc cách nhau bởi dấu chấm phẩy.
3/ Trong mỗi bươc lặp đôi khi ta cần tăng giá trị của biến nguyên n thêm p đơn
vị ta dùng lệnh sau
9
n+=p
4/ - Để in ra màn hình gia trị của một biến x ta dùng lệnh Print như sau
Print[x]
- Để in ra màn hình một đoạn ký tự ksahhoiuhlszscl ta dùng lệnh Print như
sau( trong lệnh Print đoạn ký tự được đặt trong dấu nháy kép)
Print["ksahhoiuhlszscl "]
Chú ý: -nếu in cả ký tự và biến thì giữa ký tự và biến được cách nhau bởi dấu
phẩy
Thí dụ lệnh sau sẽ in ra cả ký tự ksahhoiuhlszscl và biến x
Print["ksahhoiuhlszscl " , x]
- Kết thúc lệnh Print không có dấu chấm phẩy
VI. LUYỆN TẬP TRÊN MATHEMATICA 2.2
A/ Vẽ đồ thị
-Vẽ đồ thị các hàm trong bài tập 1, 2, 5 trang 90, 91 (giáo trình Doãn Tam Hòe)
B/ Tính giới hạn
- Làm bài tập 1 trang 112 (giáo trình Doãn Tam Hòe )
C/ Tính đạo hàm

- Làm bài tập 2 trang 112, 113 (giáo trình Doãn Tam Hòe )
- Tính các ma trận
'
( )F X

' (0) 1
( ( ))F X

, trong đó
1.
2 2 2
1
2 2
2
( , )
19( )
( )
( , )
7( )
f x y
x xy y x y
F X
f x y
x xy y x y
 
+ + − −
 
= =
 ÷
 ÷

− + − −
 
 
,
(0)
(19.0,40.0)
T
X =

10
2.
5 4 2
1
5 4 2
2
5 4 2
3
2 2
( , , )
( ) ( , , ) 2 2
( , , )
2 2
x y x y
f x y z
F X f x y z y z y z
f x y z
z x z x
 
− + −
 

 ÷
 ÷
= = − + −
 ÷
 ÷
 ÷
− + −
 
 
,
(0)
(3,1, 2)
T
X = −
D/ Tính tích phân
1.
3 2 2 2
1
( )x x a
dx
+

2.
3
2 2 2
( )
x
x a
dx



3.
4
2 2 2
( )
x
x a
dx


4.
3
3 3
x
x a
dx
+

5.
3
3 3 2
( )
x
x a
dx
+


6.
3 3 3 2

1
( )x x a
dx
+

7.
4 4 4 2
1
( )x x a
dx


8.
3
.
x
a x b
dx
+

9.
3 2
1
( . )x a x b
dx
+


10.
. .x a x b dx+


11.
3
. .x a x bdx+

12.
3
.
x
a x b
dx
+

13.
5
(2 3)
7. 2
x
x
dx
+
+

14.
5
1
(2 1) 7. 2x x
dx
+ +


15.
7
3 4
(5 8)
x
x
dx
+
+

16.
2 11
7
(4 )x
x
dx


17.
3
3 2
1
. (5 )x x
dx


18.
12
2 13
(2. 3 9)

x
x x
dx
+ −


19.
3 7
.sin (3 )x x dx

20.
9
5
sin (3 )x
x
dx

21.
5
17
sin (2 )
x
x
dx

22.
2
1
(2sin 3 4cos2 ) 7x x
dx

+ +

23.
2 3
1
(5sin 3 4)x
dx
+

24.
7
2
tan 3
cos 2
x
x
dx

25.
7
21
tan 3
sin 3
co x
x
dx

26.
7
2

3
tan 3
cos 2
x
x
x dx

27.
2
2
3
x
e
x
dx


28.
3
.ln
x
e xdx


29.
3 5
.sin 2
x
e xdx



30.
3
.ln( 2)x x dx

+

31.
3 5 2 2
.ln ( )x x x a dx

+ +

32.
6 5
.ln (3 4)x x dx

+

33.
3
2 2
2 2
.ln( )
x
x a
x x a dx


+ +


34.
5
x
3
.arcsin( )x dx


35.
5
3
x
.arcsin( )x dx


36.
5 2 7
x
3
.( 9) .arcsin( )x x dx+

37.
5 2 4
x
3
.( 9) .arcsin( )x x dx

+

38.

5 2 4
x
3
.( 9) .arccotan( )x x dx

+

39.
8 2 4
x
3
.( 9) .arctan( )x x dx

+

40.
20
1
x
x
e
dx
+

41.
/ 2
2
1
4
0

1 sin xdx
π


42.
1
arctan
0
x
x
dx

43.
2
1
0
x
e dx

44.
4
0
2x dx−

45.
1 cos 2
2
0
x
dx

π
+

46.
2
1
ln
e
x x
e
dx

47.
3
1
1 5 1
0
x x
dx
+ + +

48.
1
arcsin
1
0
x
x
dx
+


49.
2
/ 3
sin
cos
/ 3
x x
x
dx
π
π


50
2
1
ln(1 )
1
0
x
x
dx
+
+

51.
ln 2
0
1

x
e dx−

52.
2
1
cos
1
1
x x
x
dx
+


53.
2
/ 2
1 2
/ 2
sin
x
x
x dx
π
π
+


54.

3
/ 2
cos
cos sin
0
x
x x
dx
π
+

55.
2
3
2
1
sin
1
0
x
x
dx


11
56.
1
2
0
lnx xdx


57.
3
1
. ln
1
e
x x
dx

58.
1
cos(ln )
e
x dx
π

59.
2
sin
1 cos
0
x x
x
dx
π
+

60.
3

1
1
0
x
x
e
e
dx


+

61.
1
2006
0
1x xdx−


62.
2
0
sin(2002 sin )x x dx
π
+

63.
0
( sin 2006 sin 2007 )x x dx
π



64.
1
2006
0
sin( )x x dx
π


E/ Nhập vào các ma trận và thực hiện các tính toán cộng, trừ, nhân ma trận.
Tính ma trận chuyển vị, định thức và ma trận ngược của một ma trận nếu có
Bài 1: Tính A+B, A-B, 3A , 3A-3B,
T
A
,
T
B

(3A-3B)
T
1.
1 1
2 1
A

 
=
 


 
,
2 1
1 8
B

 
=
 

 
2.
1 2
2 1
A
 
=
 
 
,
3 2
4 2
B
− −
 
=
 
 
3.
6 1

2 4
3 5
A

 
 
=
 
 

 
,
1 4
1 5
1 10
B
 
 
= −
 
 
 
4.
2 1 1
1 1 4
A
 
=
 
− −

 
,
2 3 4
3 1 2
B

 
=
 
− −
 
5.
2 2 1 0 1
1 1 2 0 1
A

 
=
 
− −
 
,
1 1 1 1 0
3 4 9 6 7
B

 
=
 
− − −

 
Bài 2: Cho
2 1
1 0
3 4
A
− −
 
 
=
 
 

 

0 3
2 0
4 1
B
 
 
=
 
 
− −
 
, hãy tính
1.
3 2A B−
2.

2 4A B
+
3.
3
1
2 2
A B

+

Bài 3: Tính ma trận
ij
( ) .
m n
C c A B
×
= =
.
Dùng lệnh Sum[] để tính
12
c
,
32
c

23
c
1.
0 1 0
4 0 2

8 1 7
A

 
 
=
 
 

 
,
2 1
3 4
1 6
B
 
 
= −
 
 
 
2.
5 0 0
0 8 0
0 0 7
A
 
 
= −
 

 
 
,
1
5
1
2
0 0
0 8 0
0 0
B
 
 
= −
 
 
 
 
Bài 4: Tính định thức của các ma trận sau
1.
1
5
1
2
0 0
0 8 0
0 0
 
 


 
 
 
 
2.
6 3 7
0 0 0
4 6 3

 
 
 
 

 
3.
0.3 0.2 0.2
0.2 0.2 0.2
0.4 0.4 0.3
 
 
 
 

 
4.
6 0 3 5
4 13 6 8
1 0 7 4
8 6 0 2


 
 

 
 

 
 
12
Bài 5: Tìm ma trận nghịch đảo nếu có của các ma trận sau
1.
2 17 11
1 11 7
0 3 2

 
 
− −
 
 

 
2.
2 0 1 1
3 0 0 1
1 1 2 1
4 1 1 0
 
 

 
 
− −
 

 
3.
1 2 1 1
4 9 5 6
0 1 1 1
3 5 3 3
− − −
 
 
− − −
 
 
− −
 

 
4.
1
3
3 1 1 3
3 1 2 3
.
0 1 1 0
3 2 1 0
− −

 
 
− − −
 
 
 
− −
 
Bài 6: Tính chuẩn của các ma trận trong bài tập 5 theo chuẩn vô cùng

Bài 7: Giải hệ phương trình
.A X B=
bằng lệnh
LinearSolve[ , ]A B
1.
2 17 11
1 11 7
0 3 2
A

 
 
= − −
 
 

 
,
1
2

3
B
 
 
= −
 
 
 
2.
2 0 1 1
3 0 0 1
1 1 2 1
4 1 1 0
 
 
 
 
− −
 

 
,
1.4
2.6
1
4
B
 
 
 

=
 

 

 

3.
1 2 1 1
4 9 5 6
0 1 1 1
3 5 3 3
A
− − −
 
 
− − −
 
=
 
− −
 

 
,
3
2
2
7
B

 
 
 
=
 

 

 
4.
1
3
3 1 1 3
3 1 2 3
.
0 1 1 0
3 2 1 0
− −
 
 
− − −
 
 
 
− −
 
,
2
4
2

5
B

 
 
 
=
 

 

 
13

×