Dear all R users, I was trying to create a polymonial using the polynom() of "PolynomF" package. My problem is, if I pass coefficients as simple numerical values, it is working, but for matrix coefficient it is not. Here is my try: library(PolynomF) z <- polynom() ## following is working p1 <- 1- 4.5*z ## However following is not giving results as intended p1 <- diag(3)- matrix(1:9, 3,3)*z In the 2nd example, I was expecting a matrix polynomial in "z" upto degree 1, which I did not get. Can anyone please correct me how to address this problem? Thanks, -- View this message in context: http://r.789695.n4.nabble.com/Problem-to-create-a-matrix-polynomial-tp2993411p2993411.html Sent from the R help mailing list archive at Nabble.com.