Wolfram mathematica 求解随机微分方程

下面的代码用于在Mathematica中数值求解一个粒子的随机方程。我想知道是否有一种方法可以将它推广到一个以上粒子的情况，并将其平均化。有人知道怎么做吗

Clear["Global`*"]
{ a = Pi, , b = 2 Pi, l = 5, k = 1};
ic = x@tbegin == 1;
tbegin = 1;
tend = 400;
interval = {1, 25};
lst := NestWhileList[(# + RandomVariate@TruncatedDistribution[interval, 
StableDistribution[1, 0.3, 0, 0, 1]]) &, tbegin, # < tend &];
F[t_] := Piecewise[{{k, Or @@ #}}, -k] &[# <= t < #2 & @@@ 
Partition[lst, 2]];
eqn := x'[t] == (F@(t)) ;
sol = NDSolveValue[{eqn, ic}, x, {t, tbegin, tend}, 
MaxSteps -> Infinity];

Plot[sol@t, {t, tbegin, tend}]
First[First[sol]]
Plot[sol'[t], {t, tbegin, tend}]
Plot[F[t], {t, tbegin, tend}]

清除[“全局”*]
{a=Pi，b=2pi，l=5，k=1}；
ic=x@tbegin == 1;
tbegin=1；
倾向=400；
区间={1,25}；
lst:=嵌套列表[（#+RandomVariate@TruncatedDistribution[间隔，
稳定分布[1,0.3,0,0,1]]）和，tbegin，#F[t]：=分段[{{k，或@@@k}，-k]&[#在表达式中
eqn:=x'[t]=（F@（t）
有语法错误。你能纠正它吗？我纠正了。你现在能帮我吗？@DanielIn表达式
eqn:=x'[t]=（F@（t）
有语法错误。你能纠正吗？我纠正了。你现在能帮我吗？@Daniel