function dxdt=igminmod(t,x,x0);
Sg=0.026;Si=5e-4;p1=Sg;p2=0.025;p3=Si*p2;gb=90;ib=8;
data=[0 92 11;2 350 26;4 287 130;6 251 85;8 240 51;10 216 49;
    12 211 45; 14 205 41;16 196 35; 19 192 30;22 172 30; 27 163 27;
    32 142 30; 42 124 22; 52 105 15; 62 92 15; 72 84 11; 82 77 10;
    92 82 8;102 81 11; 122 82 7;142 82 8;162 85 8;182 90 7];
i=ig_iconc(t,data);
dxdt=[-(p1+x(2))*x(1)+p1*gb;-p2*x(2)+p3*(i-ib)];
