solve this quastion with this data and formula i send ypu 
وزن هواپیما ( W = 4,600 ) پوند
نسبت میانگین نیروی پیشران به وزن ( \left(\frac{T}{W}\right){\text{ave}} = 0.14 )
سرعت باد ( V_w = 0 ) (بدون باد)
شیب باند ( \phi = 0^\circ ) (بدون شیب)
انحراف فلپ ( \delta_f = 15^\circ )
مساحت بال ( S = 175 \text{ft}^2 )
طول بال ( b = 35 \text{ft} )
نسبت Aspect ( A = 7 )
ضریب اصطکاک زمین ( \mu_g = 0.03 )
ضریب بالابری حداکثر ( C{L_{\text{max}}} = 1.69 )
فاکتور اوزولد ( e = 0.80 )
شتاب گرانشی ( g = 32.174 \text{ft/s}^2 )
چگالی هوا در شرایط استاندارد سطح دریا ( \rho = 0.0023769 \text{slugs/ft}^3 )
:نیروی پیشران هنگام برخاست 
نیروی پیشران در (V=0)،(TV=0=2,000 پوند) 
نیروی پیشران در (65=V_R)،(TV=VR=1,200 پوند)
A 7
agave 8.7883
agV0 12.1445
agVR 6.1160
b 35
CD0 0.0620
CD_muCL 0.0572
CD_out_of_ge 0.1012
CDg 0.0839
CL_max 1.6900
CL_out_of_ge 0.8300
CLg 0.8900
CLOF 1.2779
dCDi -0.0173
dCL 0.1200
e 0.8000
g 32.2000
h 3.6000
h_screen 50
hTR 47.4389
k 0.7236
mu_g 0.0300
phi 0
qLOF 20.5697
qR 18.8199
rho 0.0024
RTR 5.7240e+03
S 175
SCL 19.7690
sigma 0.4800
SNGR 900.9520
SR 125.8400
STO 1.7820e+03
STR 735.4132
T_D_VLOF 591
t_rotate 1
theta_CL 0.1288
tNGR 14.3190
tTCL 0.1471
tTO 21.0560
tTR 5.5899
TV_V0 2000
TV_VR 1200
V50 137.2800
VLOF 131.5600
VLOF_VS_ratio 1.1500
VR 125.8400
VS 114.4000
Vw 0
W 4600
"The ground speed \left(V_G\right), which is the time rate of change of ground distance (s), is illustrated in Figure 9.6. It is defined as follows:
VG=ds dt=V-Vw#(9.1)
 to provide a reasonably accurate representation of jet engine thrust at low speeds:
FNFN0=1-k1V+k2V2#(9.2)
From the definitions of the lift and drag coefficients, given by Equations 7.1 and 7.5, respectively, the lift and drag forces can be expressed as follows:
L=qSCL=12ρSCLV2#(9.3)D=qSCDg=12ρSCDgV2#(9.4)
With the slats and flaps set for takeoff and the landing gear extended, the zero-lift drag coefficient can be determined from the following expression:
CD0=CD0clean +ΔCD0slats +ΔCD0flaps +ΔCD0gear #(9.5)
C_{D_g}=C_{D_0}+\lambda KC_L^2
\mu=\frac{F_{fr}}{R}
The total friction force that acts on the airplane can be expressed as follows:
Ffr=μRMW+μRNW#(9.8)
Ffr=μWcos⁡γG-L#(9.9)
0.0.8. Net Acceleration Force
For the sake of brevity, the net acceleration force acting on the airplane in the direction of motion is herein denoted as \sum F, where
 F=FN-D-Ffr-Wsin⁡γG#(9.10)
or {\ }\sum F=F_N-D-\mu\left(W\cos\funcapply\gamma_G-L\right)-W\sin\funcapply\gamma_G{\ } (from Equation 9.9)
As the runway gradient is a very small quantity, \ ^6 the small angle approximations \cos\funcapply\gamma_G\cong1 and \sin\funcapply\gamma_G\cong\gamma_G (for \gamma_G measured in radians) can be used without loss of accuracy.
Thus \sum F=F_N-D-\mu(W-L)-W\gamma_G
Substituting Equations 9.3 and 9.4 into Equation 9.12 results in the required expression for the net acceleration force:
 F=FN-12ρV2SCDg-μCL-μ+γGW#(9.13)
 The airplane's acceleration (with respect to the ground) can be expressed as follows:
