From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 26/CH2/EX2.4.25/2_4_25.sce | 29 +++++++++++++++++++++++++++++ 1 file changed, 29 insertions(+) create mode 100755 26/CH2/EX2.4.25/2_4_25.sce (limited to '26/CH2/EX2.4.25') diff --git a/26/CH2/EX2.4.25/2_4_25.sce b/26/CH2/EX2.4.25/2_4_25.sce new file mode 100755 index 000000000..985a344eb --- /dev/null +++ b/26/CH2/EX2.4.25/2_4_25.sce @@ -0,0 +1,29 @@ +disp('given matrix is:') +a=[1 2 0 0 0;3 5 0 0 0;0 0 2 0 0;0 0 0 7 8;0 0 0 5 6]; +disp(a) +disp('partitioning the matrix into 4 submatrices') +A11=[a(1,1:2);a(2,1:2)] +disp(A11,'A11=') +A22=[a(3,3:5);a(4,3:5);a(5,3:5)] +disp(A22,'A22=') +A12=zeros(2,3) +disp(A12,'A12=') +A21=zeros(3,2) +disp(A21,'A21=') +disp('partitioning A22 into 4 submatrices') +A221=[2] +disp(A221) +B=[A22(2,2:3);A22(3,2:3)] +disp(B,'B=') +disp(zeros(1,2)) +disp(zeros(2,1)) +disp('determinant of B=') +disp(det(B)) +disp('Hence, B is invertible') +disp('inverse of B is') +disp(inv(B)) +disp('determinant of inverse of B is:') +disp(det(inv(B))) +disp('hence the invese of A22 is:') +c=[det(inv(B)) zeros(1,2);0 3 -4;0 -2.5 3.5]; +disp(c) \ No newline at end of file -- cgit