A107389 - OEIS

This site is supported by donations to The OEIS Foundation.

A107389

Expansion of x*(1-6*x+7*x^2)/( (1-x)*(1+x)*(1-5*x+x^2)).

0, 1, -1, 2, 5, 31, 144, 697, 3335, 15986, 76589, 366967, 1758240, 8424241, 40362959, 193390562, 926589845, 4439558671, 21271203504, 101916458857, 488311090775, 2339638995026, 11209883884349, 53709780426727, 257339018249280, 1232985310819681, 5907587535849119 (list; graph; refs; listen; history; internal format)

	OFFSET	`0,4`
	LINKS	`Index to sequences with linear recurrences with constant coefficients, signature (5,0,-5,1).`
	MATHEMATICA	`m = 5 M = {{0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}, {1, m, 0, - m) Expand[Det[M - xIdentityMatrix[4]]] (-1 - 5 x + 5 x^3 + x^4) NSolve[Det[M - xIdentityMatrix[4]] == 0, x] v[1] = {0, 1, 1, 2}; v[n_] := v[n] = M . v[n - 1] digits = 50; aa = Table[Abs[v[n][[1]], {n, 1, digits}]` `Clear[M, m, v, aa] (A107389)m = 5; M = {{0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}, {1, m, 0, - m}}; Expand[Det[M - xIdentityMatrix[4]]] ; NSolve[Det[M - xIdentityMatrix[4]] == 0, x] ; v[1] = {0, 1, 1, 2}; v[n_] := v[n] = M . v[n - 1]; digits = 50; aa = Table[Abs[v[n][[1]]], {n, 1, digits}]`
	PROG	`(PARI) concat(0, Vec((1-6x+7x^2)/(1-x)/(1+x)/(1-5*x+x^2)+O(x^98))) \\ Charles R Greathouse IV, Jan 25 2012`
	CROSSREFS	`Sequence in context: A127298 A000133 A059086 * A189559 A077483 A119242` `Adjacent sequences: A107386 A107387 A107388 * A107390 A107391 A107392`
	KEYWORD	`sign,easy`
	AUTHOR	`Roger L. Bagula (rlbagulatftn(AT)yahoo.com), May 24 2005, corrected Sep 04 2008`
	EXTENSIONS	`Irregular sign at a(2) switched by R. J. Mathar, Jan 24 2012`

Content is available under The OEIS End-User License Agreement .

Last modified February 15 21:56 EST 2012. Contains 205860 sequences.