EULERINTEGRAL_f
//EULERINTEGRAL_f Scicos vector intgerator by Euler method
//Classical interface function ver1.0 - scilab-2.7
//17 novembre 2003 - IRCOM GROUP - Author : A.Layec
function [x,y,typ]=EULERINTEGRAL_f(job,arg1,arg2)
x=[];y=[];typ=[]
select job
case 'plot' then
standard_draw(arg1)
case 'getinputs' then
[x,y,typ]=standard_inputs(arg1)
case 'getoutputs' then
[x,y,typ]=standard_outputs(arg1)
case 'getorigin' then
[x,y]=standard_origin(arg1)
case 'set' then
x=arg1
graphics=arg1.graphics;exprs=graphics.exprs
model=arg1.model;
while %t do
text=['Set discrete integral parameters';'(Euler method)']
[ok,step,x0,exprs]=getvalue(text,..
['Step';'Initial state'],list('vec',1,'vec',-1),exprs)
nu=size(x0,'*')
if ~ok then break,end
graphics.exprs=exprs;
model.in=nu
model.out=nu
model.rpar=step
model.dstate=x0(:)
model.dep_ut=[%t %f] //for compatibility
x.graphics=graphics;x.model=model
break
end
case 'define' then
x0=0
step=1
nu=1
model=scicos_model()
model.sim=list('int_euler',2)
model.in=nu
model.out=nu
model.evtin=1
model.evtout=[]
model.rpar=step
model.dstate=x0
model.blocktype='d'
model.dep_ut=[%t %f]
exprs=[string(step);strcat(sci2exp(x0))]
gr_i=['thick=xget(''thickness'')'
'pat=xget(''pattern'')'
'fnt=xget(''font'')'
'xpoly(orig(1)+[0.7;0.62;0.549;0.44;0.364;0.291]*sz(1),orig(2)+[0.947;0.947;0.884;0.321;0.255;0.255]*sz(2),"'lines"')'
'xset(''font'',5,1)'
'txt=''Euler'';'
'style=5;'
'rectstr=stringbox(txt,orig(1),orig(2),0,style,1);'
'w=(rectstr(1,3)-rectstr(1,2))*%zoom;'
'h=(rectstr(2,2)-rectstr(2,4))*%zoom;'
'xstringb(orig(1)+sz(1)/2-w/2,orig(2)-h-4,txt,w,h,''fill'');'
'e=gce();'
'e.font_style=style;'
'xset(''thickness'',thick)'
'xset(''pattern'',pat)'
'xset(''font'',fnt(1),fnt(2))'
]
x=standard_define([2 2],model,exprs,gr_i)
end
endfunction