
The input flow rate is given by
f
1
(t) ¼
0:01 ft
3
=s, t 0 s
0:05 ft
3
=s,0< t 0:01 s
0:01 ft
3
=s, t > 0:01 s
8
>
<
>
:
Numerical values of the system parameters are
A ¼ 0:0055 ft
2
, k ¼ 30 lb=ft, R ¼ 10
5
lb s=ft
2
, p
0
¼ 14:7 lb=in:
2
The system is at steady state prior to the pulse input in flow.
(a) Use the RK–Fehlberg integrator to simulate the transient respon se of p(t).
(b) Find the analytical solution for p(t).
(c) Find the solution for p(t) without the accumulato r present.
(d) Plot the responses from parts (a), (b), and (c) on the same graph.
(e) Simulate the response for p(t) with Simulink using the ‘‘ode45’’ integrator, and compare
the results with those