function [ sols ] = BackSub( A,b ) % [n,m]=size(b); b(n,:)=b(n,:)/A(n,n); for i=(n-1):-1:1 b(i,:)=(b(i,:)-A(i,(i+1):n)*b((i+1):n,:))/A(i,i); end sols=b; end