• Area
  • Properties of Integrals
  • Fundamental Theorem of Calculus

# Area

曲線下的面積 (只考慮[a,b][a,b] 上的非負函數)

Let S={(x,y)∣a≤x≤b,0≤y≤f(x)}.S = \{(x,y) | a \le x \le b , 0 \le y \le f(x) \}.
Divide [a,b][a,b] into nn subintervals:
[x0,x1],[x1,x2]...[xn−1,xn],wherex0=a,xn=b[x_0,x_1],[x_1,x_2]...[x_{n-1},x_n], \; where \; x_0 =a, x_n = b
Use a rectangle to approximate each SiS_i.Define that:

  • Un=f(x1)Δx1+f(x2)Δx2+...+f(xn)Δxn=∑i=1nf(xi)ΔxiU_n = f(x_1)\Delta x_1+f(x_2)\Delta x_2+...+f(x_n)\Delta x_n = \sum^n_{i=1}f(x_i)\Delta x_i
  • Rn=f(x0)Δx1+f(x1)Δx2+...+f(xn−1)Δxn=∑i=1nf(xi−1)ΔxiR_n = f(x_0)\Delta x_1+f(x_1)\Delta x_2+...+f(x_{n-1})\Delta x_n = \sum^n_{i=1}f(x_{i-1})\Delta x_i
    whereΔxn=xi−xi−1∀i∈{1,2,...,n}.where \Delta x_n = x_i - x_{i-1} \forall i \in \{ 1,2,...,n \}.

# Define

If f:[a,b]→R+∪{0}f: [a,b] \to \mathbb{R^+} \cup \{ 0 \} is continuous,
then A=lim⁡n→∞Rn=lim⁡n→∞∑i=1nf(xi)ΔxiA = \lim_{n \to \infty}R_n = \lim_{n \to \infty} \sum^n_{i=1}f(x_i)\Delta x_i
then A=lim⁡n→∞Ln=lim⁡n→∞∑i=1nf(xi−1)ΔxiA = \lim_{n \to \infty}L_n = \lim_{n \to \infty} \sum^n_{i=1}f(x_{i-1})\Delta x_i

# Partition

P={a=x0,x1,x2,...,xn=b}P = \{ a = x_0, x_1,x_2,...,x_n =b \} is a partition of [a,b][a,b].

  1. 等分: Δxi=b−an\Delta x_i = \frac{b-a}{n}.
  2. 不等分: Ex.x0=a,x1=ar,...,xn=arn=b,r=banab≥0,a≠0x_0 =a,x_1=ar,...,x_n = ar^n = b, r = \sqrt[n]{\frac{b}{a} } ab \ge 0 , a \neq 0
    P={a,ar,ar2,...,arn},Δxi=xi−xi−1=ari−ari−1.P = \{ a,ar,ar^2,...,ar^n \},\Delta x_i =x_i - x_{i-1} = ar^i - ar^{i-1}.

# Example

P={0=x0,x1=12n,x2=12n−1,...,xn−1=14,xn=12}.P = \{0 = x_0,x_1 = \frac{1}{2^n}, x_2 = \frac{1}{2^{n-1}},...,x_{n-1} = \frac{1}{4}, x_n = \frac{1}{2} \}.
Δxi=xi−xi−1=12n+1−x−12n−x,∥P∥=12−14=14.\Delta x_i = x_i - x_{i-1} = \frac{1}{2^{n+1-x}} - \frac{1}{2^{n-x}}, \left \| P \right \| = \frac{1}{2} - \frac{1}{4} = \frac{1}{4}.

Remark
  1. We want to find a partition of [a,b][a,b] such that ∥P∥→0\left \| P \right \| \to 0 implies n→∞n \to \infty.
  2. If n→∞n \to \infty,then ∥P∥→0\left \| P \right \| \to 0 不一定成立

# Example

