#include <bits/stdc++.h>
using namespace std;
const long long MOD = 1e9 + 7;
const int MAX = 10000000;
template <typename T1, typename T2>
inline bool chmin(T1 &a, T2 b) {
  if (a > b) {
    a = b;
    return 1;
  }
  return 0;
}
template <typename T1, typename T2>
inline bool chmax(T1 &a, T2 b) {
  if (a < b) {
    a = b;
    return 1;
  }
  return 0;
}
template <typename T>
T gcd(T a, T b) {
  if (b == 0) return a;
  return gcd(b, a % b);
}
long long memo[1000001];
vector<long long> ans_vec;
long long func(long long N) {
  if (N == 1) return 1;
  if (memo[N] == -1) {
    return memo[N] = (N * func(N - 1)) % MOD;
  } else {
    return memo[N];
  }
}
int main() {
  cin.tie(0);
  ios::sync_with_stdio(false);
  long long N;
  cin >> N;
  ans_vec.resize(N - 1);
  memset(memo, -1, sizeof(memo));
  vector<long long> a(N);
  for (int i = 0; i < (int)(N); i++) cin >> a[i];
  if (N == 2) {
    cout << a[0] << endl;
    return 0;
  }
  vector<long long> b(N - 1);
  for (int i = 0; i < (int)(N - 1); i++) {
    b[i] = a[i + 1] - a[i];
  }
  ans_vec[0] = 2;
  ans_vec[1] = 3;
  for (int i = 0; i < (int)(N - 3); i++) {
    for (int j = 0; j < (int)(i + 2); j++) {
      ans_vec[j] *= i + 3;
      ans_vec[j] %= MOD;
    }
    ans_vec[i + 2] = ans_vec[i + 1] + func(i + 2);
  }
  long long ans = 0;
  for (int i = 0; i < (int)(N - 1); i++) {
    ans += b[i] * ans_vec[i];
    ans %= MOD;
  }
  cout << ans << endl;
}
