%%--------------------------------------------------------------------------%%
%% gausspp.m - drive of Gaussian Elimination wit Partial Pivoting           %%
%%--------------------------------------------------------------------------%%
fin=fopen('gaussmat.dat','r');
n=fscanf(fin,'%d',1);
A=fscanf(fin,'%f',[n n]);
A=A';
b=fscanf(fin,'%f',n);
% Y=A\b
[A b]
X=gausspivot(A,b,n)

