Main derivation
S01 · definition
Before
f_x(x,y)
After
\lim_{h\to0,\ h\neq0}\frac{f(x+h,y)-f(x,y)}{h}Conditions: y\ \text{fixed}
Rule: mko_partial_def · Uses: mko_partial_def
Leaf basis: mko_partial_def + def_unfold
Unfold the x-partial by the MKO: y stays fixed, only the first coordinate changes by h (h≠0), and the difference quotient is taken as h→0. The y-partial fixes x and moves only the second coordinate (MKO paragraph 0). A partial derivative is a limit, not the quotient at one particular h.
S02 · definition
Before
c=(c_0,c_1,c_2,c_3,c_4,c_5)
After
f(x,y)=c_0+c_1x+c_2y+c_3x^2+c_4xy+c_5y^2
Conditions: none
Rule: mko_mvf_def · Uses: G1, mko_mvf_def
Leaf basis: mko_mvf_def + def_unfold
Read the coefficient vector as a polynomial in the fixed order. The derivation below works for arbitrary coefficients, and the instance of G2 is substituted at the end.
S03 · equivalent
Before
f(x+h,y)-f(x,y)
After
h(c_1+2c_3x+c_4y)+c_3h^2
Conditions: none
Rule: L_numerator_x · Uses: S02, L_numerator_x
Leaf basis: L_numerator_x + distributive
Lemma L_numerator_x (proved in S03.1, S03.2): the numerator is quadratic in h with zero constant term. N1 checks this polynomial identity on 3^9 grid points (degree at most 2 in each variable, which decides it).
S04 · equivalent
Before
\frac{h(c_1+2c_3x+c_4y)+c_3h^2}{h}After
q_x(h)=(c_1+2c_3x+c_4y)+c_3h\quad(h\neq0)
Conditions: h\neq0
Rule: mul_inverse_nonzero · Uses: S03, mul_inverse_nonzero
Leaf basis: mul_inverse_nonzero + frac_rules + power_laws
Since h≠0, we may multiply by 1/h: h·A/h=A and c_3h²/h=c_3h. So for every nonzero h, the quotient is the affine expression A+Bh in h, with A=c_1+2c_3x+c_4y and B=c_3. At h=0 the quotient is undefined, but the limit does not look at h=0.
S05 · theorem_application
Before
\lim_{h\to0}\big(A+Bh\big)After
A
Conditions: A=c_1+2c_3x+c_4y,\ B=c_3
Rule: L_affine_limit · Uses: S04, L_affine_limit
Leaf basis: L_affine_limit + func_limit_def
Lemma L_affine_limit (proved by ε-δ in S05.1, S05.2): the affine expression A+Bh has limit A as h→0. This turns 'just substitute h=0' into an argument that can be checked.
S06 · implies
Before
f_x(x,y)=\lim_{h\to0}q_x(h)After
f_x(x,y)=c_1+2c_3x+c_4y
Conditions: none
Rule: L_affine_limit · Uses: S01, S04, S05
Leaf basis: L_affine_limit + def_unfold
S04 says the quotient equals A+Bh for h≠0, and S05 says its limit is A. By the definition in S01, the x-partial exists and equals A=c_1+2c_3x+c_4y. This is an exact formula valid at every point (x,y) and for arbitrary coefficients.
S07 · implies
Before
\frac{f(x,y+h)-f(x,y)}{h}=(c_2+c_4x+2c_5y)+c_5hAfter
f_y(x,y)=c_2+c_4x+2c_5y
Conditions: h\neq0; x\ \text{fixed}
Rule: L_affine_limit · Uses: S01, S02, S05
Leaf basis: L_affine_limit + distributive + mul_inverse_nonzero
The y-partial fixes x and moves only y. Subtracting term by term in the same way, f(x,y+h)−f(x,y)=h(c_2+c_4x+2c_5y)+c_5h² (N2 is an identity); dividing by h≠0 gives the affine A'+B'h with A'=c_2+c_4x+2c_5y and B'=c_5. By L_affine_limit, f_y=A'.
S08 · equivalent
Before
c=(1,-2,3,1,2,-1)
After
f_x=-2+2x+2y,\quad f_y=3+2x-2y
Conditions: none
Rule: eq_replace · Uses: S06, S07, G2
Leaf basis: eq_replace + numeric_evaluation
Substitute c_1=−2, c_2=3, c_3=1, c_4=2, c_5=−1: f_x=−2+2·1·x+2y=−2+2x+2y; f_y=3+2x+2·(−1)·y=3+2x−2y.
S09 · equivalent
Before
(x,y)=(1,2)
After
f_x(1,2)=4,\ \ f_y(1,2)=1,\ \ \nabla f(1,2)=(4,1)
Conditions: none
Rule: numeric_evaluation · Uses: S08
Leaf basis: numeric_evaluation
f_x(1,2)=−2+2+4=4 (N3); f_y(1,2)=3+2−4=1 (N4). The gradient just collects the two partials in coordinate order: (4,1) (MKO paragraph 1).
S10 · implies
Before
\frac{f(1+h,2)-f(1,2)}{h}After
\frac{4h+h^2}{h}=4+h\ \xrightarrow{h\to0}\ 4Conditions: h\neq0
Rule: L_affine_limit · Uses: G2, S04, S09
Leaf basis: L_affine_limit + numeric_evaluation + distributive
Check the instance directly: f(1,2)=6 (N5) and f(1+h,2)−6=4h+h² (N6 is an identity in h), so the quotient is 4+h with limit 4, agreeing with S09. This is a computation independent of the general formula.
S11 · implies
Before
q_x(h)=4+h
After
q_x(1)=5,\ \ q_x(\tfrac1{10})=\tfrac{41}{10},\ \ q_x(h)\neq4\ \text{for every}\ h\neq0Conditions: h\neq0
Rule: order_basics · Uses: S10
Leaf basis: order_basics + numeric_evaluation
h=1 gives the quotient (f(2,2)−6)/1=5 (N7); h=1/10 gives 41/10 (N8). In fact q_x(h)=4+h differs from 4 for every nonzero h. So no single nonzero-h quotient is the partial derivative; the partial is the limit as h→0 (MKO paragraph 2: a finite difference is an approximation or a sample of the quotient).
S12 · implies
Before
\frac{f(1+h,2+h)-f(1,2)}{h}After
\frac{5h+2h^2}{h}=5+2h\ \xrightarrow{h\to0}\ 5\neq4=f_x(1,2)Conditions: h\neq0
Rule: L_affine_limit · Uses: G2, S09, L_affine_limit
Leaf basis: L_affine_limit + numeric_evaluation + distributive
Moving both coordinates at once: f(1+h,2+h)−6=5h+2h² (N9 is an identity), the quotient is 5+2h, and by L_affine_limit its limit is 5. This is the rate of change along the direction (1,1), which here equals f_x+f_y=4+1, not the x-partial 4. A partial derivative must keep the other coordinate fixed (the misconception named by the MKO).
S13 · implies
Before
\text{S06 to S12}After
f_x=c_1+2c_3x+c_4y,\ f_y=c_2+c_4x+2c_5y;\ \nabla f(1,2)=(4,1)
Conditions: none
Rule: and_intro_elim · Uses: S06, S07, S09, S11, S12
Leaf basis: and_intro_elim
Summary: the partial-derivative formulas of a quadratic follow from the difference-quotient limit; the instance has gradient (4,1) at (1,2); single-h quotients (5, 41/10) and the limit with both coordinates moving (5) are not the x-partial.
Substeps and lemmas
S03.1 · equivalent
Before
c_0,\ c_1x,\ c_2y,\ c_3x^2,\ c_4xy,\ c_5y^2\ \text{at}\ (x+h,y)\ \text{minus at}\ (x,y)After
0,\ c_1h,\ 0,\ c_3(2xh+h^2),\ c_4hy,\ 0
Conditions: none
Rule: distributive · Uses: S02, power_laws, distributive
Leaf basis: distributive + power_laws + add_inverse
Subtract term by term: the constant, c_2y and c_5y² do not involve x and cancel to 0 because y is fixed; c_1(x+h)−c_1x=c_1h; c_3((x+h)²−x²)=c_3(2xh+h²), using (x+h)²=x²+2xh+h²; c_4((x+h)y−xy)=c_4hy.
S03.2 · equivalent
Before
c_1h+c_3(2xh+h^2)+c_4hy
After
h(c_1+2c_3x+c_4y)+c_3h^2
Conditions: none
Rule: distributive · Uses: S03.1, distributive
Leaf basis: distributive + mul_comm + add_comm
Factor h out of the terms of first degree in h: c_1h+2c_3xh+c_4yh=h(c_1+2c_3x+c_4y), leaving c_3h². This proves Lemma L_numerator_x.
S05.1 · implies
Before
\varepsilon>0
After
\delta=\frac{\varepsilon}{|B|+1}>0Conditions: \varepsilon>0
Rule: exists_intro · Uses: func_limit_def, abs_rules
Leaf basis: exists_intro + abs_rules + order_basics
Take any ε>0. |B|+1≥1>0, so δ=ε/(|B|+1) is positive. The +1 avoids division by 0 when B=0.
S05.2 · implies
Before
0<|h|<\delta
After
|(A+Bh)-A|=|B||h|\le|B|\delta=\frac{|B|}{|B|+1}\varepsilon<\varepsilonConditions: 0<|h|<\delta
Rule: abs_rules · Uses: S05.1, abs_rules, order_basics
Leaf basis: abs_rules + order_basics + forall_intro
|(A+Bh)−A|=|Bh|=|B||h|. |h|<δ and |B|≥0 give |B||h|≤|B|δ=|B|ε/(|B|+1). Since |B|<|B|+1, the ratio |B|/(|B|+1)<1, and multiplying by the positive ε gives <ε. ε was arbitrary, so the limit is A; S05.1 and S05.2 prove Lemma L_affine_limit.
Error controls and negative controls
Do not move both coordinates
One misconception named by the MKO is to change both coordinates at once. In S12, (f(1+h,2+h)−f(1,2))/h→5, while the x-partial is 4. The definition of the x-partial requires y to stay fixed, with h added to the first coordinate only.
One quotient is not the partial
Another misconception named by the MKO is to claim a partial exists from one nonzero-h quotient. Here q_x(h)=4+h is never 4 for nonzero h (h=1 gives 5, h=1/10 gives 41/10). Even when some finite difference happens to give the right number (for example the central difference (q_x(1)+q_x(−1))/2=(5+3)/2=4), that is a feature of quadratics, not a proof of existence: for u(x,y)=|x| at (0,0), the central difference (|h|−|−h|)/(2h)=0 for every h≠0, yet the x-partial does not exist (the right quotient is 1 and the left quotient is −1).
Existence of partials does not imply continuity
MKO assumption a2 warns that existing partials are not used to infer continuity of general functions. Example: g(x,y)=xy/(x²+y²) for (x,y)≠(0,0) and g(0,0)=0. g(x,0)=0 and g(0,y)=0 for all x, y, so both partials exist at the origin and equal 0; but along y=x, g(t,t)=1/2 for all t≠0, which does not tend to g(0,0)=0, so g is not continuous at the origin.
- q_x(1)=5, so f_x(1,2)=5 →
reject — The partial is the limit as h→0: q_x(h)=4+h→4 (S10, S11). - (f(1+h,2+h)−f(1,2))/h→5, so the x-partial is 5 →
reject — This moves y as well; with y=2 fixed the limit is 4 (S12). - If both partials of a function exist at a point, it is continuous there →
reject — g(x,y)=xy/(x²+y²), g(0,0)=0 has both partials 0 at the origin, but g(t,t)=1/2 (error control 3). - If the central difference (q(h)+q(−h))/2 gives a number for some h, the partial exists and equals it →
reject — For u(x,y)=|x| at (0,0) the central difference is always 0, yet the x-partial does not exist (error control 2).