sig=0.1;cout0=Table[If[i≤4,1,1/16],{i,1,6}];(*Forii=1,ii≤iter,ii++,*)ii=1;nder[n_,m_,f_]:=(Derivative@@UnitVector[m,n])[f]cout0=Table[If[i≤4,1,1/16],{i,1,6}];ccvar=Array[cc,6];ccvar=ccvar/.cc[i_]cc[i][t];For[n=1;eqnlist={};,n≤6,n++,eqnlist=Flatten[Append[eqnlist,{cc[n]'[t]signder[n,12,ff][Flatten[{ccvar,ccvar}]],cc[n][0]==cout0[[n]]}]]]ff[{c1_Real,c2_Real,c3_Real,c4_Real,alpha11_Real,alpha21_Real,alpha31_Real,alpha41_Real,alpha12_Real,alpha22_Real,alpha32_Real,alpha42_Real,alpha13_Real}]:=1.solDEsuper=NDSolve[eqnlist,ccvar,{t,0,teq}]
NDSolve::ndnum:Encountered non-numerical value for a derivative at t == 0..
NDSolve[t]0.1[{cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t],cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t]}],cc[1][0]1,[t]0.1[{cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t],cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t]}],cc[2][0]1,[t]0.1[{cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t],cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t]}],cc[3][0]1,[t]0.1[{cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t],cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t]}],cc[4][0]1,[t]0.1[{cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t],cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t]}],cc[5][0],[t]0.1[{cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t],cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t]}],cc[6][0],{cc[1][t],cc[2][t],cc[3][t],cc[4][t],cc[5][t],cc[6][t]},{t,0,0.1}
′
cc[1]
(1,0,0,0,0,0,0,0,0,0,0,0)
ff
′
cc[2]
(0,1,0,0,0,0,0,0,0,0,0,0)
ff
′
cc[3]
(0,0,1,0,0,0,0,0,0,0,0,0)
ff
′
cc[4]
(0,0,0,1,0,0,0,0,0,0,0,0)
ff
′
cc[5]
(0,0,0,0,1,0,0,0,0,0,0,0)
ff
1
16
′
cc[6]
(0,0,0,0,0,1,0,0,0,0,0,0)
ff
1
16