!!! New Identity Discover !!!

gon · #1 18 มิถุนายน 2004, 16:56

รามานุจัน (Ramanujan) (1887 - 1920) ครั้งหนึ่งได้โพสต์ปัญหาไว้ใน Journal of The Indian Mathemetical Society :
จงพิสูจน์ว่า (cos 2p/9)^1/3+(cos 4p/9)^1/3-(cos p/9)^1/3 = [(3/2)(9^1/3-2)]^1/3

หรือ (cos 2p/9)^1/3+(cos 4p/9)^1/3+(cos 8p/9)^1/3 = [(3/2)(9^1/3-2)]^1/3

ซึ่งก็จะได้ว่า (sec 2p/9)^1/3+(sec 4p/9)^1/3-(sec p/9)^1/3 = [6(9^1/3-1)]^1/3 ด้วย

วันนี้(18 มิถุนายน 2547) ตอนบ่าย ผมได้ค้นพบ เอกลักษณ์ใหม่ ซึ่งผมยังไม่เคยเห็นที่ไหนมาก่อน ก็หวังว่ามันคงจะไม่ปรากฏใน Journal ที่ไหนสักแห่งในโลกนะ (เคยเจอมาทีแล้ว)

(cos 2p/7)^1/3+(cos 4p/7)^1/3+(cos 6p/7)^1/3 = {(1/2)[5 - 3(7^1/3)]}^1/3
และ
(sec 2p/7)^1/3+(sec 4p/7)^1/3+(sec 6p/7)^1/3 = [8 - 6(7^1/3)]^1/3

ก็ไม่รู้ว่าจะไปโพสต์ไว้ใน Journal ที่ไหน ของไทยเราก็ไม่รู้ว่ามีหรือเปล่า เลยมาบันทึกไว้ในบอร์ดแห่งนี้.
ใครว่าง ๆ ช่วยลองจิ้มเครื่องคิดเลขอีกทีว่ามันถูกไหม. ถ้าจะลองพิสูจน์ดูก็แล้วแต่ใจชอบครับ.

warut · #2 20 มิถุนายน 2004, 04:15

ขอแสดงความยินดีกับคุณ gon ด้วยครับ เอกลักษณ์ข้างต้นทุกอันนับเป็นเอกลักษณ์ที่
สวยงามมาก ดังนั้นไม่ว่าจะเป็น "discovery" หรือ "rediscovery" ก็เป็นเรื่องที่
น่าปิติยินดีกับผู้ที่ค้นพบด้วยตนเองเป็นอย่างยิ่ง ผมคงไม่มีความสามารถที่จะพิสูจน์
เอกลักษณ์เหล่านั้นได้ ดังนั้นถ้าคุณ gon จะมาแสดงเป็นขวัญตาให้กับชาว mathcenter
ชมก็จะดีมากครับ อีกอย่างนึงถ้าเกิดคุณ gon โชคดีได้เป็นคนพบคนแรกจริงๆแล้วเพื่อ
ที่จะให้ได้ "full credit" คุณ gon ต้องแสดงการพิสูจน์ด้วยครับ เพราะในยุคปัจจุบัน
เอกลักษณ์เช่นนี้สามารถถูกค้นพบได้ด้วย "computer search" ดังนั้นการพิสูจน์
ด้วย "มือ" จึงเป็นสิ่งสำคัญ (ถ้าเป็นยุค Ramanujan คงไม่เป็นปัญหา)

