#include<iostream>
#include<algorithm>
#include<vector>
#include<string>
#include<set>
#include<queue>
#include<stack>
#include<bitset>
using namespace std;
int p = 1000000007;
#define int long long
#define vel vector<long long>
#define vvel vector<vel>
#define rep(i,n) for(long long i=0;i<n;i++)
#define sor(v) sort(v.begin(),v.end())
#define mmax(a,b) a=max(a,b)
#define mmin(a,b) a=min(a,b)
#define mkp make_pair
#define pin pair<int,int>
#define V vector
#define Endl endl
#define veb vector<bool>
#define sq(a) (a)*(a)
#define rev(s) reverse(s.begin(),s.end())
#define end_program(s) cout << s <<endl;return 0
int kai_size = 1000001;
vel kai(kai_size, 1);
vel ink(kai_size, 1);
vel dist;
int RE() {
	vel v(3, 2);
	return v.at(4); 
}
int ru(int a, int r) {
	if (r == 0) { return 1; }
	int ans = ru(a, r / 2);
	ans *= ans; ans %= p;
	if (r % 2 == 1) { ans *= a; }
	return ans % p;
}
int inv(int a) {
	return ru(a, p - 2);
}
void make_kai() {
	rep(i, kai_size-1) { kai[i + 1] = (kai[i] * (i + 1)) % p; }
	rep(i, kai_size) { ink[i] = inv(kai[i]); }
}
int com(int n, int r) {
	if (r < 0 || n < r) { return 0; }
	int ans = kai[n] * ink[r];
	ans %= p; ans *= ink[n - r]; ans %= p;
	return ans;
}
vel dis(int mid1, vvel &way) {
	int n = way.size();
	vel dist(n, -1); dist[mid1] = 0;
	queue<int> q;
	q.push(mid1);
	while (!q.empty()) {
		int st = q.front(); q.pop();
		rep(i, way[st].size()) {
			int to = way[st][i];
			if (dist[to] == -1) {
				dist[to] = dist[st] + 1;
				q.push(to);
			}
		}
	}
	return dist;
}
pin most_far(int now, int n, vvel &way) {
	vel dist1 = dis(now, way);
	pin ans = mkp(-1, 0);
	rep(i, n) {
		if (dist1[i] > ans.first) { ans = mkp(dist1[i], i); }
	}
	return ans;
}

int per(int a, int b) {
	int ans = a % b;
	if (ans < 0) { ans += b; }
	return ans;
}
vel uni(vel &v) {
	sor(v);
	vel ans(1, v[0]);
	for (int i = 1; i < v.size(); i++) {
		if (v[i] != v[i-1]) { ans.push_back(v[i]); }
	}
	v = ans;
	return v;
}
int s_gcd(int a, int b) {
	if (b == 0) { return a; }
	return s_gcd(b, a%b);
}
int gcd(int a, int b) {
	if (a < b) { swap(a, b); }
	return s_gcd(a, b);
}
bool is_prime(int i) {
	for (int j = 2; j*j <= i; j++) {
		if (i%j == 0) { return false; }
	}
	return true;
}
int get(int i, int j, vvel &ans) {
	int n = ans.size();
	if (i < 0 || n <= i || j < 0 || n <= j) { return 1; }
	return ans[i][j];
}
int lcm(int a, int b) {
	return a * b / gcd(a, b);
}
int msb(int a) {
	int ans = -1;
	int count = 0;
	while (a > 0) {
		if (a % 2 == 1) { ans = count; }
		a = a / 2; count++;
	}
	return ans;
}
int rui(int a) {
	int ans = 1;
	rep(i, a) { ans *= 2; }
	return ans;
}
signed main() {
	int n, m; cin >> n >> m;
	vvel dp(m, vel(n + 1)); dp[0][0] = 1;
	vvel ndp(m, vel(n + 1)); ndp[0][0] = 1;
	rep(i, m-1) {
		for (int j = 1; j < n; j++) {
			dp[i][j] %= p;
			int pl = dp[i][j];
			dp[i + 1][j + 1] += pl;
			dp[i + 1][j] += 2 * pl;
			dp[i + 1][j - 1] += pl;
			ndp[i][j] %= p;
			pl = ndp[i][j];
			ndp[i + 1][j + 1] += pl;
			ndp[i + 1][j] += 2 * pl;
			ndp[i + 1][j - 1] += pl;
		}
		ndp[i][0] %= p;
		int pl = ndp[i][0];
		ndp[i + 1][0] += 2 * pl;
		ndp[i + 1][1] += pl;
		if (i == 0) {
			dp[i][0] %= p;
			pl = dp[i][0];
			dp[i + 1][0] += 2 * pl;
			dp[i + 1][1] += pl;
		}
	}
	vel sum(m), nsum(m);
	rep(i, m) {
		rep(j, n) {
			sum[i] += dp[i][j];
			nsum[i] += ndp[i][j];
		}
		sum[i] -= dp[i][0];
		sum[i] %= p;
		nsum[i] %= p;
	}
	sum[0] = 1;
	int ans = 0;
	rep(i, m) {
		ans += (sum[i] * nsum[m-1 - i]) % p;
	}
	ans *= 4;
	cout << ans % p << endl;
	return 0;
}


Measured memory_bytes: 158848