#include <bits/stdc++.h>
using namespace std;
double dp[1005][5005];
double a[10005], b[10005], c[10005], x[1005], d[1005], c1[1005], d1[1005];
void solve(int n) {
  int i;
  for (i = 1; i <= n; i++) {
    if (i == 1) {
      c1[i] = c[i] / b[i];
      d1[i] = d[i] / b[i];
    } else {
      c1[i] = c[i] / (b[i] - a[i] * c1[i - 1]);
      d1[i] = (d[i] - a[i] * d1[i - 1]) / (b[i] - a[i] * c1[i - 1]);
    }
  }
  for (i = n; i >= 1; i--) {
    if (i == n) {
      x[i] = d1[i];
    } else {
      x[i] = d1[i] - c1[i] * x[i + 1];
    }
  }
}
int main() {
  int n, m, i, j, a, b;
  cin >> n >> m;
  cin >> a >> b;
  if (m == 1) {
    printf("%d\n", 2 * (n - a));
    return 0;
  }
  if (a == n)
    printf("0\n");
  else {
    for (i = n - 1; i >= a; i--) {
      for (int f = 0; f < 100; f++) {
        for (j = 1; j <= m; j++) {
          if (j == 1)
            dp[i][j] = (dp[i][j] + dp[i + 1][j] + dp[i][j + 1]) / 3 + 1;
          else if (j == m)
            dp[i][j] = (dp[i][j] + dp[i + 1][j] + dp[i][j - 1]) / 3 + 1;
          else {
            dp[i][j] =
                (dp[i + 1][j] + dp[i][j] + dp[i][j - 1] + dp[i][j + 1]) / 4 + 1;
          }
        }
      }
    }
    printf("%.12lf\n", dp[a][b]);
  }
}
