MATLAB and C Programs - Digital Signal Processing: A Practitioner's Approach - page 173

{
rnd = rand();
udf = rnd;
udf = (udf/intmax)-0.5;
sum = sum + udf;
}
return(sum);
}

void mat_p( float *some_mat, int max_row, int max_col)
{
int col,row,index;
float val;
for (row =0; row < max_row; row++)

{
for (col =0; col < max_col; col++)

{
index = max_col*row + col;
val = *(some_mat + index);
printf (" %2.3e ",val);
}
printf("\n");
}
printf("\n");
}

float m_rem(float , float *);

float m_rem(float yk, float *pk)
{

static float xk=0, xk1=0.1, xk2;
static float sk=0, sk1=0, sk2=0, sk3, xk_corkt=0;

static float Hk[3] [3] = {{1e-9,0,0},{0, 1e-9,0},{0,0,1e-9}};
static float Hk_new[3][3],delH_num[3][3],delH[3][3],delH_den;

static float ek=0,Jk=0, Jk_hat=0;

static float gk[3] = {0,0,0};