f(x)=x2f(x) = x^2 on [a,b][a,b]. Find AA.

  1. Taking a partition of P={a=x0,x1,...,xn=b}P = \{ a = x_0, x_1,...,x_n=b \},where xi=a+b−ani,Δx=b−anx_i = a + \frac{b-a}{n}i,\Delta x = \frac{b-a}{n}.
  2. Ln=f(x0)Δx+f(x1)Δx+...+f(xn−1)Δx=b−an∑i=0n−1f(xi)=b−an∑i=0n−1{a2+2a(b−a)ni+(b−a)2n2i2}=a2(b−a)+2a(b−a)2n2⋅n(n−1)2+(b−a)3n3⋅n(n−1)(2n−1)6=a2(b−a)+a(b−a)2(1−1n)+(b−a)36(1−1n)(2−1n)L_n = f(x_0)\Delta x+f(x_1)\Delta x+...+f(x_{n-1})\Delta x = \frac{b-a}{n} \sum^{n-1}_{i=0}f(x_{i}) \\= \frac{b-a}{n} \sum^{n-1}_{i=0} \{ a^2 + 2\frac{a(b-a)}{n}i + \frac{(b-a)^2}{n^2}i^2 \} = a^2(b-a) + \frac{2a(b-a)^2}{n^2} \cdot \frac{n(n-1)}{2} + \frac{(b-a)^3}{n^3}\cdot \frac{n(n-1)(2n-1)}{6} \\= a^2(b-a) + a(b-a)^2(1-\frac{1}{n})+\frac{(b-a)^3}{6}(1-\frac{1}{n})(2-\frac{1}{n})
  3. A=lim⁡n→∞Ln=a2(b−a)+a(b−a)2+13(b−a)3A = \lim_{n \to \infty}L_n = a^2(b-a)+a(b-a)^2 + \frac{1}{3}(b-a)^3

# Define

UnU_n,上和: Taking xi∗∈[xi−1,xi]x_i* \in [x_{i-1},x_i] such that ff has a maximum f(xi∗)on[xi−1,xi]f(x_i*) on [x_{i-1},x_i],then ∑i=1nf(xi∗)Δxi\sum_{i=1}^{n}f(x_i*)\Delta x_i is called the upper sum of ff on [a,b][a,b] with a partition P={x0=a,x1,...,xn=b}P = \{ x_0 =a,x_1,...,x_n =b \}.

LnL_n,下和: Taking xi∗∈[xi−1,xi]x_i* \in [x_{i-1},x_i] such that ff has a minimum f(xi∗)on[xi−1,xi]f(x_i*) on [x_{i-1},x_i],then ∑i=1nf(xi∗)Δxi\sum_{i=1}^{n}f(x_i*)\Delta x_i is called the lower sum of ff on [a,b][a,b] with a partition P={x0=a,x1,...,xn=b}P = \{ x_0 =a,x_1,...,x_n =b \}.

# The Definite Integral

Divide [a,b][a,b] into nn subintervals (子區間) of each width Δx=b−an\Delta x = \frac{b-a}{n}:
P={x0=a,x1,...,xn=b}P = \{ x_0 =a,x_1,...,x_n = b \}.
Let xi∗x_i* be a sample point of subinterval [xi−1,xi][x_{i-1},x_i] for all i=1,2,...,ni = 1,2,...,n
The definite integral of ff from aa to bb is defined by ∫abf(x)dx=lim⁡n→∞∑i=1nf(xi∗)Δx\int^b_af(x)dx = \lim_{n\to \infty}\sum_{i=1}^nf(x_i*)\Delta x if the limit exists.
Moreover, we say that ff is integrable on [a,b][a,b].

# Example

∫0π2sinxdx=?\int^{\frac{\pi}{2}}_0 sinxdx = ?

Taking xi=0+π2−0n⋅i,i={1,2,3,...,n},Δx=π2nx_i = 0+\frac{\frac{\pi}{2}-0}{n}\cdot i, \; i = \{ 1,2,3,...,n\}, \; \Delta x = \frac{\pi}{2n}.

