/*
このプログラムは, multiply_trans でかけられる行列 Mの
転置
 すなわち転置推移確率行列の転置, 転置確率行列の転置を出力する.
  Time-stamp: "2020-06-10 Wed 18:24 JST hig"
*/
#define _CRT_SECURE_NO_WARNINGS /* Visual C++ に必要なおまじない */
#include <stdio.h>

#define NS 10

int multiply_trans(double pn[], double p[],int m);
int print_dist(double p[],  int m);

int main(){
  int x;
  double p[NS];
  double pn[NS];
  int m=NS;
  int n;
  
  for(n=0;n<m;n++){
    for(x=0;x<m;x++){
      p[x]=0.0;
    }
    p[n]=1.0;

    /*    print_dist(p,m);*/
    multiply_trans(pn,p,m);
    print_dist(pn,m);

  }
  
  return 0;
}


/* int multiply_trans(double *pn, double *p) を書いてね */ 
/* p に M をかけて pn=M p を書き込む */

/* これは吸収壁の場合. 周期的の場合は特別扱いは2行で済む */
int multiply_trans(double pn[], double p[], int m){
  int y;
  
  for(y=0;y<m-1;y++){
    pn[y]=p[y+1];
  }
  pn[m-1]=p[0];

  return 0;
}

/* tとpを1行に出力*/
int print_dist(double p[], int m){
  int x;
  for(x=0;x<m;x++){
    printf("%.1f, ",p[x]);
  }
  printf("\n");
  return 0;
}