ขอเล่าตัวอย่างการค้นพบด้วย "computer search" ให้ฟังเล่นๆสักอันนึงนะครับ
แต่ก่อนเราคิดว่าการหา bit ที่ n ของ p ก็ไม่ได้ง่ายไปกว่าการหาค่าของ p ตั้งแต่
bit ที่ 1 ถึง n จนกระทั่งเมื่อไม่นานมานี้มีการค้นพบเอกลักษณ์ที่เกี่ยวกับค่า p ที่
ทำให้เราสามารถหาbit ที่ n ของ p ได้โดยตรง ไม่ต้องหากันตั้งแต่ bit ที่ 1 ถึง n
ที่สำคัญก็คือการค้นพบเอกลักษณ์ดังกล่าวเกิดขึ้นโดย "computer search" ก่อน
แล้วจึงมีการพยายามพิสูจน์อย่าง "rigorous" ด้วย "มือ" ต่ออีกทีครับ ตัวผมเองเชื่อ
ว่ายังมีเอกลักษณ์ที่สวยงามอีกจำนวนมากที่พบแล้วโดยคอมพิวเตอร์แต่ยังไม่มี
การเปิดเผยเพราะยังไม่สามารถพิสูจน์ได้

gon · #3 20 มิถุนายน 2004, 12:28

ครับ. ผมพิสูจน์ได้ แต่ยังมีข้อบกพร่องเล็กน้อยประมาณ 1% ที่ผมยังมองไม่ออกว่าผิดตรงไหน ที่จริงก็ไม่ได้ยุ่งยากมากมายอะไรหรอกครับ. ผมสนใจเรื่องเอกลักษณ์ทางตรีโกณเป็นพิเศษอยู่แล้ว อันนี้ก็ต่อยอดความรู้ผมออกไปอีก. เอกลักษณ์ง่าย ๆ อย่าง (sec 2p/7)³ + (sec4p/7)³+(sec 6p/7)³ = ... ผมก็คงไม่กล้ามาโพสต์ แต่อันนี้ผมว่าสวยมาก ด้วยความดีใจก็เลยรีบมาโพสต์ก่อนครับ. ตอนนี้ผมก็กำลังเมล์ไปถาม Dr. BRUCE C. BERNDT. อยู่ ไม่รู้ว่าจะออกหัวหรือก้อยครับ.
อ้อ. ขอเปลี่ยนนิดหน่อยเพื่อให้ทุกพจน์มีค่าเป็นบวก แบบนี้ครับ.

(cos p/7)^1/3 + (cos 3p/7)^1/3 - (cos 2p/7)^1/3 = { 0.5 [ 3ท7^1/3 - 5 ] }^1/3 and

(sec p/7)^1/3 + (sec 3p/7)^1/3 - (sec 2p/7)^1/3 = { 6ท7^1/3 - 8 }^1/3

gon · #4 24 มิถุนายน 2004, 15:33

เดี๋ยวจะมาพิสูจน์ให้ดูนะครับ. ตอนนี้กำลังสนุกสนานกับการควานหา Indentity อื่น ๆ อยู่ หมดกระดาษไปเกือบ 40 แผ่นแล้ว ชักอ่อนใจว่าจะไม่เจออีกอีกแล้ว

แต่คณิตศาสตร์นี่มันดีจริง ๆนะ ครับ. ยิ่งคิดยิ่งสนุก ได้ความรู้แนวคิดใหม่ ๆ เกิดขึ้นตลอดเวลา ยิ่งเป็นเรื่องที่เราชอบด้วยยิ่งสนุก ผมเล่นสมการกำลังสามจนคิดว่าไม่มีเหลือแล้วเอกลักษณ์สวย ๆ เขียนทั้งเป็น general case ออกมาก็ยังไม่เจอที่มันสวย ๆ ตอนนี้ขยับไปถึงกำลัง 5 กับ กำลัง 7 แล้วครับ.

ว่าแต่ว่า โดยรูปแบบ กำลัง 5 มันควรจะได้ ทำไปทำมา อ้าว สมการ พาราเมตริกซ์ไม่บังเกิดซะอีก เลยต้อง Solve ตรง ๆ Solve เสร็จอ้าวมันขัดแย้งกัน not valid เลยเป็นไปไม่ได้ สมการกำลัง 5 ง่าย ๆ ที่มองเห็นได้ก็คือ (ถ้าจำไม่ผิด) x⁵ - 5x³ + 5x - 1 = 0 ซึ่งก็จะมีรากทั้ง 5 ของสมการเป็น (ถ้าจำไม่ผิด) 2cosp/15, 2cos5p/15, 2cos7p/15, 2cos11p/15, 2cos13p/15 ใครมีสมการกำลังคี่อื่นสวย ๆ บอกผมด้วยนะครับ.