Let the right sum be ∑i=1nf(xi)Δx=∑i=1nπ2nsin(π2ni)=Rn\sum_{i=1}^nf(x_i)\Delta x = \sum_{i=1}^n \frac{\pi}{2n}sin(\frac{\pi}{2n}i) = R_n
∵2sin(π4n)sin(π2ni)=cos(π4n−π2ni)−cos(π4n+π2ni)=cos(π(1−2i)4n)−cos(π(1+2i)4n)\because 2sin(\frac{\pi}{4n})sin(\frac{\pi}{2n}i) = cos (\frac{\pi}{4n} - \frac{\pi}{2n}i) - cos (\frac{\pi}{4n} + \frac{\pi}{2n}i) = cos(\frac{\pi(1-2i)}{4n}) - cos(\frac{\pi(1+2i)}{4n})
⇒sin(π2ni)=12sin(π4n)(cos(π(1−2i)4n)−cos(π(1+2i)4n))\Rightarrow sin(\frac{\pi}{2n}i) = \frac{1}{2sin(\frac{\pi}{4n})}(cos(\frac{\pi(1-2i)}{4n}) - cos(\frac{\pi(1+2i)}{4n}))
∴Rn=12sin(π4n)⋅π2n∑i=1n(cos(π(1−2i)4n)−cos(π(1+2i)4n))=π4nsin(π4n){(cosπ4n−cos3π4n)+(cos3π4n−cos5π4n)+...+(cos(2n−1)π4n−cos(2n+1)π4n)}=π4nsin(π4n)(cos(π4n)−cos(π4(2+1n))).\therefore R_n = \frac{1}{2sin(\frac{\pi}{4n})} \cdot \frac{\pi}{2n} \sum^n_{i=1}(cos(\frac{\pi(1-2i)}{4n}) - cos(\frac{\pi(1+2i)}{4n})) = \frac{\frac{\pi}{4n}}{sin(\frac{\pi}{4n})} \{ (cos\frac{\pi}{4n} - cos\frac{3\pi}{4n})+(cos\frac{3\pi}{4n} - cos\frac{5\pi}{4n}) +...+(cos\frac{(2n-1)\pi}{4n} - cos\frac{(2n+1)\pi}{4n}) \} = \frac{\frac{\pi}{4n}}{sin(\frac{\pi}{4n})} (cos(\frac{\pi}{4n}) - cos(\frac{\pi}{4}(2+\frac{1}{n})) ).
⇒lim⁡n→∞Rn={1−0}=1\Rightarrow \lim_{n \to \infty}R_n = \{ 1-0 \} =1.

# Example

  1. lim⁡n→∞15+25+...+n5n6\lim_{n \to \infty} \frac{1^5+2^5+...+n^5}{n^6}.
  2. lim⁡n→∞1n{1+1n+1+2n+...+1+nn}\lim_{n\to \infty}\frac{1}{n} \{ \sqrt{1 + \frac{1}{n}} + \sqrt{1 + \frac{2}{n}}+...+\sqrt{1 + \frac{n}{n}} \}.
  3. lim⁡n→∞{11+n+12+n+...+12n}\lim_{n \to \infty} \{ \frac{1}{1+n} + \frac{1}{2+n} + ... + \frac{1}{2n} \}.

lim⁡n→∞15+25+...+n5n6=lim⁡n→∞1n{(1n)5+(2n)5+(3n)5+...+(nn)5}=lim⁡n→∞1n∑i=1n(in)5=lim⁡n→∞∑i=1n(in)51n=lim⁡n→∞∑i=1nf(xi)Δx=∫01x5dx=16x6∣x=0x=1=16\lim_{n \to \infty} \frac{1^5+2^5+...+n^5}{n^6} = \lim_{n \to \infty} \frac{1}{n} \{ (\frac{1}{n})^5 + (\frac{2}{n})^5 +(\frac{3}{n})^5 + ... + (\frac{n}{n})^5 \} \\= \lim_{n \to \infty} \frac{1}{n}\sum^n_{i=1}(\frac{i}{n})^5 = \lim_{n \to \infty} \sum^n_{i=1}(\frac{i}{n})^5\frac{1}{n} = \lim_{n \to \infty} \sum^n_{i=1}f(x_i)\Delta x \\= \int^1_0 x^5 dx = \frac{1}{6}x^6|^{x=1}_{x=0} = \frac{1}{6}.

lim⁡n→∞1n{1+1n+1+2n+...+1+nn}=lim⁡n→∞∑i=1n(1+in⋅1n)\lim_{n\to \infty}\frac{1}{n} \{ \sqrt{1 + \frac{1}{n}} + \sqrt{1 + \frac{2}{n}}+...+\sqrt{1 + \frac{n}{n}} \} = \lim_{n\to \infty} \sum^n_{i=1}(\sqrt{1+\frac{i}{n}} \cdot \frac{1}{n})

