我的m文件如下：
%heart.m
function [dxdt]=heart(t,x)
Ra=0.01;
Rm=0.01;
Rc=0.0398;
Rs=0.8738;
Da=1;
Dm=1;
Ls=0.001025;
Cs=2.896;
Cr=4.00;

Ep=0.06;
Vd=5.0;
Vo=15.0;
V=Vo-Vd;
%parameters

Em=3.0;   
%Max left ventricular elastance
Tc=0.8;   
%Cardiac cycle
T=0:0.001:1;
t=T./(0.2+0.1555*Tc); 
%Normalized time
En=1.553174.*(t./0.7).^(1.9)./(1+(t./0.7).^(1.9)).*1./(1+(t./1.173474).^(21.9));
%Normalized elastance function
Ea=Em*En;
E=Ea+Ep;

dxdt=zeros(size(x));
dxdt(1)=-Da.*x(2)+(Dm/Rm).*x(4)-(Dm/Rm).*(E.*(x(1)-Vd)-Ep*V);
dxdt(2)=-(Da*Ra+Rc)/Ls.*x(2)-(Da/Ls).*x(3)+(Da/Ls).*(E.*(x(1)-Vd)-Ep*V);
dxdt(3)=1/Cs.*x(2)-x(3)./(Rs*Cs)+x(4)./(Rs*Cs);
dxdt(4)=-x(3)./(Cr*Rs)-(1/(Cr*Rs)+Dm/(Cr*Rm)).*x(3)+Dm/(Cr*Rm).*(E.*(x(1)-Vd)-Ep*V);
return

function main
options=odeset('RelTol',1e-6);
X0=[68.36;0;90;5];
tspan=[1,5];
[t,X]=ode45(@heart,tspan,X0,options);
figure
hold on
plot(t,X(:,1))
return

求助大虾帮忙看看那里错了，运行后总显示
??? Input argument 'x' is undefined.
???  In an assignment  A(I) = B, the number of elements in B and
 I must be the same.

万分感谢呢~~