Dear all,

I tried unsuccessfully to plot  (in 3D) in this solution using:
imax=10, mu=5, sigma=1.5
I failed to write the initial conditions

Here is the maple code

Initial conditions

imax:=10;
mu:=5; sigma:=1.5;

y[0] := x -> exp(-(x-mu)^2/(2*sigma^2))/(sigma*sqrt(2*Pi));   

differential equation

  des:= [seq(diff(y[i](t), t) = y[i-1](t)*(sum(2^(j-i+1)*
    beta[i-1, j]*y[j](t), j = 1 .. i-1))+(1/2)*i*beta[i-1, i-1]
    *y[i-1](t)^2-y[i](t)*(sum(2^(j-i)*beta[i, j]*y[j](t),
    j = 1 .. i-1))-y[i](t)*(sum(beta[i, j]*y[j](t),
    j = 1 .. imax)), i = 1 .. imax)];

Could any one give me an help ?