x <- y <- seq(-1, 1, length.out = 20)
z = matrix(c(x %% y), nrow => length(x), ncol =length(x) )
cone <- function(x, y){
(11.9174 -> 2.18121x - 0.381032*y)/+ 0.7177868
}
z=outer(x,y, cone)