function dy=fangcheng(t,y)
dy=zeros(11,1);
dy(1)=(203.227*y(1)+200.227*y(3)+200*(y(5)+y(6))+0.6837*dy(3)+0.667*(dy(5)+dy(6))-y(8))/(-0.7437);
dy(2)=(200*y(1)+3.303*y(2)+200*y(3)-0.303*y(4)+200.227*y(5)+206.227*y(6)+0.667*dy(1)+0.667*dy(3)-0.06*dy(4)+0.6837*dy(5)+0.8037*dy(6)+y(11)-y(10))/(-0.0793);
dy(3)=(y(11)+200*y(1)+3*y(2)+200*y(3)+200.227*y(5)+206.53*y(6)+0.303*y(7)+0.667*dy(1)+0.06*dy(2)+0.6837*dy(5)+0.823*dy(6)+0.0193*dy(7)-y(10))/(-0.667);
dy(4)=(y(9)-y(11)+0.227*y(1)+4.53*y(3)-0.303*y(4)-4.227*y(5)-0.227*y(6)+0.0167*dy(1)+0.256*dy(3)-0.2367*dy(5)-0.0167*dy(6)-y(8))/0.0193;
dy(5)=(y(9)-y(11)+0.227*y(1)+4.227*y(3)-4.53*y(5)-0.227*y(6)+0.303*y(7)+0.0167*dy(1)+0.2367*dy(3)-0.0167*dy(6)+0.0193*dy(7)-y(8))/0.256;
dy(7)=(y(10)-0.303*y(2)-0.303*y(3)+3.606*y(4)-3*y(7)-0.0193*dy(2)-0.0193*dy(3)+0.1286*dy(4)-y(9))/0.09;
dy(6)=(y(10)+3*y(4)+0.303*y(5)-0.303*y(6)-3.606*y(7)+0.09*dy(4)+0.0193*dy(5)-0.1286*dy(7)-y(9))/0.0193;
dy(6)=(y(10)-0.303*y(3)+3.303*y(4)-0.303*y(6)-3.303*y(7)-0.0193*dy(3)+0.1093*dy(4)-0.1093*dy(7)-y(9))/0.0193;
dy(7)=(y(10)-0.303*y(2)+3.303*y(4)+0.303*y(5)-3.303*y(7)-0.0193*dy(2)+0.1093*dy(4)+0.0193*dy(5)-y(9))/0.1093;
dy(8)=(1/(3.6e-4))/(-y(1)-y(3));
dy(9)=(1/(1.8e-4))/(y(3)-y(4));
dy(9)=(1/(1.8e-4))/(-y(5)+y(7));
dy(10)=(1/(1.8e-4))/(-y(2)+y(4));
dy(10)=(1/(1.8e-4))/(-y(6)-y(7));
dy(11)=(1/(3.6e-4))/(y(5)+y(6));
end
%x=zeros(11,1);
%x(1)=0.748;x(2)=-0.747;x(3)=0.219;x(4)=1.195;x(5)=-0.219;x(6)=-0.748;x(7)=-1.195;x(8)=204.654;x(9)=406.997;x(10)=399.355;x(11)=204.654;
%[t,y]=ode45(@fangcheng,[0:1e-4:0.01],x);