#include <algorithm>
#include <array>
#include <cmath>
#include <cstdint>
#include <iostream>
#include <vector>
using namespace std;
using point = array<int, 2>;
int64_t operator^(point a, point b) { return a[0] * (int64_t)b[1] - a[1] * (int64_t)b[0]; }
int64_t operator&(point a, point b) { return a[0] * (int64_t)b[0] + a[1] * (int64_t)b[1]; }
point operator-(point a) { return {-a[0], -a[1]}; }
point operator-(point a, point b) { return {a[0] - b[0], a[1] - b[1]}; }
vector<point> env(const vector<point>& v) {
vector<point> e;
for (auto& p: v) {
while (e.size() > 1 && (e.back() - e.end()[-2] ^ p - e.back()) >= 0) e.pop_back();
e.push_back(p);
}
return e;
}
int main() {
cin.tie(0)->sync_with_stdio(0);
int n; cin >> n;
vector<point> v[2];
v[0].resize(n), v[1].resize(n);
for (int i = 0; i < n; ++i) cin >> v[0][i][1] >> v[1][i][1] >> v[0][i][0] >> v[1][i][0];
vector<point> e[4];
for (int c: {0, 1}) {
sort(v[c].begin(), v[c].end());
e[c * 2] = env(v[c]);
for (auto& p: v[c]) p = -p;
reverse(v[c].begin(), v[c].end());
e[c * 2 + 1] = env(v[c]);
}
auto ans = 0.L / 0;
int i[4]{};
while (1) {
int w = -1;
for (int j = 0; j < 4; ++j) if (i[j] + 1 < e[j].size() && (!~w || (e[j][i[j] + 1] - e[j][i[j]] ^ e[w][i[w] + 1] - e[w][i[w]]) < 0)) w = j;
point a = e[0][i[0]] - -e[1][i[1]], b = e[2][i[2]] - -e[3][i[3]], d;
if (~w) d = e[w][i[w] + 1] - e[w][i[w]], ++i[w];
if ((!~w || d[1] < 0) && isnan(ans)) ans = a[1] * 1.L * b[1];
if (!~w) break;
if (d[1] < 0) ans = min(ans, (d ^ a) * 1.L * (d ^ b) / d[0] / d[0]);
}
assert(isfinite(ans) && ans >= 0);
cout.precision(20);
cout << fixed << ans;
}