#include<iostream.h>
#include<iomanip.h>
#include<math.h>

void forward(int n, int *c,double *pi, double **o, double **a, double **state);
void backward(int n, int *c,double *pi, double **b, double **a, double **state);
double forbackward(int n, double **o, double **b);
void optimal_state(int n, double **r, double **o, double **b,double **a, double **state);
void viterbi(int n, int *c, double *pi, double **r, double **o, double **b, double **delta, int **psai, double **a, double **state);
void parameter_estimation(int n, int *c,double *pi,double *sum4, double **r, double **o, double **b, int **psai, double **delta, double ***ksai, double **a, double **state);
void searchmax(int n, double **r);
void searchmax1(int n, int **psai, double **delta, double ***g);
int searchmax2(int n, double **delta);	 
void compute_ksai(int n, int *c, double **o, double **b, double ***ksai, double **a, double **state);


   double start_a[3][3] = {0.9,0.05,0.05,                //a[t