地理加权回归的3D图绘制问题

向大家求教 作为刚入门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的曲面图
还请大家指教