Scribd
Upload
263 views

Example 3 - Bisection Method

This document describes using the bisection method to determine the drag coefficient (c) of a parachute design within 0.001% accuracy. The method iteratively calculates the function f(c) at values between an initial upper (xU) and lower (xL) bound, then selects the midpoint (xR) of the range as the new bounds. Over 14 iterations, the range is narrowed from 2 to 0.002, converging on a drag coefficient of 14.8011.

Uploaded by

Ryan Irfandi
Copyright
© © All Rights Reserved
Available Formats
Download as XLS, PDF, TXT or read online on Scribd
263 views

Example 3 - Bisection Method

This document describes using the bisection method to determine the drag coefficient (c) of a parachute design within 0.001% accuracy. The method iteratively calculates the function f(c) at values between an initial upper (xU) and lower (xL) bound, then selects the midpoint (xR) of the range as the new bounds. Over 14 iterations, the range is narrowed from 2 to 0.002, converging on a drag coefficient of 14.8011.

Uploaded by

Ryan Irfandi
Copyright
© © All Rights Reserved
Available Formats
Download as XLS, PDF, TXT or read online on Scribd
You are on page 1/ 2

Parachute Design (using the bisection method for finding roots)

ct

g m
1-e m - 40
f (c) =

Mass, m
Velocity, v
Gravity, g
Time, t

=
=
=
=

68.100
40.000
9.810
10.000

kg
m/sec
m/sec2
sec

Determine the drag coefficient, c, within 0.001%.


Bisection Method
Iteration
1
2

xL
14
14

f(xL)
1.611116
1.611116

xU
16
15

f(xU)
-2.23026
-0.38446

xRNEW - xROLD
ea1 =
100 (%)
NEW
xR

xR
15
14.5

e a2

f(xR)
-0.38446
0.593698

ea1

ea2

3.4483

6.6667
3.4483

xU - xL
=
100(%)
xU + xL

Parachute Design (using the bisection method for finding roots)


f (c) =

Mass, m
Velocity, v
Gravity, g
Time, t

=
=
=
=

ct

g m
1-e m - 40

68.100
40.000
9.810
10.000

kg
m/sec
m/sec2
sec

Determine the drag coefficient, c, within 0.001%.


Bisection Method
Iteration
1
2
3
4
5
6
7
8
9
10
11
12
13
14

xL
14
14
14.5
14.75
14.75
14.75
14.78125
14.7969
14.7969
14.8008
14.8008
14.8008
14.8008
14.8010

f(xL)
1.611116
1.611116
0.593698
0.09983
0.09983
0.09983
0.038775
0.008303
0.008303
0.000691
0.000691
0.000691
0.000691
0.000215

xU
16
15
15
15
14.8750
14.8125
14.8125
14.8125
14.8047
14.8047
14.8027
14.8018
14.8013
14.8013

f(xU)
-2.23026
-0.38446
-0.38446
-0.38446
-0.1435
-0.02213
-0.02213
-0.02213
-0.00692
-0.00692
-0.00311
-0.00121
-0.00026
-0.00026

xRNEW - xROLD
ea1 =
100 (%)
NEW
xR

xR
15
14.5
14.75
14.875
14.8125
14.78125
14.79688
14.80469
14.80078
14.80273
14.8018
14.8013
14.8010
14.8011

e a2

f(xR)
-0.38446
0.593698
0.09983
-0.1435
-0.02213
0.038775
0.008303
-0.00692
0.000691
-0.00311
-0.00121
-0.00026
0.000215
-2.2E-05

ea1

ea2

3.4483
1.6949
0.8403
0.4219
0.2114
0.1056
0.0528
0.0264
0.0132
0.0066
0.0033
0.0016
0.0008

6.6667
3.4483
1.6949
0.8403
0.4219
0.2114
0.1056
0.0528
0.0264
0.0132
0.0066
0.0033
0.0016
0.0008

xU - xL
=
100(%)
xU + xL

You might also like