prm=[2.5,4,1.3,0.8]
from SAU import *
xx=[0]+[1]*20 
SAUe=SAU(prm)  
yt=[]
for xt in xx:
   SAUe.zdn_zn(xt)
   SAUe.model()
   SAUe.otobraz()
   yt.append(SAUe.ypr[1])
import pylab
pylab.plot(yt)
pylab.show()