arxiv: v1 [math.nt] 16 Nov 2009

Similar documents
Expression for Restricted Partition Function through Bernoulli Polynomials

REGULARIZATION OF CERTAIN DIVERGENT SERIES OF POLYNOMIALS

GENERALIZED HARMONIC NUMBER IDENTITIES AND A RELATED MATRIX REPRESENTATION

The r-generalized Fibonacci Numbers and Polynomial Coefficients

Math 155 (Lecture 3)

1 6 = 1 6 = + Factorials and Euler s Gamma function

Harmonic Number Identities Via Euler s Transform

Recurrence Relations

The Arakawa-Kaneko Zeta Function

Series with Central Binomial Coefficients, Catalan Numbers, and Harmonic Numbers

q-lucas polynomials and associated Rogers-Ramanujan type identities

A Note on Generating Functions and Summation Formulae for Meixner Polynomials of One and Two variables

A GENERALIZATION OF THE SYMMETRY BETWEEN COMPLETE AND ELEMENTARY SYMMETRIC FUNCTIONS. Mircea Merca

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

Zeros of Polynomials

FLOOR AND ROOF FUNCTION ANALOGS OF THE BELL NUMBERS. H. W. Gould Department of Mathematics, West Virginia University, Morgantown, WV 26506, USA

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

The Riemann Zeta Function

Optimally Sparse SVMs

A symmetrical Eulerian identity

1 Generating functions for balls in boxes

Dr. Clemens Kroll. Abstract

NICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =

An enumeration of flags in finite vector spaces

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

2.4 - Sequences and Series

WHAT ARE THE BERNOULLI NUMBERS? 1. Introduction

Math 475, Problem Set #12: Answers

Some results on the Apostol-Bernoulli and Apostol-Euler polynomials

SOME TRIBONACCI IDENTITIES

ECEN 644 HOMEWORK #5 SOLUTION SET

and Genocchi Polynomials

Frequency Response of FIR Filters

The log-behavior of n p(n) and n p(n)/n

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

4 The Sperner property.

CHAPTER 4 BIVARIATE DISTRIBUTION EXTENSION

Appendix to Quicksort Asymptotics

Application to Random Graphs

#A51 INTEGERS 14 (2014) MULTI-POLY-BERNOULLI-STAR NUMBERS AND FINITE MULTIPLE ZETA-STAR VALUES

Some Basic Diophantine Equations

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

Evaluation of Some Non-trivial Integrals from Finite Products and Sums

A Note on the Symmetric Powers of the Standard Representation of S n

arxiv: v2 [math.nt] 9 May 2017

New Generalization of Eulerian Polynomials and their Applications

Random Models. Tusheng Zhang. February 14, 2013

On forward improvement iteration for stopping problems

Chapter 2: Numerical Methods

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

A q-analogue of some binomial coefficient identities of Y. Sun

A MODIFIED BERNOULLI NUMBER. D. Zagier

Lecture 7: Properties of Random Samples

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 19 11/17/2008 LAWS OF LARGE NUMBERS II THE STRONG LAW OF LARGE NUMBERS

10. Comparative Tests among Spatial Regression Models. Here we revisit the example in Section 8.1 of estimating the mean of a normal random

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

Introductions to PartitionsP

Roger Apéry's proof that zeta(3) is irrational

EXPANSION FORMULAS FOR APOSTOL TYPE Q-APPELL POLYNOMIALS, AND THEIR SPECIAL CASES

A symbolic approach to multiple zeta values at the negative integers

Determinants of order 2 and 3 were defined in Chapter 2 by the formulae (5.1)

Super congruences concerning Bernoulli polynomials. Zhi-Hong Sun

An Asymptotic Expansion for the Number of Permutations with a Certain Number of Inversions

The Random Walk For Dummies

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

Subject: Differential Equations & Mathematical Modeling-III

Some remarks for codes and lattices over imaginary quadratic

SOME TRIGONOMETRIC IDENTITIES RELATED TO POWERS OF COSINE AND SINE FUNCTIONS

The multiplicative structure of finite field and a construction of LRC

