viernes, 30 de septiembre de 2011

mi blog

soluciones de la practica y programas

SEIDEL

http://www.geronet.com.ar/?p=6


Aqui les dejo un link... que resuelve ecuaciones para SEIDEL....

PAGINA MUY INTERESANTES, en las cuales pueden encontrar muchos ejemplos y el codigo en diferentes lenguajes de Factorizacion LU... :D
http://es.wikipedia.org/wiki/Factorizaci%C3%B3n_LU
http://www.mty.itesm.mx/etie/deptos/m/ma95-843/lecturas/l843-34.pdf
http://cicia.uprrp.edu/publicaciones/docentes/Factorizacion%20LU.pdf
http://www.telefonica.net/web2/lasmatematicasdemario/Algebra/Matrices/FactLU.htm
PAGINAS DE CODIGO
http://usuarios.multimania.es/neossoftware/lu.cpp
http://www.lawebdelprogramador.com/news/1051912/Codigo_en_java_para_resolver_matrices_metodo_LU.html
http://mate.uprh.edu/~pnm/notas4031/matlab/apendice.htm

wn lenguaje zpl
http://www.saber.ula.ve/bitstream/123456789/16903/1/rey.pdf
http://www.gennio.com/search/Factorizacion-LU

Mi blogpost

Este es mi blogspot
visiten la pagina
En el sgte enlace encontraran el codigo para comprobar el MAL CONDICIONAMIENTO
pseudocodigo-mal condicionamiento
PINCHA EN EL SGTE LINK PARA INGRESAR AL BLOG DEL AUXILIAR CARLOS
http://ayu156.blogspot.com/


ENLACE DEL AYUDANTE

http:\\ayu156.blogspot.com

Este es el link del Blog del Auxiliar Carlos y cuenta con el codigo de los programas de Jacobi, Seidel y LU
http://ayu156.blogspot.com

Enlace al blog del auxi

http://ayu156.blogspot.com/

Blog del Auxiliar

Link para el blog del aux.

Blog del auxiliar

http://www.ayu156.blogspot.com

BLOG DEL AUX...

link del BLOG del AUX-156 CARLOS

enlase blog del ayudante

http://ayu156.blogspot.com/

jueves, 29 de septiembre de 2011

Codigo en Java Seidel-Jacobi

Este el codigo que desarrollé en Java de Seidel y Jacobi. Los metodos de llenar Matrices y mostrar se los dejo a ustedes
//@utor Jose Espejo Cuba
//Despejamos para que la diagonal principal sea cero
public static void despeje(double M[][],double V[],int m){
double d;
for(int i=0;i d=M[i][i];
for(int j=0;j if(i!=j)
M[i][j]=-M[i][j]/d;
else
M[i][j]=0.0;
}
V[i]=V[i]/d;
}
}

//Metodo de Jacobi
public static void Jacobi(double M[][],double v[],double v1[],double v2[],int m){
for(int i=0;i for(int j=0;j if(i!=j)
v2[i] +=M[i][j]*v1[j];
else
v2[i]+=v[i];
}
}
}
//Metodo de Seidel
public static void Seidel(double M[][],double v[],double v1[],double v2[],int m){
for(int i=0;i for(int j=0;j if(j v2[i] +=M[i][j]*v2[j];
else
v2[i] +=M[i][j]*v1[j];
}
v2[i]+=v[i];
}
}
//Esta es mi condicion de parada
public static double condicion(double v1[],double v2[],int m){
double may1,may2,error;
may1=may2=error=0.0;
for(int i=0;i error=Math.abs(v2[i]-v1[i]);
if(error>may1)
may1=error;
if(v2[i]>may2)
may2=v2[i];
}
return(may1/may2);
}
//Iteramos
public static void iterar(double M[][],double V[],int m){
int cont=0;
double cond=99999;
double x1[]=new double[m];
double x0[]=new double[m];
Arrays.fill(x1, 0.0);
Arrays.fill(x0, 0.0);
while(cond>=0.00000006){
//Seidel(M,V,x0,x1,m);
Jacobi(M,V,x0,x1,m);
cond=condicion(x0,x1,m);
x0=Arrays.copyOf(x1, m);
Arrays.fill(x1, 0.0);
cont++;
}
System.out.println("Numero de iteraciones: "+cont);
for(int i=0;i System.out.println(x0[i]);
}
System.out.println();
}