KM2016 / fem5 / fem5_0

#include <ctype.h>
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <string.h>

static double funD( double xi, double b [], double alpha [], double r [], double d1 [], double t);
static double funC( double xi, double b [],double alpha [], double r [], double d1 [], double t);
static double funA( double xi, double b [],double alpha [], double t);
static double fun_q(int si, double x [],double b [], double alpha [],double r [],double tau);
static double fun_a(int fb, double t, double alpha [], double r []);

void K_ij(int N, double K [], double x [], double b [], double t, 
	  double alpha [], double r [], double d1 [],
	  double (*funD)(double xi, double b [], double alpha [], double r [], double d1 [], double t), 
	  double (*funC)(double xi, double b [],double alpha [], double r [], double d1 [], double t),
	  double (*funA)(double xi, double b [],double alpha [], double t), 
	  double (*fun_a)(int fb, double t, double alpha [], double r []) );
	  
void q_i(int k, double q[],int knots_N, double x[], double b[],  
	 int si, double alpha[],double r[],double tau,
	 double (*fun_q)(int, double *,double *,
		    double *,double *,double) ); 
		    
 double domik(double x, int j, int mp, double *xx);
 
 double intgD(int mp, int i, double *x, double *b,
 double alpha [], double r [], double d1 [], double t,	     
 double (*fun)( double xi, double b [], double alpha [], double r [], double d1 [], double t));
 
 double intgC(int mp, int i, int mpi, double *x, double *b, 
 double alpha [], double r [], double d1 [], double t,	     
 double (*fun)( double xi, double b [], double alpha [], double r [], double d1 [], double t),
 double (*domik)(double, int, int, double *));
 
 double intgA(int, int ,int, int , int , double *, double *, double *, double,
	      double (*)( double , double *,double *, double), 
	      double (*)(double, int, int, double *));
 int i_max(int N, int p, double V[]);
 int lufact(int, double [], double [], double []);
 void SaveToFile(int, double*, const char* ); 
 void AddToFile(int, double*, const char* );
 //int output(int test, int N, double* x, double *b, double *solution, const char* filename);