fourup.m

这是我和同事门一起努力编的的孤立子程序。基本上所以孤立子的情况都考虑到了

%根据k1和k2的值计算双孤子解v的表达式
%
%

clear all;

k1=1;
k2=1;
k3=1;
k4=1;

u ='-1/2*(k1*exp(z1)+k3*exp(z3)+k2*exp(z2)+exp(z2+2*z4+A24)*k2+exp(z2+z4+A24+z3)*k2+k4*exp(z4)+exp(z1+z2+A12)*k1+exp(z1+z2+A12)*k2+exp(z1+z3+A13)*k1+exp(z1+z3+A13)*k3+exp(z1+z4+A14)*k1+exp(z1+z4+A14)*k4+exp(z2+z3+A23)*k2+exp(z2+z3+A23)*k3+exp(z2+z4+A24)*k2+exp(z2+z4+A24)*k4+exp(z3+z4+A34)*k3+exp(z3+z4+A34)*k4+exp(2*z2+z4+A24)*k4+exp(z2+z4+A24+z1)*k2+exp(z2+z3+A23+z1)*k2+exp(z2+z4+A24+z1)*k4+exp(z1+z4+A14+z2)*k4+exp(2*z1+z4+A14)*k4+exp(z1+2*z4+A14)*k1+exp(z1+z4+A14+z3)*k1+exp(z1+z4+A14+z2)*k1+exp(z1+z3+A13+z4)*k3+exp(2*z1+z3+A13)*k3+exp(z1+z3+A13+z4)*k1+exp(z1+2*z3+A13)*k1+exp(z1+z3+A13+z2)*k1+exp(z1+z2+A12+z4)*k2+exp(z1+z2+A12+z3)*k2+exp(2*z1+z2+A12)*k2+exp(z1+z2+A12+z4)*k1+exp(z1+z4+A14+z3)*k4+exp(z2+z3+A23+z1)*k3+exp(2*z2+z3+A23)*k3+exp(z2+z3+A23+z4)*k3-exp(z2+z4+A24+z1)*k1-exp(z2+z3+A23+z4)*k4-exp(z2+z3+A23+z1)*k1-exp(z1+z4+A14+z3)*k3-exp(z1+z4+A14+z2)*k2-exp(z1+z3+A13+z4)*k4+exp(z1+z2+A12+z3)*k1+exp(z2+z3+A23+z4)*k2+exp(z1+2*z2+A12)*k1-exp(z1+z3+A13+z2)*k2-exp(z1+z2+A12+z4)*k4-exp(z1+z2+A12+z3)*k3-exp(z3+z4+A34+z1)*k1-exp(z2+z4+A24+z3)*k3-exp(z3+z4+A34+z2)*k2+exp(z1+z3+A13+z2)*k3+exp(2*z3+z4+A34)*k4+exp(z3+z4+A34+z2)*k4+exp(z3+z4+A34+z1)*k4+exp(z3+2*z4+A34)*k3+exp(z3+z4+A34+z2)*k3+exp(z3+z4+A34+z1)*k3+exp(z2+z4+A24+z3)*k4+exp(z2+2*z3+A23)*k2)/(exp(z1)+exp(z2)+exp(z3)+exp(z4)+exp(z1+z2+A12)+exp(z1+z3+A13)+exp(z1+z4+A14)+exp(z2+z3+A23)+exp(z2+z4+A24)+exp(z3+z4+A34))/(1+exp(z1)+exp(z2)+exp(z3)+exp(z4))'
u=subs(u,'k1',k1);
u=subs(u,'k2',k2);
u=subs(u,'k3',k3);
u=subs(u,'k4',k4);

u=subs(u,'2*A12','q2_12a');  
u=subs(u,'A12','q_12a');
u=subs(u,'2*A13','q2_13a');
u=subs(u,'A13','q_13a');
u=subs(u,'2*A14','q2_14a'); 
u=subs(u,'A14','q_14a');
u=subs(u,'2*A23','q2_23a');
u=subs(u,'A23','q_23a');
u=subs(u,'2*A24','q2_24a');  
u=subs(u,'A24','q_24a');
u=subs(u,'2*A34','q2_34a');  
u=subs(u,'A34','q_34a');


u=simplify(u)