spectralradius

spectralradius - radius = norm(d_evalue(1,1)); for i = 2:m

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
function [convergence,radius,norm2] = spectralradius(A) % % Given matrix, return spectral radius, convergence of the spectral % raidus and the 2nd norm of the original matrix % [evector,d_evalue]=eig(A); [m,n] = size(d_evalue);
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: radius = norm(d_evalue(1,1)); for i = 2:m if(norm(d_evalue(i,i)) > radius) radius = norm(d_evalue(i,i)); end; end; convergence = 'not convergent'; if(radius<1) convergence = 'convergent'; end; norm2 = norm(A);...
View Full Document

Ask a homework question - tutors are online