Dear Robert,

I hope you are doing well this sunday morning.

I am contacting you about this code I am trying to plot in 3D ( 3D surfce plot ).

Could help me please doing it ?

How can I convert this equation into matlab ?

imax:= 50; mu := 30; sigma := 5;
y[0] :=x -> exp(-(x-mu)^2/(2*sigma^2))/(sigma*sqrt(2*Pi));
beta := (m, l) -> (m+l )^3 ;
des1:= [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)), j = 1 .. imax)];
ics := {Seq(y(j)[0]=exp(-(j-mu)^2/(2*sigma^2))/(sigma*sqrt(2*Pi)); , j = 1 .. imax)};
S := dsolve(des union ics, numeric); plots[odeplot](S, [seq([t,j, y[j](t)], j = 1 ..imax)], t = 0 .. 200, axes = box, labels = [t, j, y[j]], orientation = [250, 60]);

Best regards !

sebgo

Dear Robert,

I hope you are doing well this sunday morning.

I am contacting you about  this code I am trying to plot in 3D ( 3D surfce plot).

Could help me please doing it ?

How can I convert this equation into matlab ?

imax:= 50; mu := 30; sigma := 5;

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

beta := (m, l) -> (m+l )^3 ;

des1:= [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)), j = 1 .. imax)];

ics := {Seq(y(j)[0]=exp(-(j-mu)^2/(2*sigma^2))/(sigma*sqrt(2*Pi)); , j = 1 .. imax)};

S := dsolve(des union ics, numeric); plots[odeplot](S, [seq([t,j, y[j](t)], j = 1 ..imax)], t = 0 .. 200, axes = box, labels = [t, j, y[j]], orientation = [250, 60]);

Best regards !

sebgo

this is the initial conditions

ics := {Seq(y(j)[0] = exp(-(j-mu)^2/(2*sigma^2))/(sigma*sqrt(2*Pi)), j = 1 .. imax)}

Sebgo

this is the initial conditions

ics := {Seq(y(j)[0] = exp(-(j-mu)^2/(2*sigma^2))/(sigma*sqrt(2*Pi)), j = 1 .. imax)}

Sebgo

Dear Alec

Initial conditions are here

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

for x=1:imax

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

end

Beta is a function 

beta(l+m)=(l+m)^3

SEBGO

Dear Alec

Initial conditions are here

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

for x=1:imax

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

end

Beta is a function 

beta(l+m)=(l+m)^3

SEBGO