Write a Regula Falsi program in Matlab to find the angle of the projectile motio

Question

Write a Regula Falsi program in Matlab to find the angle of the projectile motion equation. Let the program start at [0,1.6]

Examen #01 (ME2020) Fecha: M, 2018 Nombre 1. El lanzamiento de un balón desde un jugador (el lanzador) hasta otro es modelado matemáticamente por las ecuaciones del movimiento parabólico, una de las ecuaciones que permite hallar la altura hg en la dirección del eje "y es la siguiente: 1 x2g 2 v cos2(0) hg + 3 x tan(0) _ + he+ 2 donde hg es la áltura del jugador que recibe el batón, ho es la altura.bel lanzador, vo es la velocidad inicial del balón al dejar las manos del lazador, x es la distancia horizontal que tiene que recorrer el balón para llegar al otro jugador, g es la gravedad, y 0 es el ángulo que se mide desde una linea horizontal justo cuando el balón es lanzado de la mano del jugador. Si g 32.2 pies/s2, ho 6.5 pies, hg -7 pies, x 60 pies y vo 58 pies/s usar Matlab para hacer un programa que permita hallar el ángulo por los siguentes métodos (empezar con el intervalo [0, 1.6l) (b) (a) El método de la posición falsa (80 puntos) (d) El método de la secante (10 puntos) () Usando la función tero de Matlab (10 puntos) FAXFb 0 Signs 014 tacan u.) (BONO) Hacer un diagrama de flujo para el método de la secante (5 puntos)

Explanation / Answer

(False Position or Regula Falsi Method). To approxi-
mate a root of the equation f(x) = 0 in the interval [a, b]. Proceed with
the method only if f(x) is continuous and f(a) and f(b) have opposite signs.

Put the values of a = 0 and b = 1.6

a as left and b as right end points
function [c, err, yc] = regula (f, a, b, delta, epsilon, maxl)
%Input - f is the function input as a string ’f’
% - a and b are the left and right end points
% - delta is the tolerance for the zero
% - epsilon is the tolerance for the value of f at the zero
% - max1 is the maximum number of iterations
%Output - c is the zero
% - yc = f(c)
% - err is the error estimate for c
ya=feva1(f,a);
yb=feva1(f,a);
if ya*yb>0
disp (’Note: f(a)*f(b)>0’),
break,
end
for k=1:max1
dx=yb*(b-a)/(yb-ba);
c=b-dx;
ac=c-a;
yc=feval(f,c);
if yc==0, break;
elseif yb*yc>0
b=c;
yb=yc;
end
dx=min (abs (dx), ac);
if abs (dx)<delta, break, end
if abs (yc)<epsilon, break, end
end
c;
err=abs (b-a)/2;
yc=feval(f,c);

Related Questions

Navigate

• Browse (All)
• Browse W
• Subjects

---

---

Integrity-first tutoring: explanations and feedback only — we do not complete graded work. Learn more.