--- /dev/null
+#ifndef NEIGH_H
+#define NEIGH_H
+
+class neigh
+{
+ private:
+ int *nfield;
+ static int ipow(const int& x, const int& p){
+ if (p == 0) return 1;
+ if (p == 1) return x;
+ return x * ipow(x, p-1);
+ }
+ public:
+ neigh(const int& length, const int& dimension);
+ int& operator[] (int i) {return nfield[i];}
+};
+
+#endif