Question: How do i plot in an procedure ?

Hi,

i have to make an iteration in the following way, but i have a problem to plot the result. Has someone an idea how to plot?

 

s_1 := 0.779e-1*t^2-.1345*t+.2031:
s_2 := 1.3502*t-6.8742:
s_3 := .1112*t^2-1.5108*t+11.527:
s_4 := 1.9322*t-15.099:
s_5 := -0.178e-1*t^2+2.648*t-22.312:
V_0 := .3*.15:
p_0 := 100000:

t:= 0:
tmp:= p_0:
S := 0:

while t do

if t >= 0 and t if t > 9.5 and t if t > 12.875 and t if t > 15.5 and t if t > 20.75 and t
V_1 := .3 * S:
K := tmp * 3:

p_1 := p_0 + V_1 / V_0 * K:
t1 := plot(p_0 + V_1 / V_0 * K,color=red):
p.t := display({p_0 + V_1 / V_0 * K})

tmp := p_1:
t := t + 0.125:
od;

with(plots) : display({t1});

Please Wait...