向大家求教 作为刚入门R的新手我实在是写不出来了

比如有地理加权回归模型yi=a1(u,v)x1+a2(u,v)x2
要画出各个ai以u,v为轴关于500个data的曲面图
而我现在只能print出以各个u,v为原点的矩阵Ai=(a1i,a2i)
代码如下
i=1
> while(i<=500){
+ D<-(U-U
)^2+(V-V)^2
+ h=0.019
+ w<-exp(-(D/h)^2)
+ W<-diag(w)
+ A<-solve(t(X)%*%W%*%as.matrix(X))%*%t(X)%*%W%*%Y
+ print(A)
+ i=i+1
+ }
但是不知道怎么提出其中的a1a2,以及如何累计每次循环的结果画出由500数据点拟合的a1和a2的曲面图
还请大家指教