ไม่รู้ว่าจะเป็นเพราะ Abel Impossibility Theorem อะไรนั่นหรือเปล่าด้วย ทำให้มันหลุดออกมาเป็นรากที่ 5 หรือ รากที่ 3 ไม่ได้ หลายวันมานี้ หัวมึนตึ๊บตลอดครับ. ทั้งอาบน้ำ เข้าห้องน้ำ ทานข้าว คิดตลอด มีไอเดียใหม่ ๆ วิ่งเข้ามาในสมองตลอด แต่พอลงมือทำในกระดาษ ติดนิดติดหน่อยทุกทีเลย เดี๋ยวพรุ่งนี้หยุดคิดแล้วครับ. งานอื่นไม่เดินเลย

gon · #5 25 มิถุนายน 2004, 11:47

ท.บ. ให้ x₁, x₂, x₃ เป็นรากของสมการ 8x³ - ax² + bx - 1 = 0
จะได้ว่า x₁^1/3 + x₂^1/3 + x₃^1/3 = (a + 24 +6t)^1/3/2
และ (x₁x₂)^1/3 + (x₂x₃)^1/3 + (x₃x₁)^1/3 = (b + 12 + 3t)^1/3/2
โดยที่ : t³ - 3(a + 2b + 12)t - [ ab + 12(a + 2b + 6) ] = 0

พิสูจน์ : จะได้ว่า x₁x₂x₃ = 1/8
สมมติให้ y_i = x_i^1/3 ; i = 1, 2, 3 หรือ y = x^1/3 \ y₁y₂y₃ = (x₁x₂x₃)^1/3 = (1/8)^1/3 = 1/2
\ จะมีสมการ 2y³ - py² + qy - 1 = 0 ...(*)

จัดรูปจะได้ (2y³ - 1)³ = (py² - qy)³ แล้ว (2x - 1)³ - p³x² + q³x + 3pqx(2x - 1) = 0
ฎ 8x³ - (12 + p³ - 6pq)x² + (6 + q³ - 3pq)x - 1 = 0
เทียบสัมประสิทธิ์ จะได้ว่า

a = 12 + p³ - 6pq
b = 6 + q³ - 3pq ฎ 2b = 12 + q³ - 6pq

ให้ p³ = a + 24 + 6t ก็จะได้ว่า q³ = b + 12 + 3t และ pq = t + 6
แต่ p³q³ = (pq)³ ฎ (a + 24 + 6t)(b + 12 + 3t) = (t + 6)³ ซึ่งเมื่อกระจายแล้วจัดรูปจะได้ว่า
t³ - 3(a + 2b + 12)t - [ ab + 12(a + 2b + 6) ] = 0

จาก (*) จะได้ว่า y₁ + y₂ + y₃ = p/2 ฎ x₁^1/3 + x₂^1/3 + x₃^1/3 = (a + 24 +6t)^1/3/2
และ y₁y₂ + y₂y₃ + y₃y₁ = q/2 ฎ (x₁x₂)^1/3 + (x₂x₃)^1/3 + (x₃x₁)^1/3 = (b + 12 + 3t)^1/3/2

พิจารณาสมการ 8x³ + 4x² - 4x - 1 = 0 จะได้ว่า รากของสมการนี้คือ cos 2p/7, cos 4p/7, cos 6p/7
\ เมื่อ a = -4, b = -4 ก็จะได้ว่า t = -2(7)^1/3 เมื่อนำไปแทนใน ท.บ. ข้างต้น ก็จะได้เอกลักษณ์ทั้งสอง ออกมาตามลำดับ

