vecnorm<-function(x, p=2){
mynorm<-sum(abs(x)^p)^(1/p)
return(mynorm)
}