matlab ode45 初值,关于ODE45初值问题和erf函数的问题

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

在求解一微分方程时候，采用了ode45 当使用第二种初值的时候，可以完美运行求解，当使用第一种初值的时候，就会报错

报错和主程序如下， 希望大家能帮帮我

————————————————————————————————————

错误使用 erf

输入必须为实数完全数。

出错 new_gaussisan_00way_1 (line 10)

I101A=sqrt(2/PI)*(1/sqrt(x(4)))*(x(4)*(x(7)*x(9)+x(3))+x(2)*x(5))*exp(-x(8)^2/(2*x(4)))+(x(7)*x(5)+x(8)*x(3)+x(9)*x(2)+x(7)*x(8)*x(9))*erf(x(8)/(sqrt(2*x(4))));

出错 ode45 (line 263)

f(:,4) = feval(odeFcn,t+hA(3),y+f*hB(:,3),odeArgs{:});

出错 solve_00way_1 (line 12)

[tdata,xdata]=ode45('new_gaussisan_00way_1',tspan,X0);

————————————————————————————————————

主程序如下：

clear;clc;close all

t0=cputime;

tspan=0:0.01:10;

X0=[3.6672;5.5018;2.7499;8.2541;4.1256;2.0621;-1.915;-2.873;-1.436]; 【第一种初值】

X0=[0.02,0,0