read_matrix<-function(data, str){
   n<-length(str);n
   lab1<-paste("Z", str, sep="")
   D<-matrix(99, ncol=n+1, nrow=n+1, 
      dimnames=list(c("X1", lab1), c("X1", lab1)))
   for(i in 1:(n+1))
      D[i,i]<-0

   for (k in 1:nrow(data)){
      Qi<-0
      if (data$i[k]=="X1"){
         for (i in str){
            Qi<-Qi+1 
            z<-paste("Z", i, sep="") 
            if (data$j[k]==z){
               D[1, Qi+1] <- data$km[k]
               D[Qi+1, 1] <- data$km[k]
            }
         }
      }

      Qi<-0
      for (i in str){
         z<-paste("Z", i, sep="") 
         Qi<-Qi+1
         if (data$i[k]==z){
            Qj<-0
            for (j in str){
               z<-paste("Z", j, sep="") 
               Qj<-Qj+1
               if (data$j[k]==z){
                  D[Qi+1, Qj+1] <- data$km[k]
                  D[Qj+1, Qi+1] <- data$km[k]
               }
            }
         }
      }
   }
   D
}