u1=teta''+c1 teta'+ c5 alpha'+c2 alpha+e3 v+c3 vis+c12 zak==0;
u2=-teta'+alpha'+c4 alpha+e2 v+c10 teta+c9 vis+c13 zak==0;
u3=v'+e1 v+c8 alpha+c7 teta+c14 zak==0;
u4=h'-c11 v-c6 teta+c6 alpha==0;
MatrixForm[u={u1,u2,u3,u4}];

u=ReplaceAll[u,{teta''->wz',teta'->wz}];

u5=teta'-wz==0;
u=Flatten[{u,u5}];
MatrixForm[u];

qq=Flatten[Solve[u,{wz',alpha',v',h',teta'},Sort->False]];
MatrixForm[qq];

qq=Collect[qq,{wz,alpha,v,h,teta,vis,zak}];
MatrixForm[qq]