function test_forum18

model1
t=0:0.001:20; 
u=t.^2; 
P=0.1; 
set_param('model1/Gain', 'gain', num2str(P))
[t x y]=sim('model1',20,[], [t', u']);