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);