# Wczytanie z pliku macierz.dat macierzy symetrycznej 
# zapisanej w oszczednym formacie trójkątnym
# i zapisanie jej do pliku macierz.out
A=matrix(0,20,20)
x=scan("macierz.dat",skip=2) # wczytanie danych w postaci wektora
for(i in 1:20){
  idx = i*(i-1)*0.5+1
  A[i,1:i] = x[idx:(idx+i-1)]
}
A=A+t(A)
A=A-diag(0.5*diag(A)) 

# Aby print nie wypisała numerów kolumn i wierszy
dimnames(A)=list(rep("",nrow(A)), rep("",ncol(A)))
sink("macierz.out")
print(A)
sink()

# Przy okazji: kasowanie nazwy kolumn macierzy X
# dimnames(X)=NULL
# Zapisywać macierz możemy też za pomocą
# write(t(A),"macierz2.out",ncol=ncol(A))