How can solve differential equation in Matrix form( dI/dt=c.*(M*I)+c.*(L*V) )

vas=1;vbs=2;vcs=3;var=1.5;vbr=2.5;vcr=3.5; Rs=1.06;Rr=0.8;Lm=0.0664;Ls=0.2065;Lr=0.0810;p=2;Jg=23; S=1-(Lm^2)/(Ls*Lr);Os=1;Or=0;Om=Os-Or;L=Lm/Ls;tm=pi/2; vabcs=[vas;vbs;vcs]; Tm=[cos(tm) -sin(tm); sin(tm) cos(tm)]; Tc=2/3.*[1 -1/2 -1/2;0 sqrt(3)/2 -sqrt(3)/2 ]; vabs=Tc*vabcs;vur=u*vabcr;vDQr=Tc*vur;vabr=Tm*vDQr; syms ias(t) ibs(t) iar(t) ibr(t) I=[ias;ibs;iar;ibr]; c=1/(S*Ls*Lr); c11=-Rs*Lr;c12=Om*Lm^2;c13=Rr*Lm;c14=Om*Lm*Lr; c21=-c12;c22=c11;c23=-c14;c24=c13; c31=Rs*Lm;c32=-Om*Ls*Lm;c33=-Rr*Ls;c34=-Om*Lr*Ls; c41=-c32;c42=c31;c43=-c34;c44=c33; M=[c11 c12 c13 c14 ; c21 c22 c23 c24; c31 c32 c33 c34; c41 c42 c43 c44]; L=[Lr 0 -Lm 0; 0 Lr 0 -Lm ; -Lm 0 Ls 0; 0 -Lm 0 Ls ]; vab=[vabs;vabr]; odes=diff(I)==c.*(M*I)+ c.*(L*vab) [iasSol(t),ibsSol(t),iarSol(t),ibrSol(t)] = dsolve(odes); iasSol(t)= simplify(iasSol(t)) ibsSol(t)= simplify(ibsSol(t)) iarSol(t)= simplify(iarSol(t)) ibrSol(t)= simplify(ibrSol(t)) I can not find any results. And thank you in advance. email: [e-mail]