Practica 3 de CNU

EJERCICIO 1:

function [x1,x2] = Ejer1(a,b,c,n)

// Calcula las dos raices de la ecuacion de segundo grado ax2+bx+c=0,

// realizando las operaciones con n digitos de precision.

y = prod_r(b,b,n); // calcula el cuadrado de b

z = prod_r(4,prod_r(a,c,n),n); // calcula 4ac

d = resta_r(y,z,n); // calcula el discriminante

s = redond(sqrt(d),n); // calcula la raiz

x1 = div_r(suma_r(-b,s,n),prod_r(2,a,n),n); // x1 = (-b+s)/(2*a)

x2 = div_r(resta_r(-b,s,n),prod_r(2,a,n),n); // x2 = (-b-s)/(2+a)

EJERCICIO 2:

function tabla = Ejer2(a,b,c)

p = [a b c];

x=roots(p);

tabla=[];

for n=1:16;

[x2,x1] = Ejer1(a,b,c,n);

ec1=(x(1)-x1)/x(1);

ec2=(x(2)-x2)/x(2);

tabla=[tabla;n,ec1,ec2];

end

EJERCICIO 3:

function [x1,x2] = Ejer3(a,b,c,n)

y = prod_r(b,b,n); // calcula el cuadrado de b

z = prod_r(4,prod_r(a,c,n),n); // calcula 4ac

d = resta_r(y,z,n); // calcula el discriminante

s = redond(sqrt(d),n); // calcula la raiz

x1 = div_r(prod_r(-2,c,n),suma_r(b,s,n),n); // x1 = (-2*c)/(b+s)

x2 = div_r(prod_r(-2,c,n),resta_r(b,s,n),n); // x1 = (-2*c)/(b-s)

EJERCICIO 7:

function y=Ejer7(x0,h)

y=(f(x0+h)-f(x0))/h;

EJERCICIO 8:

tabla=[];

for i=1:1:16

h(i)=10^(-i);

der(i)=(f(2+h(i))-f(2))/h(i);

err(i)=abs((der(i)+0.25)/(-0.25));

tabla=[tabla;i,der(i),err(i)];

end

EJERCICIO 9:

for i=1:16

h(i)=10^(-i);

der(i)=Ejer7(2,h(i));

err(i)=abs(der(i)-df(2))/df(2);

end

loglog(h,err)