add another method: matrix computation