#include<iostream>
#include<string>
#include<cstdio>
#include<vector>
#include<cmath>
#include<algorithm>
#include<functional>
#include<iomanip>
#include<queue>
#include<ciso646>
#include<random>
#include<map>
#include<set>
#include<complex>
#include<bitset>
#include<stack>
#include<unordered_map>
#include<utility>
#include<tuple>
#include<cassert>
using namespace std;
typedef long long ll;
typedef unsigned int ui;
const ll mod = 1000000007;
const ll INF = (ll)1000000007 * 1000000007;
typedef pair<int, int> P;
#define stop char nyaa;cin>>nyaa;
#define rep(i,n) for(int i=0;i<n;i++)
#define per(i,n) for(int i=n-1;i>=0;i--)
#define Rep(i,sta,n) for(int i=sta;i<n;i++)
#define Per(i,sta,n) for(int i=n-1;i>=sta;i--)
#define rep1(i,n) for(int i=1;i<=n;i++)
#define per1(i,n) for(int i=n;i>=1;i--)
#define Rep1(i,sta,n) for(int i=sta;i<=n;i++)
typedef long double ld;
const ld eps = 1e-8;
const ld pi = acos(-1.0);
typedef pair<ll, ll> LP;
int dx[4]={1,-1,0,0};
int dy[4]={0,0,1,-1};

int n,N=1;
ld a[15];
ld dp[100050];
int po[15];

int inv(int S){
  int T=0;
  rep(i,n){
    if(S%3!=0){
      T+=po[i]*(3-S%3);
    }
    S/=3;
  }
  return T;
}

ld get_dp(int S){
  if(dp[S]!=-1) return dp[S];
  vector<int> used(n);
  int S_=S;
  bool flag=true;
  rep(i,n){
    used[i]=S_%3;
    if(S_%3==0) flag=false;
    S_/=3;
  }
  if(flag){
    int point=0;
    rep(i,n){
      if(used[i]==1) point+=a[i];
      else point-=a[i];
    }
    if(point>0) return 1;
    else return 0;
  }
  ld ma=0;
  rep(i,n){
    if(used[i]!=0) continue;
    ld mi=1;
    rep(j,n){
      if(used[j]!=0) continue;
      ld p=a[j]/(a[i]+a[j]-1);
      ld res=p*(1-get_dp(inv(S+po[i])))+(1-p)*get_dp(S+2*po[j]);
      mi=min(mi,res);
    }
    ma=max(ma,mi);
  }
  dp[S]=ma;
  return dp[S];
}

void solve(){
  cin >> n;
  po[0]=1;
  rep(i,n){
    po[i+1]=3*po[i];
  }
  //cout << 5 << " " << inv(5) << endl;
  rep(i,n) cin >> a[i];
  rep(i,n) N*=3;
  rep(i,N) dp[i]=-1;
  // rep(S,N){
  //   int S_=S;
  //   // rep(i,n){
  //   //   cout << S_%3 << " ";
  //   //   S_/=3;
  //   // }
  //   // int T_=inv(S);
  //   // cout << "" << endl;
  //   // rep(i,n){
  //   //   cout << T_%3 << " ";
  //   //   T_/=3;
  //   // }
  //   // cout << "" << endl;
  //   //cout << get_dp(S) << endl;
  // }
  cout << get_dp(0) << endl;
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout << fixed << setprecision(50);
    solve();
}