求助：函数的循环

[EPFL_GO111]

如题，我想实现这样的循环, f=r*x^2+y^2+z, for r=1:1:5;
跪求意见，谢谢大家~

[huleigxut]

有哪位大侠知道，使用mdiag函数生成一个可变化的方块对角阵！
例如：
f1=[1 2;
3 4];
f2=[5 6;
7 8]
利用mdiag后，即mdiag(f1,f2),结果：[ 1 2 0 0
3 4 0 0
0 0 5 6
0 0 7 8]
现在我要求mdiag括号里面的方块阵是变化的，即mdiag(f1,f2,...,f(k)),函数mdiag括号里面分块阵的个数，随k的变化而变化，怎么实现，不胜感激哦！