Steps To Reproduce

If we run

Z=matrix([[0,1],[0,0]])
n=var('n')
Zn=Z^n;Zn

we get the null matrix.

Expected Behavior

I would have expected a result like
matrix([0^n,n*0^(n-1)],[0,0^n]])
to see that the matrix is not automatically zero for small n.

Actual Behavior

The actual result is
matrix([[0,0],[0,0]])

Additional Information

I highly suspect that it is caused by the .simplify_full() method which turns 0^n into 0.

Environment

- **OS**: Ubuntu 20.04
- **Sage Version**: 9.0
- Jupyter Notebook

Checklist

• I have searched the existing issues for a bug report that matches the one I want to file, without success.
• I have read the documentation and troubleshoot guide