/*
  Markov連鎖の母期待値の時間発展の計算
  Time-stamp: "2021-05-31 Mon 17:47 JST hig"
*/
/* これじゃなく markovpde03 やその完成ソースを参考にしてね */

#define _CRT_SECURE_NO_WARNINGS /* Visual C++ に必要なおまじない */
#include <stdio.h>

#define NS 41

int multiply_trans(double pn[], double p[], int d);
int measure_dist(double p[], int t, int m);


int main(){
  int t,tmax;
  int x;
  double p[NS];
  double pn[NS];
  int m=NS;
  int x0i=20;
  
  scanf("%d", &tmax);
  printf("#T=%d\n", tmax);
  printf("t,");
  for(x=0;x<m;x++){
    printf("%d,",x);
  }
  printf("\n");
  
  /* ここに初期条件書いてね */
  t=0;
  for(x=0;x<m;x++){
    p[x]=0.0;
  }
 
  measure_dist(p,t,m);
  
  for(t=1;t<=tmax;t++){
    
    multiply_trans(pn,p,m);

    for(x=0;x<m;x++){
      p[x]=pn[x];
    }

    measure_dist(p,t,m);
  
  }

  return 0;
}


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

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

}


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