GaussNaive - nb=n+1; Aug=[A b]; %forward elimination for...

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

View Full Document Right Arrow Icon
function x=GaussNaive(A,b) %Problem 9.6 check %GaussNaive: naive Gauss elimination % x=GaussNaive(A,b): Gauss elimination without pivoting. %input: % A=coefficient matrix % b=right hand side vector %output: % x=solution vector [m,n]=size(A); if m~=n, error('matrix A must be squared'); end
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: nb=n+1; Aug=[A b]; %forward elimination for k=1:n-1 for i=k+1:n factor=Aug(i,k)/Aug(k,k); Aug(i,k:nb)=Aug(i,k:nb)-factor*Aug(k,k:nb); end end %back substitution x=zeros(n,1); x(n)=Aug(n,nb)/Aug(n,n); for i=n-1:-1:1 x(i)=(Aug(i,nb)-Aug(i,i+1:n)*x(i+1:n))/Aug(i,i); end...
View Full Document

This note was uploaded on 09/27/2011 for the course EGM 3344 taught by Professor Raphaelhaftka during the Spring '09 term at University of Florida.

Ask a homework question - tutors are online