ll@(a,b,c,d,e,n); Matrixm(4,4),v(4,1); m[0][0]=0; m[0][1]=1; m[0][2]=0; m[0][3]=0; m[1][0]=d; m[1][1]=c; m[1][2]=1; m[1][3]=0; m[2][0]=0; m[2][1]=0; m[2][2]=1; m[2][3]=0; m[3][0]=1; m[3][1]=0; m[3][2]=0; m[3][3]=1; v[0][0]=a; v[1][0]=b; v[2][0]=e; v[3][0]=0; v=m**(n+1)*v; wt(v[3][0]);