#include<stdio.h> #include<stdlib.h> #include<math.h> double *ifun(double *r, int N, int T, double Pr, double tau, double kappa, double *ih, double *I); //const unsigned long seed