1 // Make sure file is only included once per object
 2 #ifndef FUNCTIONS_H_
 3 #define FUNCTIONS_H_
 4 
 5 #include <vector>
 6 #include <complex>
 7 #include "typedefs.h"
 8 
 9 
10 //// ODEs with real+complex parts.
11 //// Return the derivatives at the point x,vec(y)
12 
13 std::vector<std::complex<Floattype> >
14         func1(const std::complex<Floattype> z,
15               const std::vector<std::complex<Floattype> > &y)
16 {
17   std::vector<std::complex<Floattype> > dydz(2);
18   dydz[0].real() = y[1].real();
19   dydz[0].imag() = y[1].imag();
20   dydz[1].real() = -y[0].real();
21   dydz[1].imag() = 0.5f*y[0].imag();
22   return dydz;
23 }
24 
25 #endif