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

This note was uploaded on 09/19/2009 for the course MATH numerical taught by Professor Ford during the Spring '09 term at FAU.

Ask a homework question - tutors are online