如何用matlab求解微分方程组的数值解

如何用matlab求解微分方程组的数值解
原方程很复杂,而且有很多格式,在这里表示不出来.但总体格式如下
dx/dt=f(a,b,c,d,x,y,z);
dy/dt=g(a,b,c,d,x,y,z);
dz/dt=h(a,b,c,d,x,y,z);
其中,a,b,c,d都是系数.现在,计划给定a、b、c的值,用matlab分析x、y、z分别与t和d的关系,并画图（三维图）.需要用哪些函数?具体调用方式是什么?各参数是什么意义?（matlab帮助是英文的,很多看不懂）
原方程如下
ds=500-s-0.24*s/(80+s)*(xu+xw);
dxu=xu*(0.12*s/(80+s)-1.01)+xw*(b+60*s/(80+s)/(5500-xw))-a*xu*(5000-xw)/5000;
dxw=xw*(0.12*s*(5000-xw)/(80+s)/(5500-xw)-b-0.01)+a*xu*(5000-xw)/5000;
当a和b都固定时，s(t)、xu(t)、xw(t)的二维图都已画出。
现在就是要：先固定b=0.03，并假设a=0:0.1:1,再计算，并画出s关于a和t的三维图、xu关于a和t的三维图、xw关于a和t的三维图。怎么办？直接把a改为a=0:0.1:系统会报错。