Let Δx=1−0n,xi=0+in,f(x)=1+x\Delta x = \frac{1-0}{n}, x_i = 0+\frac{i}{n} , f(x) = \sqrt{1+x}
lim⁡n→∞∑i=1n(1+in⋅1n)=∫011+xdx=23(1+x)32∣01=342−23\lim_{n\to \infty} \sum^n_{i=1}(\sqrt{1+\frac{i}{n}} \cdot \frac{1}{n}) = \int^1_0\sqrt{1+x}dx = \frac{2}{3}(1+x)^{\frac{3}{2}}|^1_0 = \frac{3}{4}\sqrt{2} - \frac{2}{3}.

Let Δx=2−1n,xi=1+in,f(x)=x\Delta x = \frac{2-1}{n}, x_i = 1+\frac{i}{n} , f(x) = \sqrt{x}
lim⁡n→∞∑i=1n(1+in⋅1n)=∫12xdx=23(x)32∣12=342−23\lim_{n\to \infty} \sum^n_{i=1}(\sqrt{1+\frac{i}{n}} \cdot \frac{1}{n}) = \int^2_1\sqrt{x}dx = \frac{2}{3}(x)^{\frac{3}{2}}|^2_1 = \frac{3}{4}\sqrt{2} - \frac{2}{3}.

lim⁡n→∞{11+n+12+n+...+12n}=lim⁡n→∞1n{11n+1+12n+1+...+1nn+1}=lim⁡n→∞∑i=1n(1n⋅11+2n)\lim_{n \to \infty} \{ \frac{1}{1+n} + \frac{1}{2+n} + ... + \frac{1}{2n} \} = \lim_{n \to \infty} \frac{1}{n} \{ \frac{1}{\frac{1}{n}+1} + \frac{1}{\frac{2}{n}+1} + ... +\frac{1}{\frac{n}{n}+1} \} = \lim_{n \to \infty}\sum^n_{i=1}(\frac{1}{n} \cdot \frac{1}{1+\frac{2}{n}})
Let Δx=1−0n,xi=in,f(x)=11+x\Delta x = \frac{1-0}{n} , x_i = \frac{i}{n}, f(x) = \frac{1}{1+x}.
∫0111+xdx=ln∣1+x∣∣x=0x=1=ln2\int^1_0 \frac{1}{1+x} dx = ln|1+x| |^{x=1}_{x=0} = ln2

# Example

∫abxqdx,q≠1,q∈Q,0<a<b\int^b_a x^q dx, q \neq 1, q \in \mathbb{Q}, 0<a<b

Let r=(ba)1nr = (\frac{b}{a})^\frac{1}{n} and xi=ari,i=0,1,...,n,Δxi=ari−1(r−1)x_i = ar^i, i = 0,1,...,n, \Delta x_i = ar^{i-1}(r-1).
P=x0=a,x1,x2,...xn=bP = {x_0 = a,x_1,x_2,...x_n=b}.
Rn=∑i=1nf(xi−1)Δxi=∑i=1n(xi−1)q⋅ari−1(r−1)=∑i=1n(ari−1)q⋅ari−1(r−1)=aq+1(r−1)∑i=1nrq(i−1)ri−1=aq+1(r−1)∑i=1nr(q+1)(i−1)=aq+1(r−1)1⋅(r(q+1)n−1)rq+1−1=r−1rq+1−1aq+1((ba)q+1−1)=r−1rq+1−1(bq+1−aq+1)R_n = \sum^n_{i=1}f(x_{i-1})\Delta x_i = \sum^n_{i=1}(x_{i-1})^q \cdot ar^{i-1}(r-1) = \sum^n_{i=1}(ar^{i-1})^q \cdot ar^{i-1}(r-1) \\= a^{q+1}(r-1) \sum^n_{i=1}r^{q(i-1)}r^{i-1} = a^{q+1}(r-1) \sum^n_{i=1}r^{(q+1)(i-1)} = a^{q+1}(r-1)\frac{1\cdot (r^{(q+1)n}-1)}{r^{q+1}-1} \\= \frac{r-1}{r^{q+1}-1}a^{q+1}((\frac{b}{a})^{q+1}-1) = \frac{r-1}{r^{q+1}-1} (b^{q+1} - a^{q+1})
∵ddr∣q=1rq+1=q+1\because \frac{d}{dr}|_{q=1}r^{q+1} = q+1
∴lim⁡r→1rq+1−1r−1=ddr∣q=1rq+1=q+1\therefore \lim_{r \to 1} \frac{r^{q+1}-1}{r-1} = \frac{d}{dr}|_{q=1}r^{q+1} = q+1

