Formula de Viete(1593)

 disp('Cálculo de Pi por la Fórmula de Viete');
n=input('valor de n:');
d=0;
h=1;
m=1;
for i=1:n
    d=sqrt(2+d);
    m=m*d;
    h=h*2;
    p=(2*h)/m;
end;
fprintf('el valor aproximado de pi es %.8f',p);

Ejercicio 3. Cubos en sumas

disp('Cubos en sumas');
a=input('valor de a:');
S=0;
n=a*a;
for i=1:n
    S=S+a;
end
fprintf('el valor del cubo es %.f',S);

Ejercicio3 Cuadrados en sumas

disp('Cuadrados en sumas');
a=input('valor de a:');
S=0;
n=a;
for i=1:n
    S=S+a;
end
fprintf('el valor del cuadrado es %.f',S);

Ejercico 2. Primeros n cubos por la Conjetura de Nicómaco

n=input('valor de n:');
S=0;
for i=1:n
    S=S+n*(n-1)+(2*i-1);
end
fprintf('el valor del cubo es%.2f',S);

Ejercicio 1

n	a	b	Exacto	Trapezoide	Simpson	Error Trapezoide	Error Simpson
6	1.1	2	2.30258509	2.46348358	2.30794875	0.16089849	0.00536366
10	1.1	2	2.30258509	2.36521443	2.30356496	0.06262934	0.00097987
6	2	3	0.69314718	0.69487734	0.69314866	0.00173016	0.00000148
10	2	3	0.69314718	0.69377140	0.69314737	0.00062422	0.00000019

Metodo de simpson

disp('Metodo de Simpson');
a=input('valor de a:');
b=input('valor de b:');
n=input('valor de n:');
Ar=0;
h=(b-a)/n;
     for i=0:n-1
         xi=a+i*h;
         xii=a+(1+i)*h;
         xm=(xi+xii)/2;
         S=1/(1-xi)+4*(1/(1-xm))+1/(1-xii);
         if (S<0)
             S=(-1)*S;
         else
             S=S*1;
         end
            
         Ar=Ar+ h/6*(S);
     end
fprintf('el valor del area es %.8f',Ar);

Metodo del trapecio

disp('Metodo del Trapecio');
a=input('valor de a:');
b=input('valor de b:');
n=input('valor de n:');
Ar=0;
h=(b-a)/n;
     for i=0:n-1
         xi=a+i*h;
         xii=a+(1+i)*h;
         S=1/(1-xi)+1/(1-xii);
         if (S<0)
             S=(-1)*S;
         else
             S=S*1;
         end
            
         Ar=Ar+ h/2*(S);
        
     end
fprintf('el valor del area es %.8f',Ar);