equacion diferencial

Programa en Scilab

s=%s

num=2;

den=s^3*(s^2+2*s+10);

g=syslin('c',num/den);

gf=tf2ss(g);

se=pfss(gf);

Resultado en Scilab:

 se  =
 
 
       se(1)
 
                        2  
    0.2 - 0.04s - 0.012s   
    --------------------   
              3            
             s             
 
       se(2)
 
    0.064 + 0.012s   
    --------------   
                2    
     10 + 2s + s

Lineas a añadir en Scilab

r=roots(denom(se(1)));

a(3)=horner(s^3*se(1),r(1));

a(2)=horner(derivat(s^3*se(1)),r(1));

a(1)=horner(derivat(derivat(s^3*se(1))),r(1))/2;

for k=1:3,
ds1(k)=a(k)/s^k,
end;

Resultado en Scilab del desarrollo de se(1)

 ds1  =
 
 
       ds1(1)
 
            
   -0.012   
   -------  
            
     s      
 
       ds1(2)
 
           
   -0.04   
   ------  
      2    
     s     
 
       ds1(3)
 
         
   0.2   
   ----  
     3   
    s