/*    特徴    */

/*    教訓    */

#include <bits/stdc++.h>
using namespace std;

using ll = long long;
using ld  = long double;
using pint = pair<int, int>;
using pll = pair<ll, ll>;
using pld = pair<ld, ld>;
const int INF=1e9+7;
const ll LINF=9223372036854775807;
const ll MOD=1e9+7;
const ld PI=acos(-1);
const ld EPS = 1e-9; //微調整用（EPSより小さいと0と判定など）

//MODINT
//変数名.nでint型の数値を受け取る
struct mint {
  int n;
  mint(int n_ = 0) : n(n_) {}
};

mint operator+(mint a, mint b) { a.n += b.n; if (a.n >= MOD) a.n -= MOD; return a; }
mint operator-(mint a, mint b) { a.n -= b.n; if (a.n < 0) a.n += MOD; return a; }
mint operator*(mint a, mint b) { return (long long)a.n * b.n % MOD; }
mint &operator+=(mint &a, mint b) { return a = a + b; }
mint &operator-=(mint &a, mint b) { return a = a - b; }
mint &operator*=(mint &a, mint b) { return a = a * b; }

int ii() { int x; if (scanf("%d", &x) == 1) return x; }
long long il() { long long x; if (scanf("%lld", &x) == 1) return x; }
string is() { string x; cin >> x; return x; }
char ic() { char x; cin >> x; return x; }
void oi(int x) { printf("%d ", x); }
void ol(long long x) { printf("%lld ", x); }
void od_nosp(double x) { printf("%.15f", x); } // 古い問題用
void od(double x) { printf("%.15f ", x); }
// long doubleで受け取り、fをLfなどに変えて出力すると、変な数値が出る
// それをなんとかするには独自の出力を作らなければならなそう
void os(const string &s) { printf("%s ", s.c_str()); }
void oc(const char &c) { printf("%c ", c); }
#define o_map(v){cerr << #v << endl; for(const auto& xxx: v){cout << xxx.first << " " << xxx.second << "\n";}} //動作未確認
void br() { putchar('\n'); }

#define gcd __gcd //llは受け取ってくれない
int lcm(int a, int b){return a / gcd(a, b) * b;}
#define begin_end(a) a.begin(),a.end() //sort(begin_end(vec));
#define REP(i,m,n) for(ll i=(ll)(m) ; i < (ll) (n) ; i++ )
#define rep(i,n) REP(i,0,n)
#define m_p(a,b) make_pair(a,b)
#define SORT_UNIQUE(c) (sort(c.begin(),c.end()), c.resize(distance(c.begin(),unique(c.begin(),c.end()))))
#define GET_POS(c,x) (lower_bound(c.begin(),c.end(),x)-c.begin())
#define p_b push_back
#define SZ(x) ((int)(x).size) //size()がunsignedなのでエラー避けに

//4近傍（上下左右） rep(i, 2) にすると右・下だけに進む
vector<int> dx_4 = {1, 0, -1, 0};
vector<int> dy_4 = {0, 1, 0, -1};

// coutによるpairの出力（空白区切り）
template<typename T1, typename T2> ostream& operator<<(ostream& s, const pair<T1, T2>& p) {return s << "(" << p.first << " " << p.second << ")";}
// 空白区切りだけ --> return s << "(" << p.first << " " << p.second << ")";
// 見やすくしたいとき --> return s << "(" << p.first << " " << p.second << ")";
// coutによるvectorの出力（空白区切り）
template<typename T> ostream& operator<<(ostream& s, const vector<T>& v) {
  int len = v.size();
  for (int i = 0; i < len; ++i) {
    s << v[i]; if (i < len - 1) s << " "; //"\t"に変えるとTabで見やすく区切る
  }
  return s;
}
// coutによる多次元vectorの出力（空白区切り）
template<typename T> ostream& operator<<(ostream& s, const vector< vector<T> >& vv) {
  int len = vv.size();
  for (int i = 0; i < len; ++i) {
    s << vv[i] << endl;
  }
  return s;
}

//最大値、最小値を更新する。aよりbのが大きい（小さい）か等しければaを更新してtrueを返す。そうでなければ何もせずfalseを返す chmax(nowmax,x);
template<typename T>
bool chmax(T& a, T b){return (a = max(a, b)) == b;}
template<typename T>
bool chmin(T& a, T b){return (a = min(a, b)) == b;}

// -------- template end - //


// - library ------------- //

// 決め打ち二分探索 with めぐる式二分探索

// 答えを key としたときに、問題が求める条件を満たせるかどうか
bool isOK(ll key, vector<ll>& S) {
  // 可否を評価（パターン色々。評価配列を作り直して貪欲に落とし込むことが多い？）
  // okを適当に大きめの値に設定した場合、配列外参照やコーナーケースに注意

  // key文字の連続部分列で条件を満たせるか
  ll len = S.size();
  vector<ll> bucket(3, 0);
  // 1,2,3 を0-indexedに
  bool flag = false;
  rep(i, len){
    bucket[S[i]]++;
    if (bucket[0]>=1 && bucket[1]>=1 && bucket[2]>=1){
      flag = true;
      break;
    }
    else {
      // iを含めてkey-1個分残す
      // ex)i=5, key=4なら、3,4,5を残す
      // つまり i-key+1 が0以上なら削除
      if (i-key+1>=0) bucket[S[i-key+1]]--;
    }
  }

  if ( flag ){
    return true;
  }
  else{
    return false;
  }

}

// 汎用的な二分探索のテンプレ
ll BinarySearch(ll N, vector<ll>& S) {
  ll ng = 2; // 1, 2はありえない
  ll ok = N + 1;
  // okは"1つも満たさない"場合があることに注意して設定。ex) = (ll)vec.size();
  // 二分探索はO(logN)なので大きめでもいいが、配列外参照やコーナーケースは注意

  // ok と ng のどちらが大きいかわからないことを考慮
  while (abs(ok - ng) > 1) {
    ll mid = (ok + ng) / 2;

    if (isOK(mid, S)) ok = mid;
    else ng = mid;
  }
  return ok;
}

// 度数法（degree）から弧度法（radian）へ変換
ld to_rad(ld deg){
  return deg * PI / 180;
}

// 座標の回転
// 座標(x, y) を 座標(0, 0)を中心にして右向きに d 度（度数法）回転した座標
// to_rad() も貼ること！！
pair<ld, ld> rotate_xy(ld x, ld y, ld d){
  pair<ld, ld> new_xy;
  new_xy.first = x * cos(to_rad(d)) - y * sin(to_rad(d));
  new_xy.second = x * sin(to_rad(d)) + y * cos(to_rad(d));
  return new_xy;
}

// --------- library end - //



// - main() -------------- //

int main(){


  ll A, B, H, M;
  cin >> A >> B >> H >> M;
  ld a, b;
  a = ld(A);
  b = ld(B);
  
  ld hourx = 0;
  ld houry = a;
  ld minux = 0;
  ld minuy = b;

  ld hour_rot = ld(0.5 * M + (ld)H * 30);
  ld minu_rot = ld(M) * ld(6);

  pair<ld, ld> new_hourxy = rotate_xy(hourx, houry, hour_rot);
  pair<ld, ld> new_minuxy = rotate_xy(minux, minuy, minu_rot);

  od(hypot(new_hourxy.first - new_minuxy.first, new_hourxy.second - new_minuxy.second));



// ---------- main() end - //
}

Measured memory_bytes: 4056