KM2016 / Parabolic / dif6

 void dif6(double (*funk)(int i,double x[],int n,double t,double xi[],
               double hm[],double hp[],int multiindex[],double tminus,
	       double tplus),
	       double (*funa)(int i,double x[],int n,double t,
               double xi[],double hm[],double hp[],int multiindex[],
               double tminus,double tplus),
	       double (*funb)(int i,double x[],int n,double t,
               double xi[],double hm[],double hp[],int multiindex[],
               double tminus,double tplus),
	       double (*fung)(int i,double x[],int n,double t,
               double xi[],double hm[],double hp[],int multiindex[],
               double tminus,double tplus),
	       double (*funf)(int i,double x[],int n,double t,
               double xi[],double hm[],double hp[],int multiindex[],
               double tminus,double tplus),
	       double (*func)(int nbound,int i,double x[],int n,double t),
	       double (*fund)(int nbound,int i,double x[],int n,double t),
	       double (*fune)(int nbound,int i,double x[],int n,double t),
	       double (*funu0)(int i,double x[],int n,double t,
               double xi[],double hm[],double hp[],int multiindex[],
                              double tminus,double tplus),
	       void (*out)(int iter,int noftau,double t1,double *dt,
		 int n,int mi[],double solution[],double xmi[],int *iterm),
       void (*inform)(int *ind,int ifunk[][6],int ifuna[][6],int ifunb[][6],
       int ifung[][6],int ifunf[][6]),
       int n,int mi[],double xmi[],double t0,int noftau,double sig[],
       int ntime,int nspace[],int nappr,int ind,double tau,int *iterm);

   static void tstep(double (*funk)(int i,double x[],int n,double t,
               double xi[],double hm[],double hp[],int multiindex[],
                              double tminus,double tplus),
	       double (*funa)(int i,double x[],int n,double t,
               double xi[],double hm[],double hp[],int multiindex[],
                              double tminus,double tplus),
	       double (*funb)(int i,double x[],int n,double t,
               double xi[],double hm[],double hp[],int multiindex[],
                              double tminus,double tplus),
	       double (*fung)(int i,double x[],int n,double t,
               double xi[],double hm[],double hp[],int multiindex[],
                              double tminus,double tplus),
	       double (*funf)(int i,double x[],int n,double t,
               double xi[],double hm[],double hp[],int multiindex[],
                              double tminus,double tplus),
	       int ifunk[][6],int ifuna[][6],int ifunb[][6],
	       int ifung[][6],int ifunf[][6],
	       int n,int mi[],double xmi[],double t0,double tau,
	       double sig[],int nappr,double x[],double h[],double z[],
	       double bp[],double bm[],double e[],double f[],
	       double g[],double v[]);


 static void wsp(int nc,int ind,double xmi[],int mi[]);

  static double funq(int i,double x[],int n,double t,
	   double xi[],double hm[],double hp[],int multiindex[],
			      double tminus,double tplus);

  static double funp(int i,double x[],int n,double t,
	   double xi[],double hm[],double hp[],int multiindex[],
			      double tminus,double tplus);

  static void gauss(int ntime,int nspace[],int n);

  static void  wgauss(int n,double ak[],double sk[]);

 static double funfi(double delt);

 static double progco(double (*fun)(int i,double x[],int n,double t,
    double xi[],double hm[],double hp[],int multiindex[],
                              double tminus,double tplus),
 int i,int n,int ifun[][6],double t0,double tau,int j,int k,double x[],
 double xi[],double hm[],double hp[],int multiindex[],
                              double tminus,double tplus);

  static double tksi(double (*fun)(int i,double x[],int n,double t,
       double xi[],double hm[],double hp[],int multiindex[],
                              double tminus,double tplus),
  int i,double t,int n,int ifun[][6],int j,int k,double x[],double xi[],
double hm[],double hp[],int multiindex[],double tminus,double tplus);

static double tjksi(int nvar,double (*fun)(int i,double x[],int n,double t,
                  double xi[],double hm[],double hp[],int multiindex[],
                              double tminus,double tplus),
int i,double t,int n,int ifun[][6],int j,double x[],double xi[],double hm[],
double hp[],int multiindex[],double tminus,double tplus);

static void raz1(int mm,double h[],int ig[],double sig,double tau,
	double q[],double ro[],double rk[],double z[],double bp[],
	double bm[],double g[],double f[],double v[]);

static void ly(int mm,double z[],double h[],double c[],double a[],double b[]);

static void lv(int mm,double z[],double h[],double v[],double a[]);

static void betay(int mm,double bp[],double bm[],double h[],
           double c[],double a[],double b[]);

static  void betav(int mm,double bp[],double bm[],double h[],
double v[],double a[]);

static void progon(int n,double c[],double a[],double b[],double g[]);

//void n_enough_memory(void);