#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
#define REP(i, n) for(int i=0; i<n; i++)
#define REPi(i, a, b) for(int i=int(a); i<int(b); i++)
template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }
template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; }
const ll MOD = 1e9+7;

int xa, ya, xb, yb;

double dist(double x1, double y1, double x2, double y2){
    return (y1-y2)*(y1-y2) + (x1-x2)*(x1-x2);
}

double solve(double p){
    double AP = dist(0, p, xa, ya);
    double BP = dist(0, p, xb, yb);
    return sqrt(AP) + sqrt(BP);
}

int main(){
    cin >> xa >> ya >> xb >> yb;

    double left = 0.0;
    double right = 1000.0;
    REP(i,100){
        double a1 = (2*left + right) / 3;
        double a2 = (left + 2*right) / 3;
        double s1 = solve(a1);
        double s2 = solve(a2);
        if(s1 < s2){
            right = a2;
        }
        else if(s2 <= s1){
            left = a1;
        }
    }

    printf("%.10lf\n", left);
    return 0;
}