a=ddtds dt=ddtV-Vw=dV dt  if dVw/dt=0#(9.14)
a=\frac{dV}{\mathrm{\ }dt}=\frac{dV}{\mathrm{\ }ds}\frac{ds}{\mathrm{\ }dt}=\frac{dV}{\mathrm{\ }ds}\left(V-V_w\right)
The ground distance measured from the start of the takeoff to the point of rotation \left(s_R\right) can be determined from the following expression:
sR=start end   ds=VwVR  V-Vwa dV#(9.16)
At any instant during the ground run, the airplane's acceleration can be obtained from Newton's second law (see Section 2.4.5):
a= Fm=g FW#(9.17)
s_R=\int_{V_w}^{V_R}\hairsp\hairsp\frac{W\left(V-V_w\right)}{g\sum\ F}\mathrm{\ }dV=\int_{V_w}^{V_R}\hairsp\hairsp\frac{WV}{g\sum\ F}\mathrm{\ }dV-\int_{V_w}^{V_R}\hairsp\hairsp\frac{WV_w}{g\sum\ F}\mathrm{\ }dV
The mean thrust FN is calculated at a reference speed, herein denoted as V^\ast, which is determined as follows:
V*2=12VR-Vw2#(9.19) or  V*=12VR-Vw=0.707VR-Vw#(9.20)
An expression for the net accelerating force ( \sum F ) can be obtained by substituting Equation 9.6 into Equation 9.13, that is,
 F=FN-12ρV2SCD0+λKCL2-μCL-μ+γGW#(9.21)
