n,m = map(int,input().split()) dp = [[0]*(m+1) for _ in range(n+1)] dp[1][1] = 1 mod = 10**9+7 for i in range(2,n+1): dp[i][1] = dp[i-1][1] % mod for j in range(2,m+1): dp[1][j] = dp[1][j-1] % mod for i in range(2,n+1): for j in range(2,m+1): dp[i][j] = (dp[i][j-1] + dp[i-1][j]) % mod print(dp[n][m])