lilei1.m

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

syms u r w s v;
estr='D2x=(u-r)*x+w*y,D2y=s*x+(v-r)*y';
result=dsolve(estr,'t');
result.x=simplify(result.x)
result.y=simplify(result.y)