Create a function csolve2.m with input matrix A and a column vector b and no output. The function should display the rank, the pivot and free variables of A, the particular solution to Ax = b, and the special solutions for A. Test it using the system below. A = [1 0 1 -8 0 -8 -7 -3 -10 5 2 7], b = [4 -3 1], Your display should be precisely the following: >> csolve2(A, b) The rank of the coefficient matrix is 2. Pivot variables 1 3 Free variables: 2 4 The particular solution is: xp = 11.0000 0 1.0000 0 The special solutions are: XB = 8.0000 -0.3333 1.0000 0 0 0.6667 0 1.0000