S = load('matrix.dat');

i = S(:,1) + 1;
j = S(:,2) + 1;
aij = S(:,3);

A = sparse(i,j,aij);