program vaja1b " konstante modela constant p=1.0 " zacetne vrednosti constant e0=0.0, b0=0.0, c0=0.0 " casovna os constant tfin=10.0 " strukture modela ed = -1.0*e + 5.0/6.0*u e = integ(ed, e0) bd = - 3.0*b + 3.0/2.0*u b = integ(bd, b0) cd = - 4.0*c - 4.0/3.0*u c = integ(cd, c0) y = e + b + c u = step(t, p) termt t .ge. tfin cinterval ci=0.1 output 10, u, y prepar 2, u, y end
Vidimo, da je odziv enak prej"snjemu primeru, ki smo ga izra"cunali na druga"cen na"cin. S tem smo tudi potrdili pravilnost obeh metod.