You are not logged in.
There are other reasons for divergence. Poor initial conditions is one reason. The condition number for that matrix if I remember was around 500. That is not too bad.
In mathematics, you don't understand things. You just get used to them.
I agree with you regarding the satisfaction and importance of actually computing some numbers. I can't tell you how often I see time and money wasted because someone didn't bother to run the numbers.
Offline
There are other reasons for divergence. Poor initial conditions is one reason. The condition number for that matrix if I remember was around 500. That is not too bad.
The condition number for n=500 is for example 1.2550e+005 ..
So,is it right that the methods do not converge for this matrix ??
Because,for example,both of the methods do converge for an other tridiagonal matrix,that has the number 4 at the main diagonal,1 at the first diagonal below this and also 1 at the first diagonal above this.
Which is the difference??
Offline
I believe that there is a simple answer to that.
In mathematics, you don't understand things. You just get used to them.
I agree with you regarding the satisfaction and importance of actually computing some numbers. I can't tell you how often I see time and money wasted because someone didn't bother to run the numbers.
Offline
I believe that there is a simple answer to that.
Could you tell me,which it is??:
Offline
You have a 4 on the diagonal and 1 on the other 2 diagonals? That matrix according to some books is diagonally dominant. This makes convergence easy in the Gauss Seidel method.
But we are straying from the point of this thread. Did you try post #18's matrices? Did they converge for you? What did you notice about the new A?
In mathematics, you don't understand things. You just get used to them.
I agree with you regarding the satisfaction and importance of actually computing some numbers. I can't tell you how often I see time and money wasted because someone didn't bother to run the numbers.
Offline
You have a 4 on the diagonal and 1 on the other 2 diagonals? That matrix according to some books is diagonally dominant. This makes convergence easy in the Gauss Seidel method.
I understand...
But we are straying from the point of this thread. Did you try post #18's matrices? Did they converge for you? What did you notice about the new A?
No,the matrix of post #18 does not converge,neither using the jacobi,nor the gaussseidel method..
Offline
Something is wrong with your program. That converges for me. Are you setting initial conditions to [0,0,0,0]?
Check this page to see how to increase Matlabs precision:
http://www.mathworks.com/help/symbolic/ … igits.html
See you a bit later, I need some rest.
In mathematics, you don't understand things. You just get used to them.
I agree with you regarding the satisfaction and importance of actually computing some numbers. I can't tell you how often I see time and money wasted because someone didn't bother to run the numbers.
Offline
Something is wrong with your program. That converges for me. Are you setting initial conditions to [0,0,0,0]?
Check this page to see how to increase Matlabs precision:
http://www.mathworks.com/help/symbolic/ … igits.html
See you a bit later, I need some rest.
I used the command digits bit I didn't get any result!!! why??
Offline
The correct comand is DIGITS:= something.
In mathematics, you don't understand things. You just get used to them.
I agree with you regarding the satisfaction and importance of actually computing some numbers. I can't tell you how often I see time and money wasted because someone didn't bother to run the numbers.
Offline
I get this error message:Error: The expression to the left of the equals sign is not a valid target for an
assignment.
I would like to ask you also something else..Is this right that the spectral radius of the iteration matrix of the 250x250 Hilbert Matrix,when we apply the GaussSeidel method is 1,and when we apply the Jacobi method 217.3320?
Last edited by evinda (20131202 08:26:01)
Offline
http://www.mathworks.com/help/symbolic/ igits.html
The examples on that page do not work?
In mathematics, you don't understand things. You just get used to them.
I agree with you regarding the satisfaction and importance of actually computing some numbers. I can't tell you how often I see time and money wasted because someone didn't bother to run the numbers.
Offline
http://www.mathworks.com/help/symbolic/ igits.html
The examples on that page do not work?
They work...But I can't apply them at a matrix
Offline
Why not? What are you entering for the matrix?
In mathematics, you don't understand things. You just get used to them.
I agree with you regarding the satisfaction and importance of actually computing some numbers. I can't tell you how often I see time and money wasted because someone didn't bother to run the numbers.
Offline
Why not? What are you entering for the matrix?
I wrote this:DIGITS:=hilb(15) and I get Error: The expression to the left of the equals sign is not a valid target for an
assignment.
Offline
DIGITS is a global variable. It sets the precision for everything. Try DIGITS:= 25
In mathematics, you don't understand things. You just get used to them.
I agree with you regarding the satisfaction and importance of actually computing some numbers. I can't tell you how often I see time and money wasted because someone didn't bother to run the numbers.
Offline
I get this message:??? DIGITS:=25

Error: The expression to the left of the equals sign is not a valid target for an
assignment.
But...anyway...I will let it like that...I have a last question..Could you tell me which is the iteration matrix ,for example of the hilbert matrix with a specific dimension,or for another matrix using the gaussseidel method,and which it is using the jacobi method ,so I can check if my result is right???
Offline
Hi;
Did you get convergence for the matrix in post #18?
In mathematics, you don't understand things. You just get used to them.
I agree with you regarding the satisfaction and importance of actually computing some numbers. I can't tell you how often I see time and money wasted because someone didn't bother to run the numbers.
Offline