查看完整版本 : 數學謎題車輪陣 (二)

rhwlam 2018-11-12 02:34 PM

數學謎題車輪陣 (二)

[b]破陣英雄榜[/b]
第01題(#02): XMing (#03)
第02題(#12): XMing (#13)
第03題(#16): XMing (#24)
第04題(#34): XMing (#35)
第05題(#37): XMing (#38)
第06題(#42): XMing (#49)
第07題(#52): XMing (#57)
第08題(#53): XMing (#54)
第09題(#63): XMing (#67)
第10題(#64): XMing (#67)

[[i] 本帖最後由 rhwlam 於 2019-1-15 03:25 PM 編輯 [/i]]

rhwlam 2018-11-12 02:39 PM

[b]第01題[/b]
若f為一個單元函數, 而f在任何參數值也可微分, 另外:
f(x) + f(y) = f( (x + y) / (1 - xy) )
以上x和y皆為實數, 而xy不等於1.
請找出所有可以符合f要求的函數.

XMing 2018-11-12 03:48 PM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-12 02:39 PM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490424629&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
第01題
若f為一個單元函數, 而f在任何參數值也可微分, 另外:
f(x) + f(y) = f( (x + y) / (1 - xy) )
以上x和y皆為實數, 而xy不等於1.
請找出所有可以符合f要求的函數. ... [/quote]
when x=0, y=0, we get f(0)=0
Fix y, differentiate w.r.t x on both sides and use chain rule and quotient rule for differentiation , we get
f'(x)=f'((x+y)/(1-xy)) * [(1-xy)(1)-(x+y)(-y)]/[(1-xy)^2]
     =f'((x+y)/(1-xy)) * (1+y^2)/(1-xy)^2
Since y can be any real number except 1/x (since xy<>1)
Let y=-x, xy=-x^2<=0, so xy<>1, then we obtain the following differential  equation
f'(x)=f'(0)* (1+x^2)/(1+x^2)^2=f'(0)/(1+x^2)
If f'(0)=0, then f(x)=c where c is constant but f(0)=0 therefore f(x)=0
If f'(0)<>0, then solve the DE f'(x)=f'(0)/(1+x^2)
f(x)=f'(0)arctan(x)+f(0)=f'(0)arctan(x)

[[i] 本帖最後由 XMing 於 2018-11-12 10:52 PM 編輯 [/i]]

rhwlam 2018-11-12 10:05 PM

[quote]原帖由 [i]XMing[/i] 於 2018-11-12 15:48 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490427391&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
when x=0, y=0, we get f(0)=0
Fix y, differentiate w.r.t x on both sides and use chain rule and quotient rule for differentiation , we get
f'(x)=f'((x+y)/(1-xy)) * [(1-xy)(1)-(x+y)(-y)]/[(1-xy)^2]
... [/quote]
一下子就解答了. 厲害!
我是在想, 以這個推論來說, 因為y可以是任何實數, 我們也可以假設y = 1/x, 那如何處理呢?
另外, 若y不設成y = -x, 又如何推論下去呢?

XMing 2018-11-12 10:47 PM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-12 10:05 PM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490444222&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]

一下子就解答了. 厲害!
我是在想, 以這個推論來說, 因為y可以是任何實數, 我們也可以假設y = 1/x, 那如何處理呢?
另外, 若y不設成y = -x, 又如何推論下去呢? ... [/quote]

我搞錯咗, xy<>1, 不能設 y=1/x. 巳在#3修正.

XMing 2018-11-16 10:59 AM

[quote]原帖由 [i]XMing[/i] 於 2018-11-12 03:48 PM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490427391&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]

when x=0, y=0, we get f(0)=0
Fix y, differentiate w.r.t x on both sides and use chain rule and quotient rule for differentiation , we get
f'(x)=f'((x+y)/(1-xy)) * [(1-xy)(1)-(x+y)(-y)]/[(1-xy)^2]
... [/quote]

[font=新細明體][/font]其實好像有點不妥.
[font=新細明體][/font]對所有實數 x,y (xy<>1) 滿足f(x)+f(y)=f((x+y)/(1-xy))..... (1)
[font=新細明體][/font]必滿足微分方程 f(x)=f'(0)/(1+x^2) 及f(0)=0
[font=新細明體][/font]而它的解是 f(x)=f'(0)atan(x)
[font=新細明體][/font]但是反轉過來, 未必所有解都會使 (1) 成立.
[font=新細明體][/font]例如: 假如f'(0)=1, f(x)=atan(x)
[font=新細明體][/font]x=2, y=3
[font=新細明體][/font]f(2)=arctan(2)=1.107149
[font=新細明體][/font]f(3)=arctan(3)=1.249046
[font=新細明體][/font]f((2+3)/(1-2*3))=f(-1)=arctan(-1)=-0.7854
[font=新細明體][/font]好明顯,f(2)+f(3)<>f((2+3)/(1-2*3))
[font=新細明體][/font]所以, 只有f(x)=0 才可以令 (1) 成立.
[font=新細明體][/font]

[[i] 本帖最後由 XMing 於 2018-11-16 11:03 AM 編輯 [/i]]

rhwlam 2018-11-16 11:21 AM

[quote]原帖由 [i]XMing[/i] 於 2018-11-16 10:59 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490616983&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
其實好像有點不妥.
對所有實數 x,y (xy<>1) 滿足f(x)+f(y)=f((x+y)/(1-xy))..... (1)
必滿足微分方程 f(x)=f'(0)/(1+x^2) 及f(0)=0
而它的解是 f(x)=f'(0)atan(x)
但是反轉過來, 未必所有解都會使 (1) 成立.
例如: 假如f'(0)=1, f(x)=atan(x)
x=2, y=3
f( ... [/quote]
XMing兄的求真精神, 令人佩服!
另外, 可考慮為  f(x)=f'(0)atan(x) 多加一個條件: -1 < x < 1.

[[i] 本帖最後由 rhwlam 於 2018-11-16 11:42 AM 編輯 [/i]]

XMing 2018-11-16 12:09 PM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-16 11:21 AM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490617969&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]

XMing兄的求真精神, 令人佩服!
另外, 可考慮為  f(x)=f'(0)atan(x) 多加一個條件: -1 < x < 1. [/quote]

加了-1<x<1未必保證到 (x+y)/(1-xy) 在 (-1,1) 之內.
例如 x=y=0.9, (0.9+0.9)/(1-0.9*0.9)=9.47 在 (-1,1) 之外
那麽 f(x)+f(y)=f((x+y)/(1-xy)) 就不成立.

[[i] 本帖最後由 XMing 於 2018-11-16 12:10 PM 編輯 [/i]]

rhwlam 2018-11-16 12:23 PM

[quote]原帖由 [i]XMing[/i] 於 2018-11-16 12:09 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490620490&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
加了-1<x<1未必保證到 (x+y)/(1-xy) 在 (-1,1) 之內.
例如 x=y=0.9, (0.9+0.9)/(1-0.9*0.9)=9.47 在 (-1,1) 之外
那麽 f(x)+f(y)=f((x+y)/(1-xy)) 就不成立. [/quote]

atan(0.9) = 0.733
atan([(0.9+0.9)/(1 - 0.9*0.9)]) = 1.466

可能我是錯的, 但無論如何, 很享受跟XMing兄討論~

[[i] 本帖最後由 rhwlam 於 2018-11-22 07:16 PM 編輯 [/i]]

XMing 2018-11-16 12:37 PM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-16 12:23 PM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490621209&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]


atan(0.9) = 0.733
atan([(0.9+0.9)/(1 - 0.9*0.9)]) = 1.466

-1<x<1的條件是在1/(1+x^2)的integral出來的.
可能我是錯的, 但無論如何, 很享受跟XMing兄討論~ [/quote]

一個function有三部份.
(1) Domain (2) rule to assign each element of domain to range. (3) Range
既然所定義的 Domain 是 (-1,1) 那麽計算 f((0.9+0.9)/(1-0.9*0.9)) 是沒有定義的並非等於atan(9.3684) 因為 9.3684 並不在定義域(domain)範圍.

奇實, 你想過的我也有想過.

rhwlam 2018-11-16 12:50 PM

[quote]原帖由 [i]XMing[/i] 於 2018-11-16 12:37 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490621856&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
一個function有三部份.
(1) Domain (2) rule to assign each element of domain to range. (3) Range
既然所定義的 Domain 是 (-1,1) 那麽計算 f((0.9+0.9)/(1-0.9*0.9)) 是沒有定義的並非等於atan(9.3684) 因為 9.3684 並不在定義域(domain)範圍.
... [/quote]
對呀, 我也覺得奇怪, 所以我也說可能我是錯的.

rhwlam 2018-11-16 04:42 PM

[b]第02題[/b]
若[i]k[/i]是正整數, 請計算
[attach]9008773[/attach]
為一個只包括[i]k[/i]而不包括[i]x[/i]的函數.

[[i] 本帖最後由 rhwlam 於 2018-11-16 04:45 PM 編輯 [/i]]

XMing 2018-11-17 12:20 PM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-16 04:42 PM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490632515&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
第02題
若k是正整數, 請計算
9008773
為一個只包括k而不包括x的函數. [/quote]
(2^k-1)/(2^k)*(pi/2) ?

rhwlam 2018-11-17 12:35 PM

[quote]原帖由 [i]XMing[/i] 於 2018-11-17 12:20 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490667044&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
(2^k-1)/(2^k)*(pi/2) ? [/quote]
對的. :smile_o12:
若能提供推論更佳~

XMing 2018-11-18 11:01 AM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-17 12:35 PM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490667697&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]

對的. :smile_o12:
若能提供推論更佳~ [/quote]
Integral [0 inf] (sin x/x) dx = pi/2
Let u=kx. It can be easily prove that
Integral [0 inf] (sin kx/x) dx = pi/2
sin kx (cos x)^k=(sin kx cos x) (cos x)^(k-1)
=1/2*(sin(k+1)x+sin(k-1)x) (cos x)^(k-1)
=1/2*(sin(k+1)x cos x + sin(k-1)x cos x) (cos x)^(k-2)
=1/4*(sin(k+2)x+2sin(kx)+sin(k-1)x) (cos x)^(k-2)
=...
=1/(2^k)*(sin(2k)x+kC1sin(2(k-1)x)+kC2sin(2(k-2)x)+...+kC(k-1)sin2x)
Note that 1+kC1+kC2+...+kC(k-1)=2^k-1
So the answer is (2^k-1)/(2^k)*pi/2
QED

[[i] 本帖最後由 XMing 於 2018-11-18 11:04 AM 編輯 [/i]]

rhwlam 2018-11-18 09:31 PM

[b]第03題[/b]
若y^3 + y^2*x - x^2*y - x^3 = 1,
請計算:
[attach]9018112[/attach]
x和y皆為實數.

可考慮以計數機(但不是電腦)作答, 答案可以只提供前兩個有效數字(即誤差為<=5%).

[b]註:[/b] 這題難度應該較高, 推論亦可能較冗長.

[[i] 本帖最後由 rhwlam 於 2018-11-21 10:51 PM 編輯 [/i]]

XMing 2018-11-19 04:47 PM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-18 09:31 PM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490737489&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
第03題
若y^3 + y^2*x - x^2*y - x^3 = 1,
請計算:
9018112
x和y皆為實數.

可考慮以計數機(但不是電腦)作答, 答案可以只提供前兩個有效數字(即誤差為<=0.5%).

註: 這題難度應該較高, 推論亦可能較冗長. ... [/quote]
0.73-0.74 ?

rhwlam 2018-11-21 11:19 PM

[quote]原帖由 [i]XMing[/i] 於 2018-11-19 16:47 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490773553&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
0.73-0.74 ? [/quote]
這一題我沒有答案.
我沒有編程計算答案, 如果以計算機計算, 自己的答案是~1.7. 不知是否正確.

XMing 2018-11-21 11:31 PM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-21 11:19 PM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490898177&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]

這一題我沒有答案.
我沒有編程計算答案, 如果以計算機計算, 自己的答案是~1.7. 不知是否正確. [/quote]

我沒有編程, 但是我用Excel 計的. 未必乎合你的要求

[[i] 本帖最後由 XMing 於 2018-11-21 11:35 PM 編輯 [/i]]

rhwlam 2018-11-22 12:23 AM

大概計算方法是:
- 先解出y^3 + y^2*x + (-x^2)*y + (-x^3-1) = 0 中y的實數解. 當中要順便證明出當x>=0時y只有一個實數解.
- 我們現將積分轉為:
[attach]9030111[/attach][b][color=#ff0000](這個積分是錯的, 抱歉! 更改了的積分在帖#23)[/color][/b]
(如有高手如何處理以上積分, 請告知, 感謝! 那麼, 便不用以下愚笨的數值積分了...)
- 若考慮f(x) = 4/3[(x^3 + 27/16)^(1/3) - x], f(x)的收斂要比一些擁有有限積分值的函數(如1/(1+x^2))要快, 因此f(x)的積分會有一個有限的收斂值.
- 若以1/(1+x^2)作為參考, 其由x=0至無限的積分為arctan(x) (= 1.57), 若我們只由x=0至80作積分, 結果(1.558)只會有<1%的誤差.
- 因此, 我們可以考慮由x=0至80對f(x)作出數值積分, 第一次嘗試為[f(0)+f(80)]*80/2 (梯形面積), 第二次嘗試可以為[f(0)+f(40)]*40/2 + [f(40)+f(80)]*40/2 (兩個梯形面積總和).
(其實, 以下我是利excel代替了計算機)
- 直至第7次嘗試, 當我們將這個積分分拆64個梯形的總和時, 我們大概知道答案為1.7或1.8. 我們只要求誤差為<2.5%, 所以其實我們只需考慮積分至x~40.
- 第8次嘗試: 由x=0至x=40的區間, 我們將f(x)的積分計算128個梯形的總和, 我們就可得出答案為~1.7.

[[i] 本帖最後由 rhwlam 於 2018-11-22 11:36 AM 編輯 [/i]]

rhwlam 2018-11-22 12:25 AM

[quote]原帖由 [i]XMing[/i] 於 2018-11-21 23:31 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490898719&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
我沒有編程, 但是我用Excel 計的. 未必乎合你的要求 [/quote]
我也用了excel, 代勞數值積分的部份. 剛已寫上自己的計算方法, 不知對不對.

XMing 2018-11-22 01:04 AM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-22 12:23 AM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490900921&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
大概計算方法是:
- 先解出y^3 + y^2*x + (-x^2)*y + (-x^3-1) = 0 中y的實數解. 當中要順便證明出當x>=0時y只有一個實數解.
- 我們現將積分轉為:

(如有高手如何處理以上積分, 請告知, 感謝! 那麼, 便不用以下愚笨的數值積分了...)
- 若考慮f(x) = 4/3[(x^3 + 27/16)^(1/3) - x], f(x)的收斂要比一 ... [/quote]

當 x=0 代入 y^3+y^2x-yx^2-x^3=1. 得 y=1
積分函數 f(x)=y-x. 所以 f(0)=1
你的f(x)= 4/3[(x^3 + 27/16)^(1/3) - x]
f(0)=4/3(27/16)^(1/3)<>1

一步錯, 步步錯, 都不需要看下去....

[[i] 本帖最後由 XMing 於 2018-11-22 11:02 AM 編輯 [/i]]

rhwlam 2018-11-22 11:33 AM

[quote]原帖由 [i]XMing[/i] 於 2018-11-22 01:04 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490902110&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
一步錯, 步步錯, 都不需要看下去... [/quote]
謝謝提點.
我重新算了一次, 其積分可能會變成:
[attach]9031173[/attach]
那麼如之前帖#20的計算次序, 可得答案為~0.75.

XMing 2018-11-22 04:01 PM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-22 11:33 AM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490915556&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]

謝謝提點.
我重新算了一次, 其積分可能會變成:
9031173
那麼如之前帖#20的計算次序, 可得答案為~0.75. [/quote]
我計到
f(x)=(8/27x^3+1/2+sqrt(8/27x^3+1/4))^(1/3)+(8/27x^3+1/2-sqrt(8/27x^3+1/4))^(1/3)-4/3x
再用trapezoidal rule 去計算積分 由 0 至 40, delta x = 0.001
大約=0.743

如果要計算誤差, 是要計f(x)的second derivative, 似乎太煩了, 所以用一個很小的delta去計
另外過x>40的積分估計應該有小.

[[i] 本帖最後由 XMing 於 2018-11-22 04:08 PM 編輯 [/i]]

rhwlam 2018-11-22 05:18 PM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-22 11:33 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490915556&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
我重新算了一次, 其積分可能會變成:
9031173
那麼如之前帖#20的計算次序, 可得答案為~0.75. [/quote]
大概計算方法是:
- 先解出y^3 + y^2*x + (-x^2)*y + (-x^3-1) = 0 中y的實數解. 當中要順便證明出當x>=0時y只有一個實數解.
- 我們現將積分轉為:
[attach]9032168[/attach]
(如有高手如何處理以上積分, 請告知, 感謝! 那麼, 便不用以下愚笨的數值積分了...)
- 若考慮以上積分符号內的函數為f(x), 即f(x) = (...)^(1/3) - (...)^(1/3) - 4x/3, f(x)的收斂要比一些擁有有限積分值的函數(如1/(1+x^2))要快, 因此f(x)的積分會有一個有限的收斂值.
- 若以1/(1+x^2)作為參考, 其由x=0至無限的積分為arctan(x) (= 1.57), 若我們只由x=0至80作積分, 結果(1.558)只會比最後答案少<1%.
- 因此, 我們可以考慮由x=0至80對f(x)作出數值積分, 第一次嘗試為[f(0)+f(80)]*80/2 (梯形面積), 第二次嘗試可以為[f(0)+f(40)]*40/2 + [f(40)+f(80)]*40/2 (兩個梯形面積總和).
(其實, 以下我是利excel代替了計算機)
- 直至第10次嘗試, 當我們將這個積分分拆512個梯形的總和時, 我們大概知道答案為0.748861.
- 直至第11次嘗試, 當我們將這個積分分拆1024個梯形的總和時, 我們大概知道答案為0.747555.
- 我們可以知道, 至x=80積分會比0.747555 - (0.748861-0.747555) = 0.746249為大, 同時只會比x-->'無限'的答案少<1%.
- 所以答案為~0.75.

XMing 2018-11-22 05:38 PM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-22 05:18 PM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490930521&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]

大概計算方法是:
- 先解出y^3 + y^2*x + (-x^2)*y + (-x^3-1) = 0 中y的實數解. 當中要順便證明出當x>=0時y只有一個實數解.
- 我們現將積分轉為:
9032168
(如有高手如何處理以上積分, 請告知, 感謝! 那麼, 便不用以下愚笨的數值積分了...)
- 若考慮以上積分符号內的函數為f(x), 即f(x) = (...)^(1/3) - ... [/quote]

其實呢題我無mud興趣.....

[[i] 本帖最後由 XMing 於 2018-11-22 05:45 PM 編輯 [/i]]

rhwlam 2018-11-22 07:54 PM

[quote]原帖由 [i]XMing[/i] 於 2018-11-22 17:38 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490931482&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
其實呢題我無mud興趣..... [/quote]
雖說沒興趣, 但師兄一開始便答對了.
我就要靠師兄提點才能答對. :smile_35:

XMing 2018-11-22 11:29 PM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-22 05:18 PM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490930521&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]

大概計算方法是:
- 先解出y^3 + y^2*x + (-x^2)*y + (-x^3-1) = 0 中y的實數解. 當中要順便證明出當x>=0時y只有一個實數解.
- 我們現將積分轉為:
9032168
(如有高手如何處理以上積分, 請告知, 感謝! 那麼, 便不用以下愚笨的數值積分了...)
- 若考慮以上積分符号內的函數為f(x), 即f(x) = (...)^(1/3) - ... [/quote]

唔係好明你點計 f(x)
以下是我的計法
設 y=z-x/3
(z-x/3)^3+x(z-x/3)^2-x^2(z-x/3)-x^3-1=0
化簡後得到
z^3-4/3x^2z-16/27x^3-1=0
假設 z^3+pz+q=0
p=-4/3x^2,
q=-16/27x^3-1
D=q^2/4+p^3/27=8/27x^3+1/4
當 x>0, D>0 三次方程有一實根, 兩虛根
所以 z=(-q/2+sqrt(D))^(1/3)+(-q/2-sqrt(D))^(1/3)
=(8/27x^3+1/2+sqrt(8/27x^3+1/4))^(1/3)+(8/27x^3+1/2-sqrt(8/27x^3+1/4))^(1/3)
y=z-x/3=(8/27x^3+1/2+sqrt(8/27x^3+1/4))^(1/3)+(8/27x^3+1/2-sqrt(8/27x^3+1/4))^(1/3)-x/3
f(x)=y-x=(8/27x^3+1/2+sqrt(8/27x^3+1/4))^(1/3)+(8/27x^3+1/2-sqrt(8/27x^3+1/4))^(1/3)-4x/3

rhwlam 2018-11-23 12:33 AM

[quote]原帖由 [i]XMing[/i] 於 2018-11-22 23:29 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490945459&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]
唔係好明你點計 f(x)
以下是我的計法 [/quote]
我們的計法是一模一樣的.
這一題找得有點草率, 沒有先計算好就放上來. 之後發覺積分部分很麻煩, 所以自己嘗試把它計算了算了, 免得師兄傷神推出其答案數值.

XMing 2018-11-24 11:21 PM

[quote]原帖由 [i]rhwlam[/i] 於 2018-11-23 12:33 AM 發表 [url=https://www.discuss.com.hk/redirect.php?goto=findpost&pid=490947476&ptid=27845599][img]https://www.discuss.com.hk/images/common/back.gif[/img][/url]

我們的計法是一模一樣的.
這一題找得有點草率, 沒有先計算好就放上來. 之後發覺積分部分很麻煩, 所以自己嘗試把它計算了算了, 免得師兄傷神推出其答案數值. ... [/quote]Exact value 是 0.75 不需要用numerical method 用積分都可計到[attach]9040086[/attach][attach]9040087[/attach]

尾5嗰行第三條算式的upper limit 應該是sqrt(s)+1/2 而不是 sqrt(s)-1/2最尾嗰行 s-> infinity

[[i] 本帖最後由 XMing 於 2018-11-25 12:37 AM 編輯 [/i]]
頁: [1] 2 3
查看完整版本: 數學謎題車輪陣 (二)