library(raster) N <- 2500 A <- read.table("BrunelNetDelta9ML.con", sep = " ", comment.char = "", quote = "", colClasses = "numeric", header = FALSE) net <- matrix(as.matrix (A), N, N)