如何用matlab求解带约束方程的微分方程组？

类似于如下，x1(0)=1,x2(0)=x3(0)=x4(0)=0
x1'=-3*x1*x2;
x2'=8*x3+x4;
x3'=7*x4;
x4'=-4*x1*x3
x1+x2+x3+x4=1，求程序应该如何编写？

这就是用matlab做数学题，你查有本书好象是matlab在高等数学中的应用，里面有例题