Tutorial 1,  pg  .40

Question 1;
a = {{3, 4, p^2 + 1}, {2, 3, 1}, {1, 
            4, -1}};
b = {{1, -5, 1}, {2, 1, -3}, {4, -4, 
            3}};
c = {{5, 1, 2}, {4, 0, 3}, {2, 
            4, -1}};
a1 = {{3, q^2, 10}, {2, 3, 1}, {1, 
            4, -1}};
A = MatrixForm[a];
B = MatrixForm[b];
CC = MatrixForm[c];
A1 = MatrixForm[a1];
Print["A=", A];
Print["B=", B];
Print["CC=", CC];
Print["A[1,3]=", a[[1, 3]]];
Print["A1[1,3]=", a1[[1, 3]]];

    
1 (a);
psol = 
        Solve[a[[1, 3]] == a1[[1, 3]], 
          p];
Print["psol=", psol];
qsol = 
        Solve[a[[1, 2]] == a1[[1, 2]], 
          q];
Print["qsol=", qsol];
Print["Det[A]=", Det[a]];

1 (b);
D1 = a c  +  b;
Print["D1=", MatrixForm[D1]];

1 (c) ;
D2 =  c a  +  b;
Print["D2=", MatrixForm[D2]];

1 (d);
detA = Det[a]; detB = Det[b]; detC = Det[c];
Print["Det[A]=", detA];
Print["Det[B]=", detB];
Print["Det[C]=", detC];