Suppose I have a square matrix P

P <- matrix(c(.3,.7, .7, .3), ncol=2)

I know that 

> P * P 

Returns the element by element product, whereas

> P%*%P

Returns the matrix product.

Now, P^2 also returns the element by element product. But, is there a
slick way to write

P %*% P %*% P

Obviously, P^3 does not return the result I expect.

Thanks,
Harold






	[[alternative HTML version deleted]]

