1

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<iostream>
using namespace std;

long double a=-5.8522e+04,
    b=3.7519e+11,
    c=6.8944e-07,
    g0=9.7803,
    gama=1.4000,
    k=0.00046,
    P0=101325,
    P00=101325.0,
    Re=6378137.0,
    rou0=1.2250,
    rou00=1.1031865,
    T0=273.0,
    T00=303.1500,
    mu=0.027441,
    R=8.3140,
    A=58857.2530,
    B=375399623239.3531;
long double n0=pow(1+k*c*(a+b/Re),2.5);

long double phi()
{
    int n=10000;
    long double st=1,ed=1+k*rou00;
    long double delta=(ed-st)/n,ans=0;
    for(int i=1;i<=n;i++)
    {
        long double X=(st-1)*pow(T00,3.5)/k/rou0/T0;
        long double H=pow(X,0.4)+A-303.15;
        H=B/H;
        H-=6378137;
        ans+=n0*Re/st/(Re+H)/
        sqrt(1-(n0*Re/st/(Re+H))*(n0*Re/st/(Re+H)))/st*delta;
        st+=delta;
    }
    return ans;
}

long double eps=1e-8,tar=0.0098901991;

long double tdiv()
{
    long double l=0.0,r=1.0;
    while(r-l>eps)
    {
        long double mid=(l+r)/2;
        k=mid;
        long double ans=phi();
        if(ans<tar) l=mid;
        else r=mid;
    }
    return l;
}

int main()
{
    long double ans=phi();
    tdiv();
    return 0;
}