#include <stdio.h>
#include <stddef.h>
#include <stdlib.h>
#define MY_END 1

void free_dmatrix(double **m, long nrl, long nrh, long ncl, long nch)
;
void free_matrix(float **m, long nrl, long nrh, long ncl, long nch)
;
void free_dvector(double *v, long nl, long nh)
;
void free_vector(float *v, long nl, long nh)
;
double **dmatrix(long nrl, long nrh, long ncl, long nch)
;
float **matrix(long nrl, long nrh, long ncl, long nch)
;
double *dvector(long nl, long nh)
;
float *vector(long nl, long nh)
;
int *ivector(long nl, long nh)
;
void free_ivector(int *v, long nl, long nh)
;
#define SPINIT(a,i,j,k,l) g=a[i][j];h=a[k][l];a[i][j]=g-s*(h+g*tau);\
	a[k][l]=h+s*(g-h*tau);