matlab ode 实数,关于ode45中erf函数（输入必须为实数完全数的报错问题）

最近在做有关ODE45的微分方程求解，在微分方程中，运用到了erf函数，虽然之前也出现过类似的问题，已经解决，但是这次丝毫找不到问题的原因所在，毫无头绪。

主程序如下：

[code]

%——————xp vp zp 分别为20个点位————————————————————————

xp=[-1.91545000000000;-1.68455000000000;-1.51545000000000;-1.28455000000000;-1.11545000000000;-0.884550000000000;-0.715450000000000;-0.484550000000000;-0.315450000000000;-0.0845500000000000;0.0845500000000000;0.315450000000000;0.484550000000000;0.715450000000000;0.884550000000000;1.11545000000000;1.28455000000000;1.51545000000000;1.68455000000000;1.91545000000000];

vp=[-2.87317500000000;-2.52682500000000;-2.27317500000000;-1.92682500000000;-1.67317500000000;-1.32682500000000;-1.07317500000000;-0.726825000000000;-0.473175000000000;-0.126825000000000;0.126825000000000;0.473175000000000;0.726825000000000;1.07317500000000;1.32682500000000;1.67317500000000;1.92682500000000;2.27317500000000;2.52682500000000;2.873175000