Suppose we have a function consisting of a series of matrices multiplied by a vector: f(X) = A * B * b --where X is a vector containing elements that are contained within A, b, and/or b, --A is a matrix, B is a matrix, and b is a vector Each Matrix and the vector is expressed as more terms, ie... X = (x1, x2, x3) A = [ x1 + y1 y4 y7 ] [ y2 x2 + y5 y8 ] ] y3 y6 x3 + y9 ] B = [ y1 x2 + y4 x3 + y7 ] [x1 + y2 y5 y8 ] ] y3 y6 y9 ] b = [y1 y2 y3]' (' means transposed) Now we want to find the Jacobian of f - ie the p