ใครจะเอาแนวคิดดังกล่าวไปใช้ต่อ ก็เชิญตามสบายเลยครับ. ถ้าเจอของใหม่อีกยิ่งดี

TOP · #6 26 มิถุนายน 2004, 12:30

อีกไม่นานเราคงได้ รามานุจันแห่งประเทศไทยสักคนหนึ่งละ (ตอนนี้ก็ตามพิสูจน์ สูตรของรามานุจันตัวจริงไปก่อน

)

เรื่อง "computer search" ดูแล้วน่าสนใจมากครับ โดยเฉพาะเรื่องที่
ทำให้เรา "สามารถหาbit ที่ n ของ p ได้โดยตรง ไม่ต้องหากันตั้งแต่ bit ที่ 1 ถึง n" ไม่รู้ว่าค้นพบโดย computer ได้อย่างไร มีใครตั้งเป็นข้อคาดเดาบางอย่างไว้ก่อน แล้วให้ computer ตรวจสอบเล่นๆอย่างนั้นหรือครับ อีกประเด็นหนึ่งคือ ถึงแม้ว่าจะไม่ต้องเริ่มต้นหาจากบิตแรก แต่ว่าสุดท้ายแล้ว เราต้องคำนวณเยอะ พอๆกับเริ่มต้นหาจากบิตแรกรึเปล่า

#7 04 กรกฎาคม 2004, 00:39

เป็นเอกลักษณ์ที่น่าสนใจครับ ไปเจอที่เวบของแคนาดาเห็นว่าคล้ายๆกับของพี่กรก็เลยเอามาฝากครับ

tan²(p/7)+tan²(2p/7)+tan²(3p/7)=21

gon · #8 04 กรกฎาคม 2004, 14:01

noonuii นั่นเอง. ขอบคุณครับที่ไปหามาฝากอีก แต่อันนี้ (tan²p/7 + ... ) พี่เก็บไว้ใน collection ของพี่เรียบร้อยแล้ว. ถ้าเจอที่แปลก ๆ อีกก็บอกมาอีกได้นะครับ.

#9 06 กรกฎาคม 2004, 00:03

ว่าแล้วเชียวว่าพี่กรต้องมีแล้ว เอาไว้ถ้าเจออันใหม่จะหามาฝากอีกครับ

gon · #10 15 กรกฎาคม 2004, 14:40

เกมจบแล้วครับ. สำหรับวันเวลาที่ค่อนข้างยาวนานในช่วงเดือนที่ผ่านมา. ไม่ว่าจะเป็นเอกลักษณ์ของรากที่ 3 ของ cos(p/9) + .. กับ เอกลักษณ์ของรากที่ 3 ของ cos (p/7)

รามานุจันได้เคยตั้งคำถามทั้งสองไว้ใน Journal ของอินเดียแล้ว ตอนแรกผมเจอเพียงเอกลักษณ์ของรากที่ 3 ของ cos(p/9) + .. กับ ทฤษฎีบทที่ทำให้ได้ เอกลักษณ์นี้ ใน The Lost Note book ของเขา. ถึงแม้ว่า ท.บ.ที่ผมใช้พิสูจน์เอกลักษณ์ของรากที่ 3 ของ cos(p/7) จะไม่ได้เขียนบันทึกไว้ใน The Lost Note Book ก็ตาม แต่การที่รามานุจันตั้งคำถามแบบนี้ได้ ก็หมายความว่า เขารู้อยู่แล้ว ตั้งแต่เกือบ 100 ปีก่อน.

