有如下简单的程序，R老是提示出错，但我找不到错在哪里，请高手指教。
(不用循环，除了outer 函数外还有没有其他的方法？？？)

rz=diag(1:3)
funs=function(x,y){x^2+y^2+sum(rz-x)}
a=1:4
outer(a,a,funs)