⇒lim⁡n→∞Rn=lim⁡n→∞∑i=1nf(xi−1)Δxi=lim⁡n→∞r−1rq+1−1(bq+1−aq+1)=lim⁡r→1r−1rq+1−1(bq+1−aq+1)=1q+1(bq+1−aq+1)=∫abxqdx\Rightarrow \lim_{n \to \infty}R_n = \lim_{n \to \infty}\sum^n_{i=1}f(x_{i-1})\Delta x_i = \lim_{n \to \infty}\frac{r-1}{r^{q+1}-1} (b^{q+1} - a^{q+1}) = \lim_{r \to 1}\frac{r-1}{r^{q+1}-1} (b^{q+1} - a^{q+1}) = \frac{1}{q+1}(b^{q+1}-a^{q+1}) = \int^b_a x^q dx.

# Properties of Integrals

  1. ∫abf(x)dx\int^b_af(x)dx is a real number.
  2. ∫abf(x)dx=∫abf(y)dy=∫abf(t)dt\int^b_af(x)dx = \int^b_af(y)dy = \int^b_af(t)dt
  3. ∫aaf(x)dx=0\int^a_af(x)dx = 0
  4. ∫abf(x)dx=−∫baf(x)dx\int^b_af(x)dx = -\int^a_bf(x)dx
  5. ∫abf(x)dx=∫acf(x)dx+∫cbf(x)dx\int^b_af(x)dx = \int^c_af(x)dx + \int^b_cf(x)dx (c∈Rc \in \mathbb{R})
  6. ∫abkdx=k(b−a)\int^b_a k dx = k(b-a)
  7. ∫ab(kf±hg)dx=k∫abf(x)dx±h∫abg(x)dx\int^b_a(kf \pm hg)dx = k\int^b_af(x)dx \pm h\int^b_ag(x)dx
  8. ∫abf(x)dx≤∣∫abf(x)dx∣≤∫ab∣f(x)∣dx\int^b_af(x)dx \le |\int^b_af(x)dx| \le \int^b_a|f(x)|dx
  9. If ff is odd, then ∫−aaf(x)dx=0\int^a_{-a}f(x)dx = 0
    If ff is even, then ∫−aaf(x)dx=2∫0af(x)dx\int^a_{-a}f(x)dx = 2\int^a_0f(x)dx

# Fundamental Theorem of Calculus

# Thm

If ff is continuous on [a,b][a,b], then g(x)=∫axf(t)dt,a≤x≤bg(x)= \int^x_a f(t)dt ,a\le x \le b is also continuous on [a,b][a,b] and differentiable on (a,b)(a,b) and g′(x)=f(x)g'(x) = f(x).
Moreover, ddx∫axf(t)dt=f(x)\frac{d}{dx}\int^x_af(t)dt = f(x).

# Example

  1. ddx∫3x(t5sin(2t+5))dt\frac{d}{dx}\int^x_3 (t^5sin(2t+5))dt
  2. ddx∫x1(u2ln∣u∣+tanu)du\frac{d}{dx}\int^1_x (u^2ln|u|+tanu)du
  3. Find the derivative of ∫2xxsinydy\int^x_{\sqrt{2}}xsinydy

ddx∫3x(t5sin(2t+5))dt=x5sin(2x+5)\frac{d}{dx}\int^x_3 (t^5sin(2t+5))dt = x^5sin(2x+5)

ddx∫x1(u2ln∣u∣+tanu)du=−ddx∫1x(u2ln∣u∣+tanu)du=−(x2ln∣x∣+tanx)\frac{d}{dx}\int^1_x (u^2ln|u|+tanu)du = -\frac{d}{dx} \int^x_1(u^2ln|u|+tanu)du = -(x^2ln|x|+tanx)

∫2xxsinydy=x∫2xsinydy\int^x_{\sqrt{2}}xsinydy = x\int^x_{\sqrt{2}}sinydy
⇒ddx∫2xxsinydy=ddxx∫2xsinydy=∫2xsinydy+xsinx=−cosx+cos2+xsinx\Rightarrow \frac{d}{dx}\int^x_{\sqrt{2}}xsinydy = \frac{d}{dx}x\int^x_{\sqrt{2}}sinydy = \int^x_{\sqrt{2}}sinydy+xsinx = -cosx+cos\sqrt{2}+xsinx

