a=0*1:6; a[6]=1; 

p<-function(s){
1+s+s^2+s^3+s^4+s^5-4*s^6
}
n=50; 


for (i in 1:n){ # para o inverso a menor raíz em módulo 
    a[i+6]=(a[i+5]+a[i+4]+a[i+3]+a[i+2]+a[i+1]+a[i])/4
}
 
c=a[n-1]/a[n]; 1/c; 

p(1/c)

n=500; 
for (i in 1:n){ # para o inverso da maior raíz em módulo 
    a[i+6]=-(a[i+5]+a[i+4]+a[i+3]+a[i+2]-4*a[i+1]+a[i])
}
 
c=a[n-1]/a[n]; 1/c; 

p(1/c)