The Method of Least Squares. To understand least squares fitting of data.

SOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0,

Binomial transform of products

Assignment Number 3 Solutions

Research Article Some E-J Generalized Hausdorff Matrices Not of Type M

A Mean Paradox. John Konvalina, Jack Heidel, and Jim Rogers. Department of Mathematics University of Nebraska at Omaha Omaha, NE USA

PROBLEM SET 5 SOLUTIONS 126 = , 37 = , 15 = , 7 = 7 1.

Asymptotic analysis of the Bell polynomials by the ray method

Benaissa Bernoussi Université Abdelmalek Essaadi, ENSAT de Tanger, B.P. 416, Tanger, Morocco

Week 5-6: The Binomial Coefficients

Some p-adic congruences for p q -Catalan numbers

Linear Regression Demystified

Sums Involving Moments of Reciprocals of Binomial Coefficients

A combinatorial contribution to the multinomial Chu-Vandermonde convolution

subcaptionfont+=small,labelformat=parens,labelsep=space,skip=6pt,list=0,hypcap=0 subcaption ALGEBRAIC COMBINATORICS LECTURE 8 TUESDAY, 2/16/2016

INFINITE SEQUENCES AND SERIES

Introducing a Novel Bivariate Generalized Skew-Symmetric Normal Distribution

ON SOME DIOPHANTINE EQUATIONS RELATED TO SQUARE TRIANGULAR AND BALANCING NUMBERS

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Bernoulli Polynomials Talks given at LSBU, October and November 2015 Tony Forbes

Course : Algebraic Combinatorics

This section is optional.

Introduction To Discrete Mathematics

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Chapter 6 Principles of Data Reduction

Generalized Little q-jacobi Polynomials as Eigensolutions of Higher-Order q-difference Operators

Generating Functions for Laguerre Type Polynomials. Group Theoretic method

arxiv: v1 [math.co] 5 Jan 2014

Bertrand s Postulate

Lecture Overview. 2 Permutations and Combinations. n(n 1) (n (k 1)) = n(n 1) (n k + 1) =

On the distribution of coefficients of powers of positive polynomials

Transcription:

Complete Bell polyomials ad ew geeralized idetities for polyomials of higher order arxiv:0911.3069v1 math.nt] 16 Nov 2009 Boris Rubistei, Stowers Istitute for Medical Research 1000 50th St., Kasas City, MO 64110, U.S.A. November 16, 2009 Abstract The relatios betwee the Beroulli ad Euleria polyomials of higher order ad the complete Bell polyomials are foud that lead to ew idetities for the Beroulli ad Euleria polyomials ad umbers of higher order. Geeral form of these idetities is cosidered ad geeratig fuctio for polyomials satisfyig this geeral idetity is foud. 1 Itroductio The history of the Beroulli polyomials B (s) couts more tha 250 years, as L. Euler first studied them for arbitrary values of the argumet. He itroduced the Euler polyomials E (s), the geeralizatio of these polyomials to the so-called Euleria polyomials H (s,ρ) was made by Frobeius 5]. The Bell polyomials 4] also appeared to be useful i combiatorial applicatios. I 1920s N. Nörlud 7] itroduced the Beroulli B (m) (s d) ad Euler polyomials E (m) (s d) of higher order addig parameters d = {d i }, 1 i m. Similar extesio for the Euleria polyomials (s,ρ d) was made by L. Carlitz i 3]. The Beroulli ad Euleria polyomials of higher order appear to be useful for the descriptio of the restricted partitio fuctio defied as a umber of iteger oegative solutios to Diophatie equatio d x = s (see 9]). The author showed i 10] that the restricted partitio fuctio ca be expressed through the Beroulli polyomials of higher order oly. It appears that there exist relatios that coect the Beroulli, Euler ad related polyomials ad their higher order geeralizatios to the complete Bell polyomials. These relatios preseted i Sectio 3 give rise to ew idetities for the polyomials of higher order, ad the idetity for the Euleria polyomials is discussed i Sectio 4. I Sectio 5 we itroduce a more geeral form of the abovemetioed idetities ad fid the geeratig fuctios for the polyomials satisfyig these relatios. The Sectio 6 is devoted to derivatio of ew relatios for the Beroulli polyomials of higher order. 1

