Scribd
Upload
14 views2 pages

Bisection Method Code Sunil

Uploaded by

Prathamesh Sr
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
14 views2 pages

Bisection Method Code Sunil

Uploaded by

Prathamesh Sr
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

SUNIL KUMAR | P23ST013 | BISECTION METHOD

F = @ (x) (x^3-6*x^2+11*x-6) ;
a = 0; b = 1.5; imax = 20; tol = 0.0001 ;
Fa = F(a); Fb = F(b);
if Fa*Fb > 0
disp('Error:The function has the same sign at points a and b.')
else
disp('iteration a b (xNS)Solution f(xNS) Tolerance')
for i = 1:imax
xNS = (a + b)/2;
toli = (b - a )/2;
FxNS = F(xNS);
fprintf('%3i %14.6f %11.6f %11.6f %11.6f %11.6f\n', i, a, b, xNS, FxNS,
toli)
if FxNS == 0
fprintf('An exact solution x=%11.6f was found',xNS)
break
end
if toli < tol
break
end
if i == imax
fprintf('Solution was not obtained in %i iterations',imax)
break
end
if F(a)*FxNS < 0
b= xNS;
else
a = xNS;
end
end
end

RESULT:

iteration a b (xNS)Solution f(xNS) Tolerance


1 0.000000 1.500000 0.750000 -0.703125 0.75000
2 0.750000 1.500000 1.125000 0.205078 0.375000
3 0.750000 1.125000 0.937500 -0.136963 0.187500
4 0.937500 1.125000 1.031250 0.059601 0.093750
5 0.937500 1.031250 0.984375 -0.031986 0.046875
6 0.984375 1.031250 1.007812 0.015442 0.023438
7 0.984375 1.007812 0.996094 -0.007858 0.011719
8 0.996094 1.007812 1.001953 0.003895 0.005859
9 0.996094 1.001953 0.999023 -0.001956 0.002930
10 0.999023 1.001953 1.000488 0.000976 0.001465
11 0.999023 1.000488 0.999756 -0.000488 0.000732
12 0.999756 1.000488 1.000122 0.000244 0.000366
13 0.999756 1.000122 0.999939 -0.000122 0.000183
14 0.999939 1.000122 1.000031 0.000061 0.000092
SUNIL KUMAR | P23ST013 | PROJECTILE MOTION
v = 120; % v = initial velocity in m/s
angle = 30 ; % angle = trajectory angle in degree
r = angle*pi/180; % converts degree to radian
g = 10; % g = acceleration due to gravity
T = 2*v*sin(r)/g ; % time of flight in seconds
fprintf('Total time of flight T= %1.2f\n',T);
disp('time(t) X Y')
for t = 0:0.5:T % t = time in seconds
Vi = v*cos(r); % velocity in X-direction
Vj = v*sin(r); % velocity in Y-direction
X = (Vi*t); % X(t)
Y = (Vj*t-0.5*g*t^2); % Y(t)
fprintf('%3.2f %11.2f %11.2f\n',t , X, Y);
end

RESULT:

Total time of flight T= 12.00


time(t) X Y
0.00 0.00 0.00
0.50 51.96 28.75
1.00 103.92 55.00
1.50 155.88 78.75
2.00 207.85 100.00
2.50 259.81 118.75
3.00 311.77 135.00
3.50 363.73 148.75
4.00 415.69 160.00
4.50 467.65 168.75
5.00 519.62 175.00
5.50 571.58 178.75
6.00 623.54 180.00
6.50 675.50 178.75
7.00 727.46 175.00
7.50 779.42 168.75
8.00 831.38 160.00
8.50 883.35 148.75
9.00 935.31 135.00
9.50 987.27 118.75
10.00 1039.23 100.00
10.50 1091.19 78.75
11.00 1143.15 55.00
11.50 1195.12 28.75
12.00 1247.08 0.00

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy