l1

0.67

:=

l2

0.33

:=

l3

0.15

:=

u1

25

:=

u2

u1

:=

u3

u1

:=

lr1

0.1

:=

lr2

lr1

:=

lr3

lr1

:=

Given

p1' t

( )

u2 p2 t

( )

⋅

u3 p3 t

( )

⋅

+

u1 p4 t

( )

⋅

+

l2

l3

+

lr1

+

(

) p1 t

( )

⋅

−

=

p1 0

( )

1

=

p2' t

( )

l2 p1 t

( )

⋅

u1

l1

+

l3

+

(

) p2 t

( )

⋅

−

=

p2 0

( )

0

=

p3' t

( )

l3 p1 t

( )

⋅

u3

l1

+

l2

+

(

) p3 t

( )

⋅

−

=

p3 0

( )

0

=

p4' t

( )

lr1 p1 t

( )

⋅

u1

l2

+

l3

+

(

) p4 t

( )

⋅

−

=

p4 0

( )

0

=

a

b

c

d





















Odesolve

p1

p2

p3

p4





















t

, 10

,





















:=

t1

0 0.01

,

10

..

:=

t2

0 1

,

10

..

:=

Rsys t1

( )

a t1

( )

b t1

( )

+

c t1

( )

+

d t1

( )

+

:=

Rsys t2

( )

1

0.983

0.966

0.949

0.932

0.916

0.9

0.884

0.868

0.853

0.838

=

0

2

4

6

8

10

0.8

0.85

0.9

0.95

1

Rsys t1

( )

Rsys t2

( )

t1 t2

,