如何用matlab求解不同参数值的ODE？

Question

只要参数p在一定范围内变化，我想求解简单的一维信号并获得不同的求解向量

dx / dt = 25.7 + x *（0.00038 * x-0.00000014 * x ^ 2）-x *（0.162 + 0.00000006 * x ^ 2）-p * x

这是我的代码

clear all

% Write here the initial conditions
xstart = 1200.0;

tspan = 0:1000;

%parameters
% s = 25.7;
% a = 0.00038;
% b = -0.00000014;
% c = 0.162;
% d = 0.00000006;
%p = 0.01;

p = 0:0.001:0.1;
x_end=zeros(length(tspan),length(p));

for ii = 1:length(p)
[t,x] = ode45(@(t,x)lv(t,x,p(ii)),tspan,xstart);
end 

和

<FormLabel>{I18n.strings("account.nickname")}<Text style={{ color: 'red' }}>*</Text></FormLabel>

我想收集针对任何p值获得的解向量。

非常感谢您的帮助或建议