ขอขอบคุณ Prof. ทั้ง 4 ที่ผมเมล์ไปถามและยินดีที่ตอบกลับมา ไว้ที่นี้อีกครั้งด้วย
Prof. Thomas J. Osler , (osler@rowan.edu) Rowan University ที่คุยตอบกลับไปกลับมากว่า 10 ฉบับ. อีกทั้งยังให้ยืมชื่อไปอ้างเพื่อนำไปคุยกับ Prof. ท่านอื่นด้วย
Prof George Andrews (andrews@math.psu.edu) ที่ส่งข้อความไปบอก Prof. Bruce Berndt อีกทาง
Prof Richard Askey (askey@math.wisc.edu) ที่ส่งข้อความไปบอก Prof. Bruce Berndt อีกทางเช่นกัน และ
Prof Bruce Berndt (berndt@math.uiuc.edu) Univ. of Illinois. ที่ช่วยคลายความโง่เขลาและให้กำลังใจแถมมาด้วย

Note. Prof. แต่ละคนถ่อมเองตัวเป็นบ้าเลย ขนาดตัวเองเชี่ยวชาญระดับไหน ยังบอกว่าตัวเองไม่เชี่ยวชาญ (I "m not expert) กันทุกคน อย่าง Prof. Thomas J. Osler ไม่รู้แกล้งหรือเปล่า ถามผมว่าสมการ 8x³ + 4x² - 4x - 1 = 0 รู้ได้ไงว่ารากของสมการคือ cos(2p/7, cos(4p/7), cos(6p/7) บอกให้ผม derive มาให้ดูหน่อย พอผมทำไปให้ดูเสร็จ กลับสรุปรูปแบบทั่วไปมาให้ผมดูอีก (งงไหมล่ะ) สมควรเอาอย่างครับ.

gon · #11 11 สิงหาคม 2004, 17:33

คิดไปคิดมา ได้มาอีกเพียบเลย ชักไม่ตื่นเต้นแล้ว รู้สึกว่ามันจะ unlimited ซะด้วย สำหรับเอกลักษณ์ของรากที่ 3 จำนวน 3 ตัว. ให้ดูเล่น ๆ ครับ.

gon · #12 11 สิงหาคม 2004, 17:35

จะเห็นได้ว่า มันไม่สวยเรียบง่ายเหมือนกับแบบแรก ๆ เท่านั้นเอง.เรื่องน่าสนใจเกี่ยวกับเอกลักษณ์ของรากที่ 3 ที่มีตอนนี้ก็คงจะเป็น จะถอดรากที่สามที่มันบวกกันอยู่ ออกมาเป็นเทอมง่าย ๆ ได้ไหม.

nooonuii · #13 12 สิงหาคม 2004, 09:06

ลอง simplify โดยใช้โปรแกรมดูครับ แต่ไม่แน่ใจว่าคำตอบจะสวยหรือเปล่า อยากเห็นคำตอบออกมาหน้าตาน่ารักเหมือนกันครับ

gon · #14 13 สิงหาคม 2004, 14:55

ถ้า simplify หลุดได้ก็ดีครับ. จะได้มีกำลังใจไปนั่ง Solve ต่อ ที่จริงพี่ก็เคยสั่งลอง Mathematica ให้มันลองดูอยู่ราว 5 นาที แต่พี่ใจร้อนครับ. ไม่หลุดสักทีก็เลิกเลย นี่ก็เอามาให้ดูอีกอัน จะเห็นได้ว่าเขียนออกมาได้เรื่อย ๆ เพราะที่จริงคิดเป็นทฤษฎีทั่วไปออกมาแล้วครับ. เดี๋ยวว่าจะเขียนเป็นบทความออกมาดู ไม่รู้จะมีคนอ่านหรือเปล่า

#15 13 สิงหาคม 2004, 17:11

อ่านซิครับ ไม่อ่านได้งไหละครับ

หัวข้อคล้ายคลึงกัน
หัวข้อ	ผู้ตั้งหัวข้อ	ห้อง	คำตอบ	ข้อความล่าสุด
An Euler's identity	<Pol>	พีชคณิต	0	21 กรกฎาคม 2001 10:34