The required distance \left(s_R\right) is obtained by a summation of the two integrals that are given in Equation 9.18. For the sake of convenience, the two integral expressions will be evaluated separately, where
1  =VwVR  WVg F dV  and  2  =VwVR  WVwg F dV#(9.22)
For the first integral, the following variable substitution is made:
f=V2#(9.23)
thus df=2V\mathrm{\ }dV and dV=\frac{1}{2V}\mathrm{\ }df
hence \int_1\hairsp=\int_{V_w}^{V_R}\hairsp\frac{WV}{g\sum F}\mathrm{\ }dV=\int_{\sqrt{f_w}}^{\sqrt{f_R}}\hairsp\frac{WV}{g\sum F}\left(\frac{1}{2V}\right)df=\frac{1}{2g}\int_{\sqrt{f_w}}^{\sqrt{f_R}}\hairsp\frac{W}{\sum F}\mathrm{\ }df
Now substitute \sum F from Equation 9.21 into Equation 9.24:
1  =12gfwfR  1FNW-μ-γG-ρS2WCD0+λKCL2-μCLf df#(9.25)
To simplify the analysis, arbitrary constants A_2,B_2, and C_2 are now introduced, where
A_2=FNW-μ-γG  and  B2=ρS2WCD0+λKCL2-μCL
and {\ }C_2=\sqrt{\frac{A_2}{B_2}}
Thus \int_1\hairsp=\frac{1}{2g}\int_{\sqrt{f_w}}^{\sqrt{f_R}}\hairsp\frac{1}{A_2-B_2f}\mathrm{\ }df=\frac{-1}{2gB_2}\left[\ln\funcapply\left(A_2-B_2f\right)\right]_{\sqrt{f_w}}^{\sqrt{f_R}}
1  =-12gB2ln⁡A2-B2V2VwVR=12gB2ln⁡A2-B2Vw2A2-B2VR2#(9.27)
The second integral is now considered:
\begin{matrix}\ &\ \int_{2}\hairsp\hairsp=\frac{V_w}{g}\int_{V_w}^{V_R}\hairsp\hairsp\frac{1}{A_2-B_2V^2}\mathrm{\ } dV=\frac{V_W}{gB_2}\int_{V_w}^{V_R}\hairsp\hairsp\frac{1}{C_2^2-V^2}\mathrm{\ } dV\\\end{matrix}
(\int_{2}\ =\frac{V_w}{2gB_2C_2}\left[\ln{\left(\frac{C_2+V}{C_2-V}\right)}\right]_{V_w}^{V_R}=\frac{V_w}{2gB_2C_2}\ln{\left\{\frac{\left(C_2+V_R\right)\left(C_2-V_w\right)}{\left(C_2-V_R\right)\left(C_2+V_w\right)}\right\}})
The two results are now added to give the distance s_R :
sR=12gB2ln⁡A2-B2Vw2A2-B2VR2-VwC2ln⁡C2+VRC2-VwC2-VRC2+Vw#(9.29)
where A_2,B_2, and C_2 are defined in Equation 9.26 .
uniform acceleration:
sR=1aVwVR  V-VwdV=1aV22-VVwVwVR=12aVR-Vw2#(9.30)
the following integral expression is obtained:
sR=VwVR  V-Vwa dV=WgVwVR  V-Vw F dV#(9.31)
To integrate Equation 9.31 numerically, it is possible to divide the takeoff distance s_R into n segments (or velocity intervals) and to use the trapezoidal rule for incremental values of V :
sR=W2gi=0n  (V)i+1-Vw Fi+1+(V)i-Vw FiΔV#(9.32)
The time that the airplane takes to accelerate to the rotation speed \(\left(V_R\right)\) may be determined using similar approaches to those presented for the distance calculation (Sections 9.4.1 to 9.4.4). The time-velocity relationship is given by
\[
\mathrm{d}t=\frac{1}{a}\mathrm{d}V\emsp\left(\mathrm{\ based\ on\ Equation\ }9.15\right)
Thus\ t=\int_{\mathrm{start\ }}^{\mathrm{end\ }}\mathrm{d}t=\int_{V_w}^{V_R}\frac{1}{a}\mathrm{d}V)
t_R=\frac{1}{\bar{a}}\int_{V_w}^{V_R}\mathrm{d}V=\frac{V_R-V_w}{\bar{a}}
s_{RL}=\Delta t_{RL}\left(\frac{V_R+V_{LOF}}{2}-V_w\right)
s_g=s_R+s_{RL}
U\cong\int\left(F_N-D\right)\mathrm{d} s=\left(F_N-D\right)_{\mathrm{avg\ }}s_{L35}
U=\frac{W}{2g}\left(V_{35}^2-V_{LOF}^2\right)+h_{sc}W
By combining Equations 9.37 and 9.38, the following expression for the distance from liftoff to the 35 ft point is obtained (for zero wind):
s_{L35}=\left(\frac{W}{\left(F_N-D\right)_{avg}}\right)\left(\frac{V_{35}^2-V_{LOF}^2}{2g}+h_{sc}\right)
An alternative-and relatively straightforward-method to estimate the distance from liftoff to the 35 ft point is based on the product of the average time and the average ground speed for this segment of the takeoff:
s_{L35}=\Delta t_{L35}\left(\frac{V_{LOF}+V_{35}}{2}-V_W\right)
s35=V3522a#(9.41)
where s_{35} is the distance from brake release (V=0) to the 35 ft screen height (typical units: m,ft ); and
V_{35} is the TAS at 35 ft (equal to the ground speed as V_w=0 ) (typical units: m/s,ft/s ).
As the greatest part of the takeoff distance will be the ground run, Equation 9.21 can be used to derive a relationship for the mean acceleration. Assuming zero slope (i.e., \gamma_G=0 ), it follows that
a= Fm=gWFN-12ρV2SCD-μRCL-μRW#(9.42)
a=gFNW(1-k)
s_{35}=\frac{V_{35}^2W}{2g(1-k)F_{N_{0.7}}}
where 𝐹𝑁0.7 is the thrust evaluated at 0.7𝑉35 (typical units: N,lb ).
The speed term in Equation 9.44 is now eliminated using the lift coefficient at the screen height 𝐶𝐿35, that is,
$$
𝐶𝐿35=2𝐿𝜎𝜌0𝑉352𝑆=2𝑛𝑧𝑊𝜎𝜌0𝑉352𝑆#(9.45)

"