#include <bits/stdc++.h>
#pragma GCC optimize("Ofast")
using namespace std;
const int maxN = 1 << 11, mod = 1000000007;
template <class T1, class T2>
inline void addMod(T1& a, T2 b) {
  a = (a + b) % mod;
}
template <class T1, class T2>
inline void multMod(T1& a, T2 b) {
  a = a * b % mod;
}
long long qpow(long long a, long long b) {
  long long res0 = 1;
  for (; b != 0; b >>= 1) {
    if (b & 1) multMod(res0, a);
    multMod(a, a);
  }
  return res0;
}
long long fact[maxN], tcaf[maxN];
long long cf(int n, int k) {
  bool static initd = false;
  if (!initd) {
    fact[0] = 1;
    for (int i = (1); i < (int)(maxN); i++) fact[i] = fact[i - 1] * i % mod;
    tcaf[maxN - 1] = qpow(fact[maxN - 1], mod - 2);
    for (int i = maxN - 1; i > 0; i--) tcaf[i - 1] = tcaf[i] * i % mod;
    initd = true;
  }
  if (k < 0 or k > n) return 0;
  long long res0 = fact[n];
  multMod(res0, tcaf[k]);
  multMod(res0, tcaf[n - k]);
  return res0;
}
long long f(int n, int d) { return cf(n + 1, d + 1); }
int x[maxN];
int main() {
  int n;
  scanf("%d", &n);
  cf(n, 0);
  for (int i = (0); i < (int)(n); i++) scanf("%d", x + i);
  int n0 = n--;
  long long res = 0;
  for (int i = (1); i < (int)(n); i++)
    for (int j = (0); j < (int)(i); j++) {
      int d = i - j - 1;
      long long tak = f(n - 1, d + 1);
      multMod(tak, d + 1);
      multMod(tak, mod - 1);
      addMod(tak, f(n - 2, d) * n);
      multMod(tak, fact[d]);
      if (i != n0 - 1) multMod(tak, fact[n - d - 2]);
      multMod(tak, x[i] - x[j]);
      addMod(res, tak);
    }
  int i = n0 - 1;
  for (int j = (0); j < (int)(i); j++) {
    int d = i - j - 1;
    long long tak = f(n - 1, d);
    multMod(tak, fact[d]);
    multMod(tak, fact[n - d - 1]);
    multMod(tak, x[i] - x[j]);
    addMod(res, tak);
  }
  printf("%lld\n", res);
  return 0;
}
