对角矩阵，全1矩阵，（ones,diag）

这段代码看不懂，有哪位可以帮忙解释一下嘛，p0+2这个地方也不是很明白，感谢，鞠躬鞠躬！
（1）能明白ones那个地方是要生成一个3×3的矩阵，下面对角矩阵diag那个地方不是很懂。
（2）这个程是三个方程，我要改成5个方程的程序，是不是应该ones那里生成5×5的矩阵呢？p0+2那里我又应该怎么处理呢？