![4. Escribir dos funciones MATLAB para la resoluci´on de ecuaciones de segundo grado ax2 +bx+c =
0, una de ellas (que llamaremos mala.m) implementando la f´ormula:
x =
−b ± b2 − 4ac
2a
y otra (buena.m) utilizando el modo alternativo:
x1 = L/2a , x2 = 2c/L
donde L = −b − signo(b)
√
b2 − 4ac.
La sintaxis de la llamada a las funciones ha de ser
>> [x1,x2]=buena(a,b,c);
y similarmente para mala.m.
Discutir la ventaja de buena.m sobre mala.m escribiendo un fichero con un ejemplo ilustrativo.](https://image.slidesharecdn.com/prac1-170524204620/85/Prac1-2-320.jpg)
Prac1
- 1. PR´ACTICAS DE C´ALCULO NUM´ERICO I PR´ACTICA 1 Para la realizaci´on de las pr´acticas vamos a utilizar MATLAB. La raz´on de utilizar este software comercial (y no lenguajes como Fortran o C) es que nos va a permitir realizar de forma integrada la im- plementaci´on de los algoritmos y la visualizaci´on gr´afica de los resultados derivados de ellos. Asimismo, cualquier trabajo de ´ındole matricial se simplifica enormemente trabajando en MATLAB (como vere- mos). Disponemos de unas notas sobre MATLAB que conviene tener siempre a mano. Estas notas no pretenden ser, sin embargo, una gu´ıa exhaustiva sino un resumen de los principales elementos que ir´an asociados a la realizaci´on de las pr´acticas. En cualquier caso, siempre tenemos el comando help en caso de duda. Se recomienda que todos los ejercicios se realicen en ficheros MATLAB (es decir, ficheros de texto con extensi´on .m), en lugar de introducir los comandos en l´ınea. En este primera pr´actica realizaremos una serie de ejercicios sencillos con el objetivo de repasar conceptos b´asicos del an´alisis num´erico y refrescar nuestros conocimientos de programaci´on en el entorno Matlab. 1 Ejercicios: Se sugieren los siguientes ejercicios: 1. Comprobar que el epsilon-m´aquina es 2−52 = 2.2204 10−16, tecleando en la l´ınea de comandos: >> a=1+2^(-53);b=a-1 y comparando con >> a=1+2^(-52);b=a-1 2. Escribir la secuencia de comandos: x=0; while x~=10 x=x+0.1 end en un fichero (con extensi´on .m) y ejecutarlo en MATLAB. Para interrumpir la ejecuci´on, pulsar CTRL+C. ¿Qu´e ocurre si en lugar de incrementarse la variable en 0.1 lo hace en 0.125?. ¿Por qu´e?. 3. Obtener el mayor y el menor n´umero positivo en punto flotante (n´umeros de overflow y underflow). Para obtener el n´umero de overflow escribir un bucle que vaya calculando las sucesivas potencias de 2 y que finalice cuando se produce overflow. El n´umero de underflow se puede obtener calculando las sucesivas potencias negativas de 2 hasta obtener un n´umero indistinguible del cero en punto flotante.
- 2. 4. Escribir dos funciones MATLAB para la resoluci´on de ecuaciones de segundo grado ax2 +bx+c = 0, una de ellas (que llamaremos mala.m) implementando la f´ormula: x = −b ± b2 − 4ac 2a y otra (buena.m) utilizando el modo alternativo: x1 = L/2a , x2 = 2c/L donde L = −b − signo(b) √ b2 − 4ac. La sintaxis de la llamada a las funciones ha de ser >> [x1,x2]=buena(a,b,c); y similarmente para mala.m. Discutir la ventaja de buena.m sobre mala.m escribiendo un fichero con un ejemplo ilustrativo.