Algorithmic Manipulation of Fibonacci Identities

Size: px
Start display at page:

Download "Algorithmic Manipulation of Fibonacci Identities"

Transcription

1 Algorthmc Manpulaton of Fbonacc Identtes 1 Stanley Rabnowtz 1 Vne Brook Road Westford, MA USA 1 Introducton Methods for manpulatng trgonometrc expressons, such as changng sums to products, changng products to sums, expandng functons of multple angles, etc, are wellknown [1] In fact, the process of verfyng trgonometrc denttes s algorthmc (see [] or [5] Roughly speakng, all trgonometrc denttes can be derved from the basc dentty sn x + cos x =1 Methods for manpulatng expressons nvolvng Fbonacc and Lucas numbers are also known; however, they do not seem to be systematcally collected n one place The fact that verfyng Fbonacc denttes s algorthmc seems to be less well-known, especally consderng the strong relatonshp there s between Fbonacc numbers and trgonometrc functons (see equaton (7 below For example, the problem of establshng the dentty F 14r (F 7 n+4r + F 7 n 7F 10r (F 6 n+4rf n + F n+4r F 6 n+1f 6r (F 5 n+4rf n + F n+4rf 5 n 35F r (F 4 n+4rf 3 n + F 3 n+4rf 4 n=f 7 4rF 7n+14 (1 s stll consdered hard as wtnessed by the fact that ths problem appeared n the advanced problem secton of The Fbonacc Quarterly [10] Yet ths dentty (and all such denttes can be proven usng a straghforward algorthm that we shall descrbe below Roughly speakng, all Fbonacc denttes can be derved from the basc dentty L n 5Fn =4( 1 n The Basc Algorthms The key dea to algorthmcally provng denttes nvolvng polynomals n F an+b and L an+b s to frst reduce them to polynomals n F n and L n Todothat, we need reducton formulas for F m+n and L m+n These are well-known In fact, all the formulas presented n ths paper can be found n the lterature (see, for example, [9] Most of these formulas were known to Lucas n the 19th century [6] ALGORITHM FbEvaluate TO NUMERICALLY EVALUATE F k AND L k : Gven an nteger constant k, toevaluate F k or L k numercally, apply the followng algorthm: STEP 1: [Make subscrpt postve] If k<0, apply algorthm FbNegate (gven below STEP : [Recurse] If k>1, apply the recurson: F n = F n 1 + F n L n = L n 1 + L n Reprnted from Applcatons of Fbonacc Numbers, volume 6 (edted by G E Bergum, et al, Kluwer Academc Publshers, Dordrecht, The Netherlands: 1996, pp (

2 Repeat step untl k < Step reduces the subscrpt by 1, so the recurson must eventually termnate STEP 3: [Intal values] If k s 0 or 1, then use the followng ntal values: F 0 =0, F 1 =1; L 0 =, L 1 =1 (3 Reprnted from Applcatons of Fbonacc Numbers, volume 6, (edted by G E Bergum, et al, pp Corrected 6/98 NOTE: Whle ths may not be the fastest way to evaluate F k and L k,tsnevertheless an effectve algorthm It s our purpose here to present effectve methods to show that the manpulatons descrbed are all algorthmc and t s not our purpose to fnd the most effcent algorthms possble Faster algorthms can be obtaned by usng the double argument formulas gven below n dsplay (1 For a descrpton and analyss of fast methods of numercally evaluatng F k and L k, see chapter 7 of [7] ALGORITHM FbNegate TO REMOVE NEGATIVE SUBSCRIPTS: F n =( 1 n+1 F n L n =( 1 n L n (4 ALGORITHM FbReduce TO REMOVE SUMS IN SUBSCRIPTS: F m+n = F ml n + L m F n L m+n = 5F mf n + L m L n These are also called the addton formulas 3 The Fundamental Identty Connectng F and L The Fbonacc and Lucas numbers are connected by the well-known dentty: (5 L n 5F n =4( 1 n (6 In the same way that all trgonometrc denttes are consequences of the fundamental dentty sn x + cos x =1,all Fbonacc/Lucas denttes are consequences of the fundamental dentty (6 That s, equaton (6 s the unque dentty connectng F n and L n The proof s a consequence of the fact that F n = 5 n+1 sn( n ln L n = n cos( n ln or F n = ( n n snh 5 ln ( n L n = n cosh ln (7

3 so that the result s equvalent to the result for trgonometrc polynomals, whose proof can be found n [] The fundamental dentty allows us to take any polynomal n F n and L n and remove all powers of L n (of degree or hgher by usng the followng algorthm ALGORITHM RemovePowersOfL : Use the dentty: L n =5Fn +4( 1 n (8 Contnue applyng ths substtuton untl no L n term has an exponent larger than 1 Smlarly, we could remove powers of F n nstead ALGORITHM RemovePowersOfF : Use the dentty: Fn = L n 4( 1 n (9 5 Contnue applyng ths substtuton untl no F n term has an exponent larger than 1 4 The Smplfcaton Algorthm Let us be gven a polynomal functon of elements of the form F x and L x, where the subscrpts of F and L are of the form a 1 n 1 + a n + + a k n k + b where b and the a are nteger constants and the n are varables To put ths expresson n canoncal form, we apply the followng algorthm: ALGORITHM FbSmplfy TO TRANSFORM AN EXPRESSION TO CANON- ICAL FORM STEP 1: [Remove sums n subscrpts] Apply algorthm FbReduce to remove any sums (or dfferences n subscrpts STEP : [Make multplers postve] All subscrpts are now of the form cx where c s an nteger For any term n whch the multpler c s negatve, apply algorthm FbNegate After performng ths step, you must take care to replace any expressons of the form ( 1 an+b where a and b are constants by a properly evaluated [( 1 a ] n ( 1 b,toavod wndng up wth expressons lke ( 1 n n your fnal canoncal form STEP 3: [Remove multplers] All subscrpts are now of the form cx where c s a postve nteger For any term n whch the multpler c s not 1, apply algorthm FbReduce successvely untl all subscrpts are varables STEP 4: [Evaluate constants] If any term nvolves an expresson of the form Fc k or L k c where c s an nteger constant, use algorthm FbEvaluate to replace F c and L c by ther numercal equvalents STEP 5: [Remove Fbonacc Powers] If any term nvolves an expresson of the form Fn k where k > 1 and n s a varable, apply algorthm RemovePowersOfF to leave only lnear terms n F n PROVING IDENTITIES To prove that an expresson s dentcally 0, apply algorthm FbSmplfy The expresson s dentcally 0 f and only f algorthm FbSmplfy transforms t to 0 3

4 Ths s because algorthm FbSmplfy leaves us wth a polynomal n varables of the form F n and L n wth the degree of F n beng 0 or 1 Convertng to trgonometrc form usng formula (7, we get a polynomal n varables of the form sn cx and cos cx wth the sne terms havng degree 0 or 1 Such polynomals cannot be dentcally 0 snce t s known that all trgonometrc denttes are consequences of the fundamental dentty sn θ+cos θ =1(see, for example, [5] and so some sne term would have to have degree or more Example Let us see how to prove F n+w F n+w L w F n+w F n w F n w F n w =(L 3w L w F n wth w odd, an dentty that comes from [11] Algorthmc proof Snce we are gven that w s odd, we make the substtuton w =m+1 before we begn We must show that F n+4m+ F n+m+1 L m+1 F n+m+1 F n m 1 F n m 1 F n 4m Fn(L 6m+3 L m+1 sdentcally 0 Frst we apply algorthm FbReduce (step 1 to get: (Fn(80F m +16L m 80F 6m 3L 6m ((F n ( 5F 4m +3L 4m +(3F 4m L 4m L n (5F m F n F n L m F m L n + L m L n (5F m + L m (5F m F n F n L m F m L n + L m L n (5F m F n + F n L m + F m L n + L m L n +(5F m F n + F n L m + F m L n + L m L n (F n (5F 4m +3L 4m +(3F 4m + L 4m L n /16 Then we apply algorthm FbNegate (step to get: (80F m Fn + 15FmF 3 n + 50F m F 4m Fn 80F 6m Fn +16FnL m +75FmF nl m +10F 4m FnL m +15F m FnL m + FnL 3 m+30f m FnL 4m +6FnL m L 4m 3FnL 6m 5FmL 3 n+6f m F 4m L n 11FmL m L n+ 6F 4m L m L n 7F m L ml n L 3 ml n +F m L 4m L n +L m L 4m L n/16 Applyng algorthm FbReduce agan (step 3 to get rd of scalar multples n subscrpts yelds: (80FmF n 15FmF 6 n + 160F m FnL m 50FmF 5 nl m +16FnL m + 5FmF 4 nl m + 100FmF 3 nl 3 m +5FmF nl 4 m 10F m FnL 5 m FnL 6 m/3 Havng evaluated constants lke F 1 all along, we proceed to step 5 and apply algorthm RemovePowersOfF Upon expandng and gatherng lke terms together, we fnd that the result s 0 Thus our example s n fact an dentty 5 Other Algorthms Sometmes we want to transform an expresson from one form to another, rather than put t n canoncal form The followng algorthms can be used for such purposes ALGORITHM ConvertToF TO REMOVE LUCAS NUMBERS: Use the dentty: L n = F n 1 + F n+1 (10 ALGORITHM ConvertToL TO REMOVE FIBONACCI NUMBERS: Use the dentty: F n = L n 1 + L n+1 (11 5 4

5 5 THE DOUBLE ARGUMENT FORMULAS Lettng m = n n formula (5 gves us the followng formulas F n = F n L n L n = 5F n + L n These correspond to the double angle formulas n trgonometry TO REMOVE SCALAR MULTIPLES OF ARGUMENTS IN SUBSCRIPTS: Repeatedly apply the reducton formulas successvely In ths manner, we obtan the denttes: F 3n = 3F nl n +5Fn 3 4 L 3n = L3 n +15FnL n 4 F 4n = F nl n (5Fn + L n L 4n = 5F (13 n 4 +30FnL n + L 4 n 8 F 5n = 5F n(5fn 4 +10FnL n + L 4 n 16 SHORTCUT: Apply the followng recurrences: L 5n = 15F nl 4 n +50FnL 3 n + L 5 n 16 (1 F kn = F (k 1nL n + L (k 1n F n L kn = 5F (k 1nL n + L (k 1n F n Even more straghtforward s to use a drect formula: (14 Example: F kn = 1 k 1 L kn = 1 k 1 k 1 k/ ( k +1 ( k 5 Fn L k n 5 Fn +1 L k 1 n (15 56F 9n = 65F 9 n F 7 nl n F 5 nl 4 n + 40F 3 nl 6 n +9F n L 8 n 56L 9n = 565F 8 nl n F 6 nl 3 n F 4 nl 5 n + 180F nl 7 n + L 9 n

6 Because of the fundamental dentty, these can also be put n other forms, such as the followng, whch converts drectly to our canoncal form 6 F kn = F n L kn = k/ k 1 k k ( k 1 ( k ( 1 (n+1 L k 1 n ( 1 (n+1 L k n (16 In general, F kn cannot be replaced by sums of powers of F n alone, so no formula of ths type s gven Also of nterest are the followng compact formulas F kn = L kn = k ( k F nfn 1 k F k ( k F nfn 1 k L (17 ALGORITHM FbExpand TO TURN PRODUCTS INTO SUMS: F m F n = L m+n ( 1 n L m n 5 L m L n = L m+n +( 1 n L m n F m L n = F m+n +( 1 n F m n (18 To remove products of more than two terms, the above expanson formulas can be used repeatedly, expandng out the results after each step Removng powers s the same as removng products For example, L n = L n L n = L n +( 1 n L 0 Smlarly for F nsowehave F n = L n ( 1 n 5 L n = L n +( 1 n (19

7 7 Contnung n ths way, we fnd: F 3 n = F 3n 3( 1 n F n 5 L 3 n = L 3n +3( 1 n L n F 4 n = L 4n 4( 1 n L n +6 5 L 4 n = L 4n +4( 1 n L n +6 F 5 n = F 5n 5( 1 n F 3n +10F n 5 L 5 n = L 5n +5( 1 n L 3n +10L n F 6 n = L 6n 6( 1 n L 4n +15L n 0( 1 n 15 L 6 n = L 6n +6( 1 n L 4n +15L n + 0( 1 n (0 To remove hgher powers, the expanson formulas (18 can be used repeatedly, expandng out the results after each step CHANGE OF BASIS (Shft Formulas ALGORITHM FbShft TO TRANSFORM AN EXPRESSION INVOLV- ING F n, L n INTO ONE INVOLVING F n+a, L n+b : ( Fn L n ( Lb F = a L a L b 5F a F b 5F b L a ( Fn+a L n+b (1 Proof Solve the lnear equatons F n+a = 1 (F al n + L a F n L n+b = 1 (5F bf n + L b L n for F n and L n In a smlar manner, we fnd ( ( Fn Fb F = a L n F a L b L a F b L b L a ( Fn L n = ( Fn+a F n+b ( ( ( Lb L a Ln+a (3 5(L a F b F a L b 5F b 5F a L n+b To change from an arbtrary bass to another, apply algorthm FbReduce to transform the gven expresson to the bass (F n,l n Then use one of the above shft formulas

8 ALGORITHM ConvertToAlphaBeta TO EXPRESS F n AND L n IN TERMS OF α AND β: Use the well-known Bnet forms: F n = αn β n 5 (4 L n = α n + β n to express F n and L n n terms of α and β, the roots of the characterstc equaton x = x+1 ALGORITHM RemoveAlphaBeta TO REMOVE α s AND β s: Sometmes a formula nvolves the quanttes α and β Totransform such expressons nto canoncal form, apply the denttes: α n = L n + F n 5 β n = L n F n 5 Expand out the resultng polynomal and express t n the form x + y 5 If y s not dentcally 0, then the resultng expresson cannot be expressed n terms of Fbonacc and Lucas numbers alone The quanttes α and β can be rentroduced nto the resultng expresson by means of the formula 5=α β (6 8 (5 and ether α or β may be removed by use of one of the denttes α + β =1 or αβ = 1 (7 whchever s preferred TO REMOVE POWERS OF α AND β: The followng denttes are frequently useful: α n = αf n + F n 1 β n = βf n + F n 1 (8 THE SUBTRACTION FORMULAS: Combnng the reducton and negaton formulas gves us the followng: F m n =( 1 n FmL n F n L m L m n =( 1 n LmL n 5F m F n

9 9 6 Arbtrary Startng Condtons Let {H n } be asequence that satsfes H n+ = H n+1 + H n (9 wth ntal condtons H 0 = a and H 1 = b We can formally manpulate ths sequence by frst convertng to Fbonacc numbers and then convertng back ALGORITHM ConvertHToF TO EXPRESS H n IN TERMS OF F n AND F n 1 : Use the dentty: H n = af n 1 + bf n (30 To prove an dentty nvolvng H n, use equaton (30 to convert H s to F s, and then apply algorthm FbSmplfy After applyng algorthms ConvertHToF and FbSmplfy to obtan a canoncal form, use the followng algorthm to transform the result (expressed wth F s and L s back to an expresson nvolvng H s ALGORITHM ConvertToH TO EXPRESS F n and L n IN TERMS OF H n AND H n+1 : where e = b a ab F n = 1 e (bh n ah n+1 L n = 1 e [(b ah n (a bh n+1 ] ALGORITHM HExpand TO EXPRESS H n+m IN TERMS OF H n AND H n+1 : (31 A more symmetrcal form s: H n+m = F m 1 H n + F m H n+1 (3 H n+m = 1 (F mg n + L m H n where G n = H n 1 + H n+1

10 CHANGE OF BASIS ALGORITHM HShft TO TRANSFORM AN EXPRESSION INVOLVING H n, H n+1 INTO ONE INVOLVING H n+c, H n+d : ( Hn = H n+1 1 F c 1 F d F c F d 1 ( Fd F c F d 1 F c 1 ( Hn+c H n+d 10 (33 To convert expressons nvolvng arbtrary subscrpts of H to expressons nvolvng H n+c and H n+d, apply algorthm ConvertHToF to put everythng n terms of F Then apply algorthm FbReduce to get everythng n terms of F n and L n Use algorthm ConvertToH to get the expresson n terms of H n and H n+1 ; and then perform an HShft to transform to the desred bass 7 Generalzed Fbonacc Numbers We may consder the sequences {u n } and {v n } defned by the recurrences: u 0 =0, u 1 =1, u n+ = Pu n+1 Qu n v 0 =, v 1 = P, v n+ = Pv n+1 Qv n (34 where P and Q are gven constants (usually ntegers wth Q 0 Let r 1 and r be the roots of the characterstc equaton so that and let so that x = Px Q (35 r 1 + r = P and r 1 r = Q, (36 r 1 = P + D D = P 4Q (37 and r = P D (38 We wll also assume that P and Q are chosen so that D 0 Ths assures that the roots are dstnct The Bnet form for u n and v n s u n = rn 1 r n r 1 r, v n = r n 1 + r n (39 Powers of r 1 and r can be removed from an expresson by means of the formulas: r1 n = v n + u n D r n = v n u n D (40

11 If a term of the form D s present, r 1 and r can be rentroduced, f desred, by the formula r 1 r = D (41 Also of nterest are the denttes: 11 r1 = Pr 1 Q r = Pr Q (4 and r n 1 = r 1 u n Qu n 1 r n = r u n Qu n 1 For these sequences, we have the followng algorthms (from [6]: ALGORITHM LucasNegate TO REMOVE NEGATIVE SUBSCRIPTS: u n = u n Q n v n = v (43 n Q n ALGORITHM LucasReduce TO REMOVE SUMS IN SUBSCRIPTS: u m+n = u mv n + u n v m v m+n = v mv n + Du m u n We also have the subtracton formulas : u m n = u mv n u n v m Q n v m n = v mv n Du m u n Q n (44 (45 THE FUNDAMENTAL IDENTITY: The fundamental dentty that connects u n and v n s: vn Du n =4Q n (46 Ths can be used to gve us the followng three algorthms ALGORITHM RemovePowersOfV : Apply the substtuton vn = Du n +4Q n (47

12 1 repeatedly untl the exponent of all v n terms s less than ALGORITHM RemovePowersOfU : Apply the substtuton u n = v n 4Q n D repeatedly untl the exponent of all u n terms s less than (48 ALGORITHM RemovePowersOfQ : Apply the substtuton Q n = v n Du n (49 4 ALGORITHM LucasExpand TO TURN PRODUCTS INTO SUMS: u m u n = v m+n Q n v m n D v m v n = v m+n + Q n v m n u m v n = u m+n + Q n u m n (50 Ths algorthm can be repeated to turn products of any number of u and v terms nto smple sums of such terms ALGORITHM ConvertToU TO REMOVE ALL v s: Use the dentty: v n = u n+1 Qu n 1 (51 ALGORITHM ConvertToV TO REMOVE ALL u s: Use the dentty: u n = v n+1 Qv n 1 (5 D RELATIONSHIP TO TRIGONOMETRIC EXPRESSIONS: To convert to trgonometrc form, use the denttes: ( n u n =Q n/ sn log r 1 / D r ( n v n =Q n/ cos log r 1 r (53

13 13 TO REMOVE SCALAR MULTIPLES IN SUBSCRIPTS: Repeatedly apply algorthm LucasReduce or use the denttes: u kn = u n v kn = k/ k 1 k k ( k 1 ( k ( 1 Q n v k 1 n ( 1 Q n v k n (54 These formulas convert drectly to our canoncal form In partcular, the double argument formulas are: u n = u n v n v n = v n Q n (55 Other formulas of nterest: u kn = 1 k 1 v kn = 1 k 1 k 1 k/ ( k +1 ( k D u n vn k D u +1 n vn k 1 u kn+s = v kn+s = k ( k u n( Qu n 1 k u +s k ( k u n( Qu n 1 k v +s TO REMOVE POWERS: Repeatedly apply algorthm LucasExpand or use the denttes: u k n = 1 v k n = 1 1 D k/ k k ( { k ( 1 Q n u(k n, f k s odd, v (k n, f k s even; ( k Q n v (k n (56 Any negatve subscrpts ntroduced can be removed by usng algorthm LucasNegate, f desred

14 14 CHANGE OF BASIS (Shft Formulas ALGORITHM LucasShft TO TRANSFORM AN EXPRESSION INVOLV- ING u n, v n INTO ONE INVOLVING u n+a, v n+b : ( ( ( un vb u = a un+a (57 v n v a v b Du a u b Du b v a v n+b PROVING IDENTITIES To effectvely prove an dentty nvolvng u s and v s, perform the followng algorthm whch s analogous to algorthm FbSmplfy ALGORITHM LucasSmplfy : Use algorthm LucasReduce to remove any sums n subscrpts Make all the multplers n subscrpts postve by usng algorthm LucasNegate Use LucasReduce agan to remove any constant multplers n the subscrpts Evaluate any numercal terms usng the recurrence (34 and ts ntal values (applyng algorthm LucasNegate frst for negatve subscrpts Fnally, use algorthm RemovePowersOfU and collect terms to arrve at a canoncal form Note that f the varables P, Q, and D all occur n the fnal expresson, you should replace D by ts equvalent value, P 4Q, from formula (37 The expresson s dentcally 0 f and only f ths canoncal form s 0 8 Arbtrary Second-Order Lnear Recurrences Let w n be an arbtrary second-order lnear recurrence wth constant coeffcents That s, w n = Pw n 1 Qw n, n (58 wth arbtrary ntal values w 0 and w 1, not both 0 Agan we requre that Q 0and P 4Q Snce u n and v n form a bass among all such sequences, w n can be expressed n terms of u n and v n Ths can be accomplshed by the followng formula: w n =(w 1 P w 0 u n + w 0 v n (59 Thus, to prove any dentty nvolvng w s, frst convert the w s to u s and v s and then apply algorthm LucasSmplfy An equvalent formula s w n =(w 1 Pw 0 u n + w 0 u n+1 = Qw 1 u n + w 0 u n+1 The reducton formula takes the form w n+m = Qw m 1 u n + w m u n+1 = Qw m u n 1 + w m+1 u n

15 9 Summatons We can perform ndefnte summatons of expressons nvolvng u n and v n any tme we can perform such summatons wth x n nstead, snce from the Bnet forms, these terms are actually exponentals wth bases r 1 and r ALGORITHM LucasSum : Frst, the expresson s converted to exponental form usng equaton (39 Then t s summed The result s converted back to u s and v s by usng equaton (40 Then r 1 and r are converted to expressons nvolvng P and Q usng equatons (38 and (36 The reader who has not kept abreast of developments n the area of symbolc algebra may be surprsed at the large class of functons that can now be automatcally summed Gosper s algorthm [3] can be used to sum complcated expressons nvolvng hypergeometrc seres Ths ncludes most of the common elementary functons See [4] pp 4 30 for an elementary exposton Example Let us see how to sum n k=1 L k cos kx algorthmcally Expandng L k by the Bnet form (4, and expandng cos kx by the formula cos x =(e x + e x /, turns the sum nto four terms each of whch s a smple geometrc progresson These are easly summed The result s turned back nto a nce form by usng formula (5 and the formula e x = cos x + sn x The result s the followng, whch we beleve to be new 15 n L k cos kx = k=1 1 [ 5+cosx 4cos x F n ( cos x 5 cos(n +1x cos x 3 ] + F n+1 (( cos x 3 cos nx +cos(n +x 4F n+ cos x cos nx (60 10 Dscoverng Identtes For the most part, the algorthms descrbed n ths paper can only be used to prove an dentty once the dentty s known or suspected It stll requres human ngenuty to dscover new and nterestng denttes There are, however, several ways these algorthms can be used to ad n dscoverng new denttes Algorthm LucasSum and ts Fbonacc counterpart, algorthm FbSum, are effectve algorthms for summng expressons Thus, they can be used, as n the prevous example, to dscover new results If the form of an dentty s suspected, algorthm FbSmplfy may frequently be used to advantage to dscover a new result For example, suppose that you suspect that F n+3 can be wrtten as a lnear combnaton of F n+, F n+1, and F ninthat case, you can use algorthm FbSmplfy to convert F n+3 af n+ bf n+1 cf n nto canoncal form The result s ( 1 n ( a + b +4c/5+L n(18 7a 3b c/10 + F n L n (8 3a b/

16 Snce the orgnal expresson s an dentty f and only f ths canoncal form s dentcally 0, we must have the three equatons 9a + b +4c =16 7a +3b +c =18 3a + b =8 Solvng these equatons yelds a =,b =,and c = 1 Ths reveals the dentty F n+3 F n+ F n+1 + F n =0 (61 11 Implementaton of the Algorthms All the algorthms descrbed n ths paper have been mplemented n Mathematca TM They are avalable from the author by emal Smlar algorthms exst for thrd-order lnear recurrences [8] and hgher-order recurrences 1 Phlosophcal Implcatons Now that these effectve algorthms for provng Fbonacc denttes are known, the phlosophcal queston arses as to whether one should study and/or publsh newly dscovered denttes The answer s yes Although formula (1 can now be proven by computer, the computer proof gves no nsght nto how the formula came about What s the sgnfcance of the bnomal coeffcents appearng n the formula? Does the formula generalze? Can a smlar proof be used to dscover new results? Should one stop submttng pretty new denttes to problem columns? The answer s no There s always the challenge nvolved n fndng an elegant proof For example, an algorthmc proof of the trgonometrc dentty sn 17x + cos 17x =1would not recognze that t follows from the Pythagorean Theorem Instead, t would expand out by the multple angle formula and get a horrendous mess whch t would then attempt to show s dvsble by (sn x + cos x 1 Elegant new denttes wll always be welcome to journals such as The Fbonacc Quarterly Dscoverng new denttes s stll mportant Ingenous proofs are always apprecated What these algorthms mply, s that the next tme you need a complcated dentty as a lemma n the mddle of a research paper, you can safely state the dentty and say that Ths dentty s straghtforward to prove References [1] Davd E Dobbs, Provng Trg Identtes to Freshpersons, The MATYC Journal 14 (1980:39 4 [] Davd E Dobbs and Robert Hanks, AModern Course on the Theory of Equatons Passac, NJ: Polygonal Publshng House, 1980 [3] R Wllam Gosper, Jr, Decson Procedure for Indefnte Hypergeometrc Summaton, Proceedngs of the Natonal Academy of Scences of the Unted States of Amerca 75 (1978:

17 [4] Ronald L Graham, Donald E Knuth, and Oren Patashnk, Concrete Mathematcs Readng, MA: Addson-Wesley Publshng Company, 1989 [5] M S Klamkn, On Provng Trgonometrc Identtes, Mathematcs Magazne 56 (1983:15 0 [6] Edouard Lucas, Théore des Fonctons Numérques Smplement Pérodques, Amercan Journal of Mathematcs 1 (1878:184 40, 89 31; reprnted as The Theory of Smply Perodc Numercal Functons, Santa Clara, CA: The Fbonacc Assocaton, 1969 [7] Roman Maeder, The Mathematca Programmer Boston, MA: Academc Press, 1994 [8] Stanley Rabnowtz, Algorthmc Manpulaton of Thrd-Order Lnear Recurrences, The Fbonacc Quarterly 34 (1996: [9] S Vajda, Fbonacc & Lucas Numbers, and the Golden Secton West Sussex, England: Ells Horwood Lmted, 1989 [10] Gregory Wulczyn, Problem H-34, The Fbonacc Quarterly 19 (1981:93 [11] Gregory Wulczyn, Problem B-464, The Fbonacc Quarterly 0 (198:370 AMS Classfcaton Numbers: 11Y16, 11B39, 11B37 17

arxiv: v1 [math.co] 12 Sep 2014

arxiv: v1 [math.co] 12 Sep 2014 arxv:1409.3707v1 [math.co] 12 Sep 2014 On the bnomal sums of Horadam sequence Nazmye Ylmaz and Necat Taskara Department of Mathematcs, Scence Faculty, Selcuk Unversty, 42075, Campus, Konya, Turkey March

More information

arxiv: v1 [math.ho] 18 May 2008

arxiv: v1 [math.ho] 18 May 2008 Recurrence Formulas for Fbonacc Sums Adlson J. V. Brandão, João L. Martns 2 arxv:0805.2707v [math.ho] 8 May 2008 Abstract. In ths artcle we present a new recurrence formula for a fnte sum nvolvng the Fbonacc

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

COMBINATORIAL IDENTITIES DERIVING FROM THE n-th POWER OF A 2 2 MATRIX

COMBINATORIAL IDENTITIES DERIVING FROM THE n-th POWER OF A 2 2 MATRIX COMBINATORIAL IDENTITIES DERIVING FROM THE n-th POWER OF A MATRIX J Mc Laughln 1 Mathematcs Department Trnty College 300 Summt Street, Hartford, CT 06106-3100 amesmclaughln@trncolledu Receved:, Accepted:,

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

= z 20 z n. (k 20) + 4 z k = 4

= z 20 z n. (k 20) + 4 z k = 4 Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5

More information

Remarks on the Properties of a Quasi-Fibonacci-like Polynomial Sequence

Remarks on the Properties of a Quasi-Fibonacci-like Polynomial Sequence Remarks on the Propertes of a Quas-Fbonacc-lke Polynomal Sequence Brce Merwne LIU Brooklyn Ilan Wenschelbaum Wesleyan Unversty Abstract Consder the Quas-Fbonacc-lke Polynomal Sequence gven by F 0 = 1,

More information

THE SUMMATION NOTATION Ʃ

THE SUMMATION NOTATION Ʃ Sngle Subscrpt otaton THE SUMMATIO OTATIO Ʃ Most of the calculatons we perform n statstcs are repettve operatons on lsts of numbers. For example, we compute the sum of a set of numbers, or the sum of the

More information

h-analogue of Fibonacci Numbers

h-analogue of Fibonacci Numbers h-analogue of Fbonacc Numbers arxv:090.0038v [math-ph 30 Sep 009 H.B. Benaoum Prnce Mohammad Unversty, Al-Khobar 395, Saud Araba Abstract In ths paper, we ntroduce the h-analogue of Fbonacc numbers for

More information

Binomial transforms of the modified k-fibonacci-like sequence

Binomial transforms of the modified k-fibonacci-like sequence Internatonal Journal of Mathematcs and Computer Scence, 14(2019, no. 1, 47 59 M CS Bnomal transforms of the modfed k-fbonacc-lke sequence Youngwoo Kwon Department of mathematcs Korea Unversty Seoul, Republc

More information

A combinatorial proof of multiple angle formulas involving Fibonacci and Lucas numbers

A combinatorial proof of multiple angle formulas involving Fibonacci and Lucas numbers Notes on Number Theory and Dscrete Mathematcs ISSN 1310 5132 Vol. 20, 2014, No. 5, 35 39 A combnatoral proof of multple angle formulas nvolvng Fbonacc and Lucas numbers Fernando Córes 1 and Dego Marques

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

The internal structure of natural numbers and one method for the definition of large prime numbers

The internal structure of natural numbers and one method for the definition of large prime numbers The nternal structure of natural numbers and one method for the defnton of large prme numbers Emmanul Manousos APM Insttute for the Advancement of Physcs and Mathematcs 3 Poulou str. 53 Athens Greece Abstract

More information

Determinants Containing Powers of Generalized Fibonacci Numbers

Determinants Containing Powers of Generalized Fibonacci Numbers 1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 19 (2016), Artcle 1671 Determnants Contanng Powers of Generalzed Fbonacc Numbers Aram Tangboonduangjt and Thotsaporn Thanatpanonda Mahdol Unversty Internatonal

More information

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of

More information

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017) Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed

More information

Beyond Zudilin s Conjectured q-analog of Schmidt s problem

Beyond Zudilin s Conjectured q-analog of Schmidt s problem Beyond Zudln s Conectured q-analog of Schmdt s problem Thotsaporn Ae Thanatpanonda thotsaporn@gmalcom Mathematcs Subect Classfcaton: 11B65 33B99 Abstract Usng the methodology of (rgorous expermental mathematcs

More information

Combinatorial Identities for Incomplete Tribonacci Polynomials

Combinatorial Identities for Incomplete Tribonacci Polynomials Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015, pp. 40 49 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM Combnatoral Identtes for Incomplete

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Affine transformations and convexity

Affine transformations and convexity Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/

More information

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016 U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

Hyper-Sums of Powers of Integers and the Akiyama-Tanigawa Matrix

Hyper-Sums of Powers of Integers and the Akiyama-Tanigawa Matrix 6 Journal of Integer Sequences, Vol 8 (00), Artcle 0 Hyper-Sums of Powers of Integers and the Ayama-Tangawa Matrx Yoshnar Inaba Toba Senor Hgh School Nshujo, Mnam-u Kyoto 60-89 Japan nava@yoto-benejp Abstract

More information

8.6 The Complex Number System

8.6 The Complex Number System 8.6 The Complex Number System Earler n the chapter, we mentoned that we cannot have a negatve under a square root, snce the square of any postve or negatve number s always postve. In ths secton we want

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

Fundamental loop-current method using virtual voltage sources technique for special cases

Fundamental loop-current method using virtual voltage sources technique for special cases Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41, The greatest common dvsor of two ntegers a and b (not both zero) s the largest nteger whch s a common factor of both a and b. We denote ths number by gcd(a, b), or smply (a, b) when there s no confuson

More information

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices Internatonal Mathematcal Forum, Vol 11, 2016, no 11, 513-520 HIKARI Ltd, wwwm-hkarcom http://dxdoorg/1012988/mf20166442 The Jacobsthal and Jacobsthal-Lucas Numbers va Square Roots of Matrces Saadet Arslan

More information

THE GENERALIZED (s, t)-fibonacci AND FIBONACCI MATRIX SEQUENCES

THE GENERALIZED (s, t)-fibonacci AND FIBONACCI MATRIX SEQUENCES TJMM 7 205), No 2, 37-48 THE GENERALIZED s, t)-fibonacci AND FIBONACCI MATRIX SEQUENCES AHMET İPEK, KAMIL ARI, AND RAMAZAN TÜRKMEN Abstract In ths paper, we study the generalzatons of the s, t)-fbonacc

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

1 Generating functions, continued

1 Generating functions, continued Generatng functons, contnued. Exponental generatng functons and set-parttons At ths pont, we ve come up wth good generatng-functon dscussons based on 3 of the 4 rows of our twelvefold way. Wll our nteger-partton

More information

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

COMPLEX NUMBERS AND QUADRATIC EQUATIONS COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not

More information

5 The Rational Canonical Form

5 The Rational Canonical Form 5 The Ratonal Canoncal Form Here p s a monc rreducble factor of the mnmum polynomal m T and s not necessarly of degree one Let F p denote the feld constructed earler n the course, consstng of all matrces

More information

Bezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0

Bezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0 Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the

More information

The KMO Method for Solving Non-homogenous, m th Order Differential Equations

The KMO Method for Solving Non-homogenous, m th Order Differential Equations The KMO Method for Solvng Non-homogenous, m th Order Dfferental Equatons Davd Krohn Danel Marño-Johnson John Paul Ouyang March 14, 2013 Abstract Ths paper shows a smple tabular procedure for fndng the

More information

Computing Correlated Equilibria in Multi-Player Games

Computing Correlated Equilibria in Multi-Player Games Computng Correlated Equlbra n Mult-Player Games Chrstos H. Papadmtrou Presented by Zhanxang Huang December 7th, 2005 1 The Author Dr. Chrstos H. Papadmtrou CS professor at UC Berkley (taught at Harvard,

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Digital Signal Processing

Digital Signal Processing Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

Some congruences related to harmonic numbers and the terms of the second order sequences

Some congruences related to harmonic numbers and the terms of the second order sequences Mathematca Moravca Vol. 0: 06, 3 37 Some congruences related to harmonc numbers the terms of the second order sequences Neşe Ömür Sbel Koaral Abstract. In ths aer, wth hels of some combnatoral denttes,

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

MA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials

MA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials MA 323 Geometrc Modellng Course Notes: Day 13 Bezer Curves & Bernsten Polynomals Davd L. Fnn Over the past few days, we have looked at de Casteljau s algorthm for generatng a polynomal curve, and we have

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

Randić Energy and Randić Estrada Index of a Graph

Randić Energy and Randić Estrada Index of a Graph EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL

More information

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0 Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

12. The Hamilton-Jacobi Equation Michael Fowler

12. The Hamilton-Jacobi Equation Michael Fowler 1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

SL n (F ) Equals its Own Derived Group

SL n (F ) Equals its Own Derived Group Internatonal Journal of Algebra, Vol. 2, 2008, no. 12, 585-594 SL n (F ) Equals ts Own Derved Group Jorge Macel BMCC-The Cty Unversty of New York, CUNY 199 Chambers street, New York, NY 10007, USA macel@cms.nyu.edu

More information

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem H.K. Pathak et. al. / (IJCSE) Internatonal Journal on Computer Scence and Engneerng Speedng up Computaton of Scalar Multplcaton n Ellptc Curve Cryptosystem H. K. Pathak Manju Sangh S.o.S n Computer scence

More information

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order: 68A Solutons to Exercses March 05 (a) Usng a Taylor expanson, and notng that n 0 for all n >, ( + ) ( + ( ) + ) We can t nvert / because there s no Taylor expanson around 0 Lets try to calculate the nverse

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

PHASED TILINGS AND GENERALIZED FIBONACCI IDENTITIES

PHASED TILINGS AND GENERALIZED FIBONACCI IDENTITIES appeared n: Fbonacc Quarterly 38(2000), pp. 282-288. PHASED TILINGS AND GENERALIZED FIBONACCI IDENTITIES Arthur T. Benjamn Dept. of Mathematcs, Harvey Mudd College, Claremont, CA 91711 benjamn@hmc.edu

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

Linear Algebra and its Applications

Linear Algebra and its Applications Lnear Algebra and ts Applcatons 4 (00) 5 56 Contents lsts avalable at ScenceDrect Lnear Algebra and ts Applcatons journal homepage: wwwelsevercom/locate/laa Notes on Hlbert and Cauchy matrces Mroslav Fedler

More information

Min Cut, Fast Cut, Polynomial Identities

Min Cut, Fast Cut, Polynomial Identities Randomzed Algorthms, Summer 016 Mn Cut, Fast Cut, Polynomal Identtes Instructor: Thomas Kesselhem and Kurt Mehlhorn 1 Mn Cuts n Graphs Lecture (5 pages) Throughout ths secton, G = (V, E) s a mult-graph.

More information

1 GSW Iterative Techniques for y = Ax

1 GSW Iterative Techniques for y = Ax 1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP C O L L O Q U I U M M A T H E M A T I C U M VOL. 80 1999 NO. 1 FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP BY FLORIAN K A I N R A T H (GRAZ) Abstract. Let H be a Krull monod wth nfnte class

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

The path of ants Dragos Crisan, Andrei Petridean, 11 th grade. Colegiul National "Emil Racovita", Cluj-Napoca

The path of ants Dragos Crisan, Andrei Petridean, 11 th grade. Colegiul National Emil Racovita, Cluj-Napoca Ths artcle s wrtten by students. It may nclude omssons and mperfectons, whch were dentfed and reported as mnutely as possble by our revewers n the edtoral notes. The path of ants 07-08 Students names and

More information

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,

More information

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars

More information

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Perron Vectors of an Irreducible Nonnegative Interval Matrix Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of

More information

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results. Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson

More information

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve

More information

a b a In case b 0, a being divisible by b is the same as to say that

a b a In case b 0, a being divisible by b is the same as to say that Secton 6.2 Dvsblty among the ntegers An nteger a ε s dvsble by b ε f there s an nteger c ε such that a = bc. Note that s dvsble by any nteger b, snce = b. On the other hand, a s dvsble by only f a = :

More information

Random Walks on Digraphs

Random Walks on Digraphs Random Walks on Dgraphs J. J. P. Veerman October 23, 27 Introducton Let V = {, n} be a vertex set and S a non-negatve row-stochastc matrx (.e. rows sum to ). V and S defne a dgraph G = G(V, S) and a drected

More information

An efficient algorithm for multivariate Maclaurin Newton transformation

An efficient algorithm for multivariate Maclaurin Newton transformation Annales UMCS Informatca AI VIII, 2 2008) 5 14 DOI: 10.2478/v10065-008-0020-6 An effcent algorthm for multvarate Maclaurn Newton transformaton Joanna Kapusta Insttute of Mathematcs and Computer Scence,

More information

The Expectation-Maximization Algorithm

The Expectation-Maximization Algorithm The Expectaton-Maxmaton Algorthm Charles Elan elan@cs.ucsd.edu November 16, 2007 Ths chapter explans the EM algorthm at multple levels of generalty. Secton 1 gves the standard hgh-level verson of the algorthm.

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

From Biot-Savart Law to Divergence of B (1)

From Biot-Savart Law to Divergence of B (1) From Bot-Savart Law to Dvergence of B (1) Let s prove that Bot-Savart gves us B (r ) = 0 for an arbtrary current densty. Frst take the dvergence of both sdes of Bot-Savart. The dervatve s wth respect to

More information

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy Comparatve Studes of Law of Conservaton of Energy and Law Clusters of Conservaton of Generalzed Energy No.3 of Comparatve Physcs Seres Papers Fu Yuhua (CNOOC Research Insttute, E-mal:fuyh1945@sna.com)

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys

More information

Quantum Mechanics for Scientists and Engineers. David Miller

Quantum Mechanics for Scientists and Engineers. David Miller Quantum Mechancs for Scentsts and Engneers Davd Mller Types of lnear operators Types of lnear operators Blnear expanson of operators Blnear expanson of lnear operators We know that we can expand functons

More information

332600_08_1.qxp 4/17/08 11:29 AM Page 481

332600_08_1.qxp 4/17/08 11:29 AM Page 481 336_8_.qxp 4/7/8 :9 AM Page 48 8 Complex Vector Spaces 8. Complex Numbers 8. Conjugates and Dvson of Complex Numbers 8.3 Polar Form and DeMovre s Theorem 8.4 Complex Vector Spaces and Inner Products 8.5

More information

Graph Reconstruction by Permutations

Graph Reconstruction by Permutations Graph Reconstructon by Permutatons Perre Ille and Wllam Kocay* Insttut de Mathémathques de Lumny CNRS UMR 6206 163 avenue de Lumny, Case 907 13288 Marselle Cedex 9, France e-mal: lle@ml.unv-mrs.fr Computer

More information

Math 261 Exercise sheet 2

Math 261 Exercise sheet 2 Math 261 Exercse sheet 2 http://staff.aub.edu.lb/~nm116/teachng/2017/math261/ndex.html Verson: September 25, 2017 Answers are due for Monday 25 September, 11AM. The use of calculators s allowed. Exercse

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information