數一B 88018100 翁婉珣
1.廣義積分 是否可積?
> Int((3+sin(x))/x,x=4..infinity);
> value(%);
>
ans:NO
2.求 x2ex2 以 x=0 為參考點的六階泰勒多項式。
> taylor(x^2*exp(x^2),x=0,6);
3. 求以下導函數:
(1)
> Diff(Int((t^4/(4!))*exp(-t),t=0..x),x);
> value(%);
(2)
> Diff(Int((t^4/4!)*exp(-t),t=0..4*x),x);
> value(%);
4. 求以下反導函數:
(1)
> Int(x^2*sin(x),x);
> value(%);
(2)
> Int(exp(x)*sin(x),x);
> value(%);
5. 令 v 是氧分子在室溫下的速度,以 m/sec 為單位。在為數非常多的氧分子中,各個速度不同。麥斯威爾 (Maxwell) 認為這些分子的速度分佈函數 (distribution function) 是 , 其中 是波茲曼 (Boltzmann) 常數, T 是凱氏溫標 (所謂室溫就是 293 度 K ),kg 是氧分子的質量。試問
> S:=v->v^2*exp((-m*v^2)/(2*k*T));
> p:=v->a*v^2*exp((-m*v^2)/(2*k*T));
> k:=1.4*10^(-23):
> T:=293:
> m:=5*10^(-26):
> a:=1/int(S(v),v=0..infinity);
(2). 室溫中氧分子的平均速度大約是多少?寫出四位有效數字。
> Int(v*(a*v^2*exp((-m*v^2)/(2*k*T))),v=0..infinity);
>
> value(%);
(3). p(v) 的最大值大約是多少?寫出四位有效數字。
> H:=diff(p(v),v);
> solve(H=0,v);
> p(-405.0678956);
> p(0);
> p(405.0678956);
6. 定義珈瑪函數
(1) 求 (1)、
(2)、 (3)、
(4)、 (5),觀察到什麼現象?
(2) 試求 (1/2)、
(1/4)、 (0)?
> Gamma(x)=Int(t^(x-1)*exp(-t),t=0..infinity);
> Gamma:=x->Int(t^(x-1)*exp(-t),t=0..infinity);
(1)
> value(Gamma(1));
> value(Gamma(2));
> value(Gamma(3));
> value(Gamma(4));
> value(Gamma(5));
ans:
> Gamma(x+1)=x*Gamma(x);
(2)
> value(Gamma(1/2));
> value(Gamma(1/4));
> value(Gamma(0));
7. 畫出 y=ex/4+0.5*sin(x), x 介於 0 與 2 之間,繞 y 軸旋轉 360 度所形成的曲面。 (提示:我不知道有沒有更好的方法,不過我會去看看圓柱坐標的繪圖指令。)
> with(plots);
>
cylinderplot(exp(1)^(x/4)+0.5*sin(x),theta=0..2*Pi,x=0..2*Pi);
8. 假設某加油站在地下橫置一支直徑 2 公尺,高 3 公尺的油桶 (實際的油桶可能會大得多)。油桶之頂側 (top side) 距加油嘴 (nozzle) 1 公尺,如下圖所示。
> (1^2*pi*3)*10^3 ;
(2)
high=1+h
> R:=Int(sqrt(1-x^2),x=0..sqrt(1-h^2))-h*sqrt(1-h^2);
> solve(R=5,h);
(3)
>
int((1+h)*6*r*((1^2-(1-h)^2)^(1/2)),h=0..2);