2 Beroulli ad Euleria polyomials ad their geeralizatio to higher order The Beroulli B (s) ad Euler E (s) polyomials are defied through the correspodig geeratig fuctios te st e t 1 = B (s) t!, 2e st e t + 1 = E (s) t!. (1) =0 Frobeius 5] studied the so-called Euleria polyomials H (s,ρ) satisfyig the geeratig fuctio =0 =0 e st 1 ρ e t ρ = H (s,ρ) t, (ρ 1), (2)! which reduces to defiitio of the Euler polyomials at fixed value of the parameter E (s) = H (s, 1). The Apostol-Beroulli 1] ad Apostol-Euler 6] are the particular cases of the Euleria polyomials. Itroduce the geeralized Euleria polyomials through the geeratig fuctio It is easy to see that e st1 αρ e t ρ = =0 H (s,ρ,α) t!. (3) B (s) = H 1 (s,1,0), E (s) = H (s, 1,1), H (s,ρ) = H (s,ρ,1). Itroduce the geeralized Euleria polyomials of higher order H (m) (s,ρ,α d) as follows e st m 1 α i ρ i e d it ρ i = =0 H (m) (s,ρ,α d) t!. (4) The geeratig fuctio for the Beroulli polyomials of higher order B (m) (s d) reads 2]: m e st d i t e dit 1 = B (m) (s d)t!. (5) It ca be checed by compariso that =0 B (m) (s d) =!π 1 ( m)! H(m) m (s,1,0 d), π 1 = m d i. It is easy to see that the Euler E (m) (s d) 2] ad Euleria polyomials of higher order H (m) (s, ρ d) 3] ca be writte as E (m) (s d) = H(m) (s, 1,1 d), (6) (s,ρ d) = H (m) (s, ρ, 1 d). 1=1 2

3 Beroulli ad Euleria polyomials of higher order as Bell polyomials The complete Bell polyomials B (a 1,a 2,...) 8] are defied through the geeratig fuctio ( ) a i B i (a 1,a 2,...) exp i! ti = t i. (7) i! Cosider Taylor series i=0 dt l e dt 1 = ( 1) i 1di B i t i i i!, where B i is the Beroulli umber. Usig the above expasio i (5) we have ( ) e st exp ( 1) i 1B iσ i t i = B i (s + a 1,a 2,...) ti i i! i!, (8) where a i = ( 1) i 1 B i σ i /i ad σ = σ (d) = m d i deotes a power sum, so that Settig here m = 1, d = 1 we fid i=0 B (m) (s d) = B (s + a 1,a 2,...). (9) B (s) = B (s + a 1,a 2,...), Cosider Taylor series for the Euleria polyomials a i = ( 1) i 1B i i. l 1 ρ e dt ρ = dt 1 ρ i=2 ρd i H i 1 (ρ) t i 1 ρ i!. where H i (ρ) = H i (0,ρ) is the Euleria umber. I a particular case with ρ i = ρ i (4) usig the above expasio we have ( exp st σ ) ( ) 1t ρh i 1 (ρ)σ i t i exp = B i (s σ 1 + a 1,a 2,...) ti 1 ρ 1 ρ i! i!, (10) where a i = ρh i 1 (ρ)σ i /(1 ρ). Comparig (10) to (7) we obtai i=2 i=0 (s,ρ d) = B (s σ 1 + a 1,a 2,...). (11) The relatios (9,11) are the particular cases of a geeral idetity (s,ρ,α d) = B (s σ 1 + a 1,a 2,...) (12) with a i = ρh i 1 (ρ,α)σ i /(1 αρ) ad H i (ρ,α) = H i (0,ρ,α). Note that the value of the polyomial is equal to the complete Bell polyomial of the same order with parameters proportioal to the value of the correspodig polyomial at zero. 3

