jacobi - end disp('MATRIZ T DE LAS ITERACIONES') disp(a)...

Info iconThis preview shows pages 1–4. Sign up to view the full content.

View Full Document Right Arrow Icon
Sheet1 Page 1 % PROGRAMA: METODO DE JACOBI % Descripcion:El programa encuentra las raices de un conjunto de ecuaciones lineales de la forma % A x X = B, utilizando el metodo iterativo de Jacobi. % curso: Calculo Numerico % programacion: Rony Yupanqui G. % fecha: Setiembre 2006 %------------VARIABLES---------% % numero de ecuaciones lineales n=3 % matrices utilizadas a=zeros(n,n) b=zeros(n,1) x0=zeros(n,1) x1=zeros(n,1) % tolerancia tol=0.00001 % valores de la matriz A a=[3 1 0 -1 3 1 % segundos miembros de las ecuaciones b=[1 -1 % aproximacion inicial x0=[0 0 %----------ALGORITMO DEL PROGRAMA-----% disp(' PROGRAMA: METODO DE JACOBI') for i=1:n, val=-1*a(i,i) for j=1:n, a(i,j)=a(i,j)/val end a(i,i)=0 b(i,1)=b(i,1)/val*-1
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 2
Background image of page 3

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 4
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: end disp('MATRIZ T DE LAS ITERACIONES') disp(a) disp('VECTOR C DE LAS ITERACIONES') disp(b) error=1 num=1 disp(' ITERACIONES: ') disp(' N X(1) X(2) X(3) ERROR') while error>tol, for i=1:n, val=0 for j=1:n, val=val+a(i,j)*x0(j,1) end x1(i,1)=val+b(i,1) end % tipo de error error=0 Sheet1 Page 2 for i=1:n, error=error+(x1(i,1)-x0(i,1))*(x1(i,1)-x0(i,1)) end error=sqrt(error) % mostrar el resultado de la iteracion con los elementos del vector x linea=sprintf(' %d %5.8f %5.8f %5.8f %5.8f',num,x0(1,1),x0(2,1),x0(3,1),error) disp(linea) x0=x1 num=num+1 end disp('SOLUCIONES X DE LAS ECUACIONES LINEALES') disp(x0) Sheet1 Page 3 0 2 3] 0] 0] Sheet1 Page 4...
View Full Document

Page1 / 4

jacobi - end disp('MATRIZ T DE LAS ITERACIONES') disp(a)...

This preview shows document pages 1 - 4. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online