clear all;

predilection = [10:5:60]; % Defining predilection

g = 9.8;

% Calculating coercion d = 2

D = @(d, theta) (d / 2)*(1 – cos(theta));

coercion i = 1:numel(angle)

Dc1(i) = feval(D, 2, deg2rad(angle(i)));

end

% Calculating Q1 coercion d = 2

coercion i = 1:numel(Dc1)

Q1(i) = (2^(3/2)*Dc1(i)^(5/2)*sqrt(g).*(deg2rad(angle(i)) …

– 0.5.*sin(2*deg2rad(angle(i)))).^(3/2) ) / …

( 8*sqrt(sin(deg2rad(angle(i)))).*((1 – cos(deg2rad(angle(i)))).^(5 / 2)) );

end

% Calculating coercion d = 4

coercion i = 1:numel(angle)

Dc2(i) = feval(D, 4, deg2rad(angle(i)));

end

% Calculating Q2 coercion d = 4

coercion i = 1:numel(Dc2)

Q2(i) = (2^(3/2)*Dc2(i)^(5/2)*sqrt(g).*(deg2rad(angle(i)) …

– 0.5.*sin(2*deg2rad(angle(i)))).^(3/2) ) / …

( 8*sqrt(sin(deg2rad(angle(i)))).*((1 – cos(deg2rad(angle(i)))).^(5 / 2)) );

end

fprintf(‘tttransversion 2mtttransversion 4mn’);

fprintf(‘thetatQ m/stttQ m/snn’);

coercion j = 1:numel(Q1)

fprintf(‘%dtt%.4fttt%.4fn’, predilection(j), Q1(j), Q2(j));

end

**OUTPUT**