4 Geeralized Euleria polyomials idetities Cosider a ormalized distributio P(t) that ca be characterized by the sequece of momets µ or cumulats κ. The momets geeratio fuctio is P(t) itself ad cumulats geeratio fuctio is l P(t), i.e., P(t) = 1 + =1 µ! t, l P(t) = =1 κ! t. (13) It is worth to ote that i the theory of symmetric fuctios µ play role of the complete homogeeous symmetric polyomials, while κ stad for the symmetric power sums. The above defiitios imply the relatio betwee the momets ad cumulats 1 ( ) 1 µ = κ + κ µ. (14) 1 =1 O the other had it follows from (13) ad (7) that so that µ = B (κ 1,κ 2,...), (15) 1 ( ) 1 B (κ 1,κ 2,...) = κ + κ B (κ 1,κ 2,...). (16) 1 =1 Usig the relatio (12) ad the defiitio (7) oe fids: l =0 (s,ρ,α d) t = (s σ 1 )t! The last relatio leads to the followig idetity: H (m) (s,ρ,α d) = (s σ 1 ) 1 (s,ρ,α d) ρ 1 αρ that is rewritte as +1 (s,ρ,α d) = (s σ 1) (s,ρ,α d) ρ 1 αρ ρ 1 αρ =1 =0 H i 1 (0,ρ,α)σ i i! t i. (17) ( ) 1 H 1 (0,ρ,α)σ 1 (s,ρ,α d), ( ) H (0,ρ,α)σ +1 (s,ρ,α d). (18) Settig here m = 1, d 1 = 1 we fid for the geeralized Euleria polyomials H +1 (s,ρ,α) = (s 1)H (s,ρ,α) ρ ( ) H (0,ρ,α)H (s,ρ,α). (19) 1 αρ 5 Geeral case Cosider a more geeral problem of costructio of polyomials P (s,r) that ca be represeted by the Bell complete polyomials with parameters a i beig proportioal to the values of the same polyomial at some s = r, that satisfy the geeralized versio of (17) =0 l =0 P (s,r) t = Q 1 (s,r)t + Q 2 (s,r)! P i 1 (r,r) i! t i, (20) 4

where Q 1 (s,r) ad Q 2 (s,r) are polyomials i s. The we fid P +1 (s,r) = Q 1 (s,r)p (s,r) + Q 2 (s,r) =0 ( ) P (r,r)p (s,r). (21) We wat to fid the coditio o the geeratig fuctio G(s,r,t) for the polyomials P (s,r). Multiplyig the equatio (21) by t /! ad summig over we obtai P +1 (s,r) ( ) ] t t = Q 1 (s,r)p (s,r) + Q 2 (s,r) P (r,r)p (s,r)!!. =0 =0 The r.h.s. of the above equatio evaluates to Q 1 (s,r)g(s,r,t)+q 2 (s,r)g(r,r,t)g(s,r,t), while the l.h.s. is obtaied by differetiatio of G(s,r,t) w.r.t. t, so that G(s,r,t) Cosider a particular case of (22) i the form G(s,r,t) =0 = Q 1 (s,r)g(s,r,t) + Q 2 (s,r)g(r,r,t)g(s,r,t). (22) = Q 1 (s)g(s,r,t) + Q 2 (r)g(r,r,t)g(s,r,t). (23) The polyomials P (s,r) geerated by G(s,r,t) satisfy the idetity P +1 (s,r) = Q 1 (s)p (s,r) + Q 2 (r) F(r,t) =0 ( ) P (r,r)p (s,r). (24) We loo for solutio i the form G(s,r,t) = H(s,t)F(r,t) that leads to the equatio = Q 1 (s) l H(s,t) ] F(r,t) + Q 2 (r)h(r,t)f 2 (r,t). (25) The term i the square bracets i the last equatio should be idepedet of s, so that the followig coditio holds l H(s,t) Q 1 (s) = M (t) for some fuctio M(t) defied up to arbitrary costat. The coditio leads to relatio ad the equatio (25) reduces to the Beroulli equatio H(s,t) = e Q 1(s)t M(t), (26) which has the solutio The solutio of (23) reads F(r,t) = M (t)f(r,t) + Q 2 (r)e Q 1(r)t M(t) F 2 (r,t), (27) F(r,t) = Q 1 (r)e M(t) CQ 1 (r) Q 2 (r)e Q 1(r)t. (28) G(s,r,t) = H(s,t)F(r,t) = 5 Q 1 (r)e Q 1(s)t CQ 1 (r) Q 2 (r)e Q 1(r)t, (29)

