Errors in bit***.h files fixed

3 files changed, 44 insertions, 0 deletions

diff --git a/src/c/differential_calculus/includes/ode.h b/src/c/differential_calculus/includes/ode.h
new file mode 100644
index 00000000..46220a84
--- /dev/null
+++ b/src/c/differential_calculus/includes/ode.h
@@ -0,0 +1,9 @@
+// FOSSEE IIT Bombay
+#ifdef __INT_ODE_H__
+#define __INT_ODE_H__
+
+double ODE(double initial_value, double start_time, double end_time, \
+			char *ode_function, double nequs, double eps_abs, double eps_rel \
+			double step_size)
+
+#endif /*__INT_ODE_H__*/
\ No newline at end of file
diff --git a/src/c/differential_calculus/interfaces/int_ode.h b/src/c/differential_calculus/interfaces/int_ode.h
new file mode 100644
index 00000000..5000ed5e
--- /dev/null
+++ b/src/c/differential_calculus/interfaces/int_ode.h
@@ -0,0 +1,9 @@
+// FOSSEE IIT Bombay
+
+#ifdef __INT_ODE_H__
+#define __INT_ODE_H__
+
+#define d0d0d0g2ODEd0(in1, in2, in3, in4) ODE(in1, in2, in3, in4, 1, 1.0e-2,\
+											 1.0e-2, 1.0e-6)
+
+#endif /*__INT_ODE_H__*/
\ No newline at end of file
diff --git a/src/c/differential_calculus/ode/ode.c b/src/c/differential_calculus/ode/ode.c
new file mode 100644
index 00000000..b815e367
--- /dev/null
+++ b/src/c/differential_calculus/ode/ode.c
@@ -0,0 +1,26 @@
+// FOSSEE IIT Bombay
+
+#include "ODE.h"
+#include "types.h"
+
+
+double ODE(double initial_value, double start_time, double end_time, \
+			char *ode_function, double nequs, double eps_abs, double eps_rel \
+			double step_size)
+{
+	double out = 0;
+	int status;
+	//Setup ODE related parameters
+	gsl_odeiv2_system sys = {ode_function, NULL, 2, NULL};
+
+	gsl_odeiv2_step    *s = gsl_odeiv2_step_alloc (gsl_odeiv2_step_rkf45, nequs);
+  	gsl_odeiv2_control *c = gsl_odeiv2_control_y_new (eps_abs, eps_rel);
+  	gsl_odeiv2_evolve  *e = gsl_odeiv2_evolve_alloc (nequs);
+  
+  	status = gsl_odeiv2_evolve_apply_fixed_step (e, c, s, &sys, &t, stepsize, x_dot);
+	
+	if (status == GSL_SUCCESS) {
+	return out;
+	}
+	else return 0;
+}