Proyecto Metodos

INSTITUTO POLITÉCNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERÍA MECANIA Y ELECTRICA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA
ACADEMIA DE COMPUTACIÓN
MÉTODOS NUMÉRICOS

TAREA 1: PROYECTO FINAL

Nombre:


Boleta:

Grupo: 3EM5

Carrera: INGENIERIA ELECTRICA

Profesor:ENRIQUE MARTINEZ ROLDAN

Fecha: 2 DICIEMBRE 2011

PROYECTO

OBJETIVO:

Hacer un menúimplementando los métodos vistos en clase.

LOCALIZACION DE RAICES
BISECCION

#include <stdio.h>
#include <conio.h>
#include <math.h>

int main()
{
float xant, m, m2, crit, raiz, tol, xi, xf, yxf, xn, yxn, cont, yxi, xnact;
xn=0;

printf("\n\n\t\t\t =Metodo de la Biseccion=");
cont=0;
printf("\n\n\n Inserta el valor de X inicial: \n\nXi= ");
scanf("%f",&xi);
printf("\n\n Inserta el valor de X final: \n\n Xf= ");
scanf("%f",&xf);
printf("\n\n Inserta el rango de tolerancia: \n\n Tol= ");
scanf("%f",&tol);
do
{
cont=cont+1;
xant=xn;
xn=(xi+xf)/2;
yxi=(xi*xi*xi)+(4*(xi*xi))-10;
yxf=(xf*xf*xf)+(4*(xf*xf))-10;yxn=(xn*xn*xn)+(4*(xn*xn))-10;
m=yxi*yxn;
m2=yxf*yxn;
if(m<0)
{
xf=xn;
}
else
{
if(m2<0)
{
xi=xn;
xnact=xn;
crit=fabs((xnact-xant)/xnact);
}
}
}
while(fabs(yxn)>=tol||crit>=tol);
raiz=xn;printf("\n\n La raiz es: %f",raiz);
printf("\n\n\n El numero de iteraciones realizadas es: %.2f\n\n",cont);

return 0;

}

SECANTE

#include <math.h>
#include <stdio.h>
#include <stdlib.h>

double xi0, xi1, error;
int max_iter;
void lee_datos()
{
printf("===========================INSTITUTO POLITECNICO NACIONAL=======================\n");printf("\n\t\t\t\tMETODO DE LA SECANTE");
printf("\n\n\tecuacion:x^2+2x");
printf("\n\nDame el valor de xi1:\n");
scanf("%lf",&xi1);
printf("\n\nDame el valor de xi0:\n");
scanf("%lf",&xi0);
printf("\n\nCual es el error permitido\n");
scanf("%lf",&error);
printf("\n\nCual es el maximo de ciclos\n");
scanf("%lf",&max_iter);
}
double funcion(double x)
{
return(pow(x,2)+(2*x));
}

int main()
{
double error_aprox, prod,xr,xn;
int ciclos=0;
lee_datos();
xr=((funcion(xi1)*(xi0-xi1))/(funcion(xi0)-funcion(xi1)));
error_aprox=1;
printf("\n ciclos xi0 xi1 xn error \n");
printf("\n %3d %10.4f %10.4f %10.4f %f\n", ciclos, xi0, xi1, xn, error );
while (ciclos<=max_iter&&error<error_aprox){
prod=funcion(xi1)*funcion(xr);
if (prod==0)
printf("\nla raiz es: %10.4f\n",xr);
xi0=xi1;
xi1=xr;
xn=((funcion(xi1)*(xi0-xi1))/(funcion(xi0)-funcion(xi1)));
ciclos++;
error_aprox=fabs((xn-xr)/(xn*100));
printf("\n %3d %11.4lf %17.4lf %17.4lf %17.4lf\n", ciclos, xi0, xi1, xn, error_aprox );
xr=xn;
}
if(ciclos<max_iter)
{
printf("\nLa raiz de la funcion es: %lf\n", xn);
printf("\n Se encontro en %d iteraciones\n\n", ciclos);
}
else
{
printf("\n No se encontro la raiz en %d iteraciones\n", ciclos);
}

system("pause");

return 0;
}

SISTEMA DE ECUACIONES LINEALES
ELIMINACION DE GAUSS
#include<stdio.h>
#include<stdlib.h>#include<conio.h>

int main()
{
int i,j,n,k;
float A[20][20],B[20],piv,cte;
printf("Eliminacion de Gauss\n\n");
printf("Numero de Ecuaciones: ");
scanf("%d",&n);
printf("\nDigite los Coeficientes:\n\n");
for(i=1;i<=n;i++)
{
for(j=1;j<=n;j++)
{
printf("A(%d,%d) = ? ",i,j);
scanf("%f",&A[i][j]);
}...