# Thm

If ff is continuous on [a,b][a,b],∫abf(x)dx=F(b)−F(a)\int^b_af(x)dx = F(b) - F(a),where F(x)F(x) is an antiderivative of ff.

# Example

∫f′(x)g(x)dx=f(x)g(x)−∫f(x)g′(x)dx\int f'(x)g(x)dx = f(x)g(x) - \int f(x)g'(x)dx
⇒∫abf′(x)g(x)dx=f(x)g(x)∣ab−∫abf(x)g′(x)dx=f(b)g(b)−f(a)g(a)−∫abf(x)g′(x)dx\Rightarrow \int^b_a f'(x)g(x)dx = f(x)g(x) |^b_a - \int^b_a f(x)g'(x)dx = f(b)g(b) -f(a)g(a) - \int^b_a f(x)g'(x)dx.

# Application

  1. ddx∫ag(x)f(t)dt=f(g(x))⋅g′(x)\frac{d}{dx}\int^{g(x)}_a f(t)dt = f(g(x)) \cdot g'(x)
  2. ddx∫h(x)g(x)f(t)dt=f(g(x))⋅g′(x)−f(h(x))⋅h′(x)\frac{d}{dx}\int^{g(x)}_{h(x)} f(t)dt = f(g(x)) \cdot g'(x) - f(h(x))\cdot h'(x)
  1. ddx∫ax2sintdt\frac{d}{dx}\int^{x^2}_a sint dt
  2. ddx∫sinxx2ln∣t2+2t∣dt\frac{d}{dx}\int^{x^2}_{sinx} ln|t^2+2t|dt
  3. lim⁡x→0∫0x21x6sint2dt\lim_{x\to 0}\int^{x^2}_0 \frac{1}{x^6}sint^2dt

ddx∫ax2sintdt=sinx2⋅ddxx2=2xsinx2\frac{d}{dx}\int^{x^2}_a sint dt = sinx^2 \cdot \frac{d}{dx}x^2 = 2xsinx^2

ddx∫sinxx2ln∣t2+2t∣dt=ln∣x2+2(x2)2∣⋅ddxx2−ln∣sinx+2sin2x∣⋅ddxsinx=2xln∣x2+2x4∣−cosxln∣sinx+2sin2x∣=ln∣x2+2x4∣2x−ln∣sinx+2sin2x∣cosx=ln∣ln∣x2+2x4∣2xln∣sinx+2sin2x∣cosx∣\frac{d}{dx}\int^{x^2}_{sinx} ln|t^2+2t|dt = ln|x^2+2(x^2)^2| \cdot \frac{d}{dx} x^2 - ln|sinx + 2sin^2x| \cdot \frac{d}{dx} sinx = 2xln|x^2+2x^4| - cosxln|sinx + 2sin^2x| = ln|x^2+2x^4|^{2x} - ln|sinx + 2sin^2x|^{cosx} = ln|\frac{ln|x^2+2x^4|^{2x}}{ln|sinx + 2sin^2x|^{cosx}}|

lim⁡x→0∫0x21x6sint2dt=lim⁡x→0∫0x2sint2dtx6=lim⁡x→0sinx4⋅2x6x5=lim⁡x→0cosx4⋅4x312x3=13\lim_{x\to 0}\int^{x^2}_0 \frac{1}{x^6}sint^2dt = \lim_{x \to 0} \frac{\int^{x^2}_0 sint^2 dt}{x^6} = \lim_{x \to 0}\frac{sinx^4 \cdot 2x}{6x^5} = \lim_{x \to 0} \frac{cosx^4 \cdot 4x^3}{12x^3} = \frac{1}{3}.

# Thm

If ff is continuous on [a,b][a,b], then

  1. ddx∫axf(t)dt=f(x)\frac{d}{dx}\int^x_af(t) dt = f(x)
  2. ∫abf(x)dx=F(b)−F(a)\int^b_af(x)dx = F(b) - F(a) if F′(x)=f(x)F'(x) = f(x)

# Reference

  • 蘇承芳老師 - 微積分甲(一)109 學年度 - Calculus (I) Academic Year 109
更新於 閱讀次數 次

用實際行動犒賞爆肝的我😀

Zrn Ye LinePay

LinePay