% see hexch1.m for comments
function [sys,x0] = hexch52(t,x,u,flag)
Indexch=8;
N = 10;
if abs(flag) == 1
  sys = mathisen(t,x,u,N,Indexch);
elseif abs(flag) == 3
  sys(1,1) = x(N);       
  sys(2,1) = x(3*N);    
elseif flag == 0
  x0 = steadyvar(Indexch,N);
  sys = [3*N,0,2,4,0,0]; 
else
  sys = [];  
end
end