where the itegratio costat C is determied from the coditio Fially we obtai G(s,r,0) = G 0 (s,r) = G(s,r,t) = Q 1 (r) CQ 1 (r) Q 2 (r). (30) G 0 (s,r)q 1 (r)e Q 1(s)t Q 1 (r) + G 0 (s,r)q 2 (r)(1 e Q 1(r)t ). (31) Settig i the last equatio r = 0, Q 1 (s) = s 1, Q 2 (0) = ρ/(1 αρ) ad G 0 (s,0) = (1 αρ)/(1 ρ) we obtai the geeratig fuctio i the l.h.s. of (3). It is ot difficult to geeralize the above geeratig fuctio to the case of polyomials of higher order. We defie the geeratig fuctio G(s,r,t,d) for the polyomials of higher order P (m) (s, r d) as follows G(s,r,t,d) = e Q 1(s) Q 1 (0)]t The polyomial idetity reads: P (m) +1 (s,r d) = (Q 1(s) Q 1 (0)(1 σ 1 ))P (m) (s,r d) + Q 2 (r) The case r = s leads to the geeral solutio G(s,t) = m e Q 1(0)d 1 t G 1 0 (s,r) + Q 2(r)(1 e Q 1(r)d i t )/Q 1 (r). (32) =0 ( ) P (r,r)σ +1 P (m) (s,r d). (33) G 0 (s)q 1 (s)e Q 1(s)t Q 1 (s) + G 0 (s)q 2 (s)(1 e Q 1(s)t ). (34) It should be oted that the solutios (31,34) are valid oly whe Q 1 (s) is idetically ozero. Otherwise the equatio (22) reduces to that ca be writte i the form G(s,r,t) F(r,t) = Q 2 (s,r)g(r,r,t)g(s,r,t), (35) = G(r,r,t), F(r,t) = l G(s,r,t) Q 2 (s,r). From the defiitio of F(r,t) we fid G(s,r,t) = expq 2 (s,r)f(r,t)] ad obtai G(r,r,t) = expq 2 (r,r)f(r,t)]. Thus we arrive at the followig equatio for F(r,t) which has the solutio F(r,t) = exp(q 2 (r,r)f(r,t)), F(r,t) = lq 2(r,r)(C t), Q 2 (r,r) ad we obtai the solutio of equatio (35) G(s,r,t) = exp Q ] 2(s,r) Q 2 (r,r) l Q 2(r,r)(C t) = Q 2 (r,r)(c t)] Q 2(s,r)/Q 2 (r,r). (36) 6

Settig i (36) t = 0 we obtai the itegratio costat C = 1 Q 2 (r,r) G Q 2(r,r)/Q 2 (s,r) 0, G 0 = G(s,r,0). Fially, we arrive at the solutio ] G(s,r,t) = G 0 1 Q 2 (r,r)tg Q 2(r,r)/Q 2 (s,r) Q2 (s,r)/q 2 (r,r). (37) For r = s equatio (35) reduces to G(s, t) The solutio i this case is foud from (34) i the limit Q 1 (s) 0 G(s,t) = 0 = Q 2 (s)g 2 (s,t). (38) G 0 (s) 1 G 0 (s)q 2 (s)t. (39) 6 Idetities for Beroulli polyomials ad umbers Usig the relatio (9) ad the defiitio (7) oe fids: l =0 B (m) (s d) t = st +! The last relatio leads to the followig idetity: B (m) (s d) = sb (m) 1 (s d) 1 =1 ( 1) i 1 B iσ i i! i t i. (40) ( ) ( 1) B σ B (m) (s d). (41) Settig here m = 1, d 1 = 1 we fid for the Beroulli polyomials B (s) = sb 1 (s) 1 ( ) ( 1) B B (s). (42) =1 =1 Fially settig i the above relatio s = 0 we arrive at the idetity for the Beroulli umbers ( ) ( 1) B B = B. (43) As all odd Beroulli umbers (except B 1 ) vaish, oe ca fid ( ) 1 ( 1) ]B σ B (m) (s d) = σ 1B (m) 1 (s d), (44) =1 from what it follows that B (m) (s d) = (s σ 1)B (m) 1 (s d) 1 B (s) = (s 1)B 1 (s) 1 ( B = B 1 1 =1 7 =1 ( ) B σ B (m) (s d), (45) =1 ( ) B B (s), (46) ) B B, (47)

where the relatio (47) is the well-ow Euler idetity ad (46) is a particular case of the sum idetity for the Beroulli polyomials where =0 ( ) B (s)b (t) = (s + t 1)B 1 (s + t) ( 1)B (s + t). Rewrite the idetity (45) usig the defiitio of the power sum as follows: B (m) (s d) = (s σ 1)B (m) 1 (s d) 1 = (s σ 1 )B (m) 1 (s d) 1 =1 ( ) B d i B(m) (s d) B (m) (s d) + = (s σ 1 )B (m) 1 (s d) + m B(m) (s d) 1 The relatio (48) produces two ew idetities ad B (m+1) d i = d d i. (s d i ) = (s σ 1 )B (m) 1 =0 ( ) ] B σ B (m) (s d) B (m+1) (s d i ), (48) (s d) ( m)b(m) (s d), (49) B m (m+1) (s d i ) = m(s σ 1 )B (m) m 1 (s d). (50) Recallig the idetity for the Beroulli polyomials of higher order 7]: we have B (m+1) B (m+1) (s + d i d i ) B (m+1) (s d i ) = B (m) 1 (s d), (s + d i d i ) = sb (m) 1 (s d) ( m)b(m) (s d). (51) Cosider a particular case of all d i = 1, i.e., d = 1; otig that σ 1 = m we fid from (49) the idetity derived by Nörlud i 7] mb (m+1) (s 1) = (m )B (m) (s 1) + (s m)b (m) 1 (s 1). (52) It is easy to show that for arbitrary ozero λ the followig relatio is valid: mb (m+1) (s λ 1) = (m )B (m) From the last relatio for m = it follows that (s λ 1) + (s λm)b(m) (s λ 1). (53) B m (m+1) (s λ 1) = (s λm)b (m) m m 1 (s λ 1) = (s λi). (54) 1 8

Acowledgemet The author thas L. Fel for umerous fruitful discussios. Refereces 1] Apostol T.M.: O the Lerch Zeta Fuctio, Pacific J. Math. 1 161-167 (1951) 2] Batema, H., Erdelýi, A.: Higher Trascedetal Fuctios, vol. 1. McGraw-Hill, New Yor (1953) 3] Carlitz, L.: Euleria Numbers ad Polyomials of Higher Order. Due Math. J. 27, 401-423 (1960) 4] Comtet, L.: Advaced Combiatorics, Chapter 2. Reidel, Dordrecht (1974) 5] Frobeius, F.G.: Über die Beroullische Zahle ud die Eulerische Polyome. Sitzugsberichte der Köiglich Preußische Aademie der Wisseschafte zu Berli, 809-847 (1910) 6] Luo, Q.-M., Srivastava H.M.: Some Geeralizatios of the Apostol-Beroulli ad Apostol-Euler Polyomials, J. Math. Aal. Appl. 308, 290-302 (2005) 7] Nörlud, N.E.: Mémoire sur les Polyomes de Beroulli. Acta Math. 43, 121-196 (1922) 8] Riorda, J.: A Itroductio ot Combiatorial Aalysis. Wiley Publicatios i Mathematical Statistics. Wiley/Chapma ad Hall, Lodo (1958) 9] Rubistei, B.Y., Fel, L.F.: Restricted Partitio Fuctio as Beroulli ad Euleria Polyomials of Higher Order. Ramauja J. 11, 331-347 (2006) 10] Rubistei B.Y.: Expressio for Restricted Partitio Fuctio through Beroulli Polyomials, Ramauja J. 15, 177-185 (2008) 9