[file head]:

#include <fcntl.h>
#include <sys/file.h>
#include <tina/all_tina.h>


/*_______________________ runerror __________________________*/
void runerror(error_text)
char error_text[];
{
/*
void exit();
*/

fprintf(stderr, "Run-time fatal error ... \n");
fprintf(stderr, ">s\n", error_text);
fprintf(stderr, "... now exiting to system.\n");
exit(1);
}

/*_____________________ int_vector _________________________*/
int *int_vector(int n)
{
int *v;

v = (int *)malloc((unsigned) (n*sizeof(int)));
if(!v) runerror("alloaction failure in int_vector()");
return v;
}

/*_________________________ free_int_vector ____________________*/
void free_int_vector(int *v)
{
free((int*) v);
}
/*_____________________ float_vector _________________________*/
float *fl_vector(int n)
{
float *v;

v = (float *)malloc((unsigned) (n*sizeof(float)));
if(!v) runerror("alloaction failure in fl_vector()");
... ...

[file tail]:
... ...
-ascending numerical order-------*/
void sorting(data,n)
float data[]; int n;
{
int i, j, min;
float dmin;

for(i=0; i<n-1; i++)
{
min = i;
dmin = data[i];
for(j=i+1; j<n; j++)
if(data[j]<dmin)
{
min = j;
dmin = data[j];
}
data[min] = data[i];
data[i] = dmin;

}
}

int get_median(data,n)
float data[]; int n;
{
int i;
float sorted[n];
int median = 0;

for(i = 0;i < n;i ++)
sorted[i] = data[i];
sorting(sorted,n);
if(fmod(n,2) == 1)
median = sorted[((n-1)/2)];
else
median = ((sorted[n/2] + sorted[n/2-1])/2);
//printf(">d\n",median);
return median;

}

int get_mode(data,n)
float data[]; int n;
{
int i,t,max;
int gray[256];
int freqgray;
for(i = 0; i < 256; i ++)
gray[i] = 0;
for(i = 0; i < n; i ++)
{
t = (int)data[i];
gray[t] ++;
}
max = gray[0];
freqgray = 0.0;
for(i = 1; i < 256; i ++)
if(gray[i] > max)
{
max = gray[i];
freqgray = i;
}
//printf(">d,>d\n",max,freqgray);
return freqgray;
}