trapintegral_f.sci

numerical modulator with matlab

//TRAPINTEGRAL_f Scicos vector intgerator by Trapezoidal method
//Classical interface function ver1.0 - scilab-2.7
//17 novembre 2003 - IRCOM GROUP - Author : A.Layec
function [x,y,typ]=TRAPINTEGRAL_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';'(trapezoidal method)']
   [ok,step,x0_in,x0_out,exprs]=getvalue(text,...
      ['Step';'Initial input state';'Initial output sate'],...
      list('vec',1,'vec',-1,'vec',-1),exprs)

    nu=size(x0_in,'*')
    if ~ok then break,end
    if size(x0_in)==size(x0_out) then
         nu=size(x0_in,'*')
         if ok then
           graphics.exprs=exprs;
           model.in=nu
           model.out=nu
           model.rpar=step
           model.dstate=[x0_out(:);x0_in(:)]
           x.graphics=graphics;x.model=model
           break
         end
    else
         message('Vector Initial input state and vector Initial output sate must have the same size');
    end
  end

case 'define' then
  x0_in=0
  x0_out=0
  step=1
  nu=1

  model=scicos_model()
  model.sim=list('int_trap',2)
  model.in=nu
  model.out=nu
  model.evtin=1
  model.evtout=[]
  model.rpar=step
  model.dstate=[x0_out(:);x0_in(:)]
  model.blocktype='d'
  model.dep_ut=[%f %f]

  exprs=[string(step);string(x0_in);string(x0_out)]
  //gr_i=['xstringb(orig(1),orig(2),''  1/s  '',sz(1),sz(2),''fill'');']
  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(''thickness'',thick)'
        'xset(''pattern'',pat)'
        'xset(''font'',fnt(1),fnt(2))'
       ]
  x=standard_define([2 2],model,exprs,gr_i)
  x.graphics.id="Trapezoidal"
end
endfunction