> fun<-function(f,a,b)
+ {
+   l1<-length(a)
+   l2<-length(b)
+   m<-matrix(NA,ncol=l1,nrow=l2)
+   rownames(m)<-as.character(b)
+   colnames(m)<-as.character(a)
+   for(i in 1:l1)
+   {
+     for(j in 1:l2)
+     {
+       m[j,i]<-f(a[i],b[j])
+     }
+   }
+   return(m)
+ }
> f<-function(a,b) 5*a+b
> a<-c(3,4,6,9,10)
> b<-c(2,4,5,6)
> fun(f,a,b)
   3  4  6  9 10
2 17 22 32 47 52
4 19 24 34 49 54
5 20 25 35 50 55
6 21 26 36 51 56

a<- 1:5
b<- 1:5
c<- expand.grid(a,b)
c$d<- 5*c$Var1+c$Var2
mat<- matrix(c$d,nrow = 5)
mat就是结果