KobayashiRui

#include<stdio.h>
#include<Python.h>
#include<string>
#include<iostream>


int main(){
    PyObject *pName,*pModule, *pTmp, *pFunc;
    char *sTmp;
    int data;

mlohry

#include <iostream>
#include <math.h>
#include <Eigen/Dense>
#include <petscts.h>

typedef struct {
  //  Vec         solution; /* global exact solution vector */
  Eigen::MatrixXd soln,rhs;
  const double y0 = 1.0e-3; // epsilon in original ODE
} AppCtx;