blob: 312e81693c5b103bc460b8c5e6ca33493fa3baf3 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
% Define problem constants
g = 9.81;
mu = 0.55;
F = 150;
m = 25;
% Define the function to be solved for. Is the angle specified in radians or degrees?
%angle=@(-asin(g/(sqrt(mu^2)))+atand(mu)+((4*n2+1)*pi)/2
func=@(angle) mu*m*g/(cosd(angle)+mu*sind(angle))-F;
% THINK, what makes range sense for angle?
lower_bound=-90;
upper_bound=90;
% Plot your function. Does plotting give you an idea about where the root is?
hold on
fplot(func)
fplot(0)
% Finaly solve for the root using the bisection script given to you, which can be called as:
%[root, fx, ea, iter] = bisect(func, lower_bound, upper_bound);
angle = bisect(func,lower_bound, upper_bound)
% These variables will be checked for your final answer:
%angle = % the angle in degree that solves this problem
fx = func(angle) % the function value at your solved angle
ea = ea % the bisection error estimate
iter = iter % the number of bisection iterations
|
