long mdfunc (double b,
	     signed short *x,
	     signed short *y,
	     double *a,
	     int n,
	     signed short *tmp);

void l1fit (signed short *x,
	    signed short *y,
	    int n,
	    double *a,
	    double *b);
