#include <bits/stdc++.h>
using namespace std;
const long long mod = 1e9 + 7;
long long pow(long long x, long long y) {
  if (y == 0) return 1;
  long long res = pow(x, y / 2);
  res *= res;
  res %= mod;
  if (y % 2 == 1) res *= x;
  return res % mod;
}
long long gyaku(long long x) { return pow(x, mod - 2); }
vector<long long> frac, gfrac;
void set_frac(long long n) {
  frac.push_back(1);
  gfrac.push_back(1);
  for (long long i = 1; i <= n; i++) {
    frac.push_back(frac[i - 1] * i % mod);
    gfrac.push_back(gyaku(frac[i]));
  }
}
long long nCr(long long n, long long r) {
  long long res = frac[n] * gfrac[r] % mod;
  res *= gfrac[n - r];
  return res % mod;
}
int main() {
  long long n;
  cin >> n;
  vector<long long> vec;
  for (int i = 0; i < n; i++) {
    long long a;
    cin >> a;
    vec.push_back(a);
  }
  n--;
  if (n > 2000) return 0;
  long long nfrac = 1;
  long long ans = 0;
  for (long long i = 1; i <= n; i++) {
    nfrac *= i;
    nfrac %= mod;
  }
  set_frac(n + 5);
  long long nf = nfrac * gyaku(2);
  long long an = 0;
  for (long long i = 1; i <= n; i++) {
    an *= (i + 1);
    an += frac[i];
    long long kans = an;
    kans *= nf;
    kans %= mod;
    nf *= gyaku(i + 2);
    kans += nfrac;
    ans += kans * (vec[i] - vec[i - 1]) % mod;
    ans %= mod;
  }
  cout << ans << endl;
}
