Wolfram Cloud Document

ClearAll["Global`*"];

α1=1;α2=1;α3=1;α4=1;α5=1;

g0=0.9;b10=0.03;b11=0.001;Topt=24;gama0=150;ds0=4;A0=0.2;w=2.10;z1=3.5;Oc=1.10;d1=0.66;d=0.4;gama1=2;d2=0.2;z=0.0019;β=1.024;γ=4;z0=10.10;T10=14.50;delta1=0.1;T0=20;b=1.30;Tmaks=35;

K=1;M=1201.;T=60.05;

x[t_]

0,0

:=10;

h1[t_]

0,0

:=28;

h2[t_]

0,0

:=1;

y[t_]

0,0

:=1.2;

z[t_]

0,0

:=0.25;

x[λ_]

m_

:=

K

∑

i=0

x[t]

m,i

*

i

λ

;

h1[λ_]

m_

:=

K

∑

i=0

h1[t]

m,i

*

i

λ

;

h2[λ_]

m_

:=

K

∑

i=0

h2[t]

m,i

*

i

λ

;

y[λ_]

m_

:=

K

∑

i=0

y[t]

m,i

*

i

λ

;

z[λ_]

m_

:=

K

∑

i=0

z[t]

m,i

*

i

λ

;

N

1

[λ_]

m_

:=g0*Exp-b*

h1[λ]

m

-Topt

Tmaks-Topt

+

z[λ]

m

-b10+b11*

h1[λ]

m

-Topt

z[λ]

m

+1

*

x[λ]

m

-

g0*

x[λ]

m

^2

gama0

;

N

2

[λ_]

m_

:=w*

h2[λ]

m

-z1*

h1[λ]

m

-T10+γ*z0-

y[λ]

m

;

N

3

[λ_]

m_

:=A0-delta1*

h2[λ]

m

;

N

4

[λ_]

m_

:=Oc-d1*

y[λ]

m

-d*

y[λ]

m

*

h2[λ]

m

;

N

5

[λ_]

m_

:=gama1*

h1[λ]

m

-T0

β

*

ds0

1+

h1[λ]

m

-Topt

-

z[λ]

m

-d2*

z[λ]

m

*

x[λ]

m

-z*

h1[λ]

m

-Topt;

A

1

m_,0

:=

N

1

[0]

m

;

A

2

m_,0

:=

N

2

[0]

m

;

A

3

m_,0

:=

N

3

[0]

m

;

A

4

m_,0

:=

N

4

[0]

m

;

A

5

m_,0

:=

N

5

[0]

m

;

DoFori=0,i<K,i++,

A

1

m,i

=Simplify

1

i!

*D

N

1

[λ]

m

,{λ,i}/.λ0;

A

2

m,i

=Simplify

1

i!

*D

N

2

[λ]

m

,{λ,i}/.λ0;

A

3

m,i

=Simplify

1

i!

*D

N

3

[λ]

m

,{λ,i}/.λ0;

A

4

m,i

=Simplify

1

i!

*D

N

4

[λ]

m

,{λ,i}/.λ0;

A

5

m,i

=Simplify

1

i!

*D

N

5

[λ]

m

,{λ,i}/.λ0;Do

x[t_]

m,k

=

1

Gamma[α1]

*Integrate

α1-1

(t-τ)

*

A

1

m,k-1

,{τ,0,t};

h1[t_]

m,k

=

1

Gamma[α2]

*Integrate

α2-1

(t-τ)

*

A

2

m,k-1

,{τ,0,t};

h2[t_]

m,k

=

1

Gamma[α3]

*Integrate

α3-1

(t-τ)

*

A

3

m,k-1

,{τ,0,t};

y[t_]

m,k

=

1

Gamma[α4]

*Integrate

α4-1

(t-τ)

*

A

4

m,k-1

,{τ,0,t};

z[t_]

m,k

=

1

Gamma[α5]

*Integrate

α5-1

(t-τ)

*

A

5

m,k-1

,{τ,0,t},{k,1,K};

f1[t_]

m_

:=

K

∑

i=0

x[t-(T/M)*m]

m,i

;

f2[t_]

m_

:=

K

∑

i=0

h1[t-(T/M)*m]

m,i

;

f3[t_]

m_

:=

K

∑

i=0

h2[t-(T/M)*m]

m,i

;

f4[t_]

m_

:=

K

∑

i=0

y[t-(T/M)*m]

m,i

;

f5[t_]

m_

:=

K

∑

i=0

z[t-(T/M)*m]

m,i

;

x[t_]

m+1,0

=

f1[(T/M)*(m+1)]

m

;

h1[t_]

m+1,0

=

f2[(T/M)*(m+1)]

m

;

h2[t_]

m+1,0

=

f3[(T/M)*(m+1)]

m

;

y[t_]

m+1,0

=

f4[(T/M)*(m+1)]

m

;

z[t_]

m+1,0

=

f5[(T/M)*(m+1)]

m

,{m,0,M-1};

X[t_]=X[t]=PiecewiseTable

f1[t]

m