Diferencia entre revisiones de «Método de bisección (Grupo 30)»

De MateWiki
Saltar a: navegación, buscar

Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/mat/public_html/w/includes/diff/DairikiDiff.php on line 434
(Programa)
Línea 25: Línea 25:
  
 
{{matlab|codigo=
 
{{matlab|codigo=
% Este programa dibuja la gráfica de la función f(x)=-1/2+1/4*x en el intervalo [-2,4]
+
% Este programa calcula el punto donde la función  
x=-2:0.01:4;                                     % coordenadas x de los puntos
+
f=@(x) sin(x)-(x/3);
y=-1/2+1/4*x;                                     % imágenes
+
ei=0;
figure(1)                                         % abrimos una pantalla para dibujar
+
ed=2*pi;
hold on                                          % para que no borre lo ya dibujado
+
while((ed-ei)>1.e-3)
plot(x,y)                                     % Dibuja la gráfica
+
        if (f(ei)*f((ed+ei)/2)<0)
 +
          ed =(ei+ed)/2;
 +
        else
 +
        ei =(ed+ei)/2;
 +
        end
 +
end
 +
Sol=((ed+ei)/2)
 
}}
 
}}
  

Revisión del 14:29 13 dic 2019

Trabajo realizado por estudiantes
Título Aproximación de raíces por el método de bisección. Grupo 30
Asignatura Matemáticas I
Curso Curso 2019-20
Autores Guillermo Izquierdo, Sara Guadalix, Enrique Adrados
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura


Se busca encontrar la intersección entre dos funciones

1 Planteamiento

f=@(x) sin(x)= (x/3)


2 Método

Explicamos brevemente en qué consiste el método de bisección ...

3 Aplicación

Explicamos cómo adaptar el método a nuestro problema. Decimos cual es la función, el intervalo, el error máximo que vamos a admitir, el criterio de parada, etc.

Además damos el valor de la aproximación con el error que hemos prefijado. Por ejemplo: El valor de la aproximación es ... con un error ...

4 Programa

Aquí incluimos el programa en Matlab, como en el ejemplo de abajo:

% Este programa calcula el punto donde la función 
f=@(x) sin(x)-(x/3);
ei=0;
ed=2*pi;
while((ed-ei)>1.e-3)
        if (f(ei)*f((ed+ei)/2)<0)
          ed =(ei+ed)/2;
        else
         ei =(ed+ei)/2;
        end
end
Sol=((ed+ei)/2)