Advanced Search

Question: Gausseidel's method

Question:Gausseidel's method

mvcr15 10 Maple
homework
August 09 2013
1
HI, I have a doubt to how to execute the program with an example. gaussseidel:= proc (A,b,x,p)  local i, j, m, soma, e, g, C,   xk:= Vector(n),  k:=0;  for i from 1 to n do   g[i]:= b[i]/a[i,i];  for j from 1 to n do   if i=j then     C[i,j]:= 0;  else   C[i,j]:= - a[i,j]/a[i,i];  end if;  end do;  end do ;  while  (x - x0) > p do  for i from 1 to n do   soma:= 0;  for j from 1 to n do   soma := soma + C[i,j]*x[j];  end do;   x[i]:= soma + g[i]; end do; k := k + 1; end do;  print (x); end; A:=Matrix([[4, 1, 1], [-2, 5, 1], [3, 1, 6]]); b:= Vector([5,0,-6.5]); x:=Vector([0,0,0]); gaussseidel(A,b,x,3);This is my program but does not it execute the example.
Please Wait...