let f be PartFunc of REAL ,REAL ; for x0 being Real st f is_right_differentiable_in x0 & ex r0 being Real st
( r0 > 0 & ( for g being Real st g in dom f & g in [.x0,(x0 + r0).] holds
f . g <> 0 ) ) holds
( f ^ is_right_differentiable_in x0 & Rdiff (f ^ ),x0 = - ((Rdiff f,x0) / ((f . x0) ^2 )) )
let x0 be Real; ( f is_right_differentiable_in x0 & ex r0 being Real st
( r0 > 0 & ( for g being Real st g in dom f & g in [.x0,(x0 + r0).] holds
f . g <> 0 ) ) implies ( f ^ is_right_differentiable_in x0 & Rdiff (f ^ ),x0 = - ((Rdiff f,x0) / ((f . x0) ^2 )) ) )
A1:
0 + x0 = x0
;
assume A2:
f is_right_differentiable_in x0
; ( for r0 being Real holds
( not r0 > 0 or ex g being Real st
( g in dom f & g in [.x0,(x0 + r0).] & not f . g <> 0 ) ) or ( f ^ is_right_differentiable_in x0 & Rdiff (f ^ ),x0 = - ((Rdiff f,x0) / ((f . x0) ^2 )) ) )
then consider r2 being Real such that
A3:
0 < r2
and
A4:
[.x0,(x0 + r2).] c= dom f
by Def3;
given r0 being Real such that A5:
r0 > 0
and
A6:
for g being Real st g in dom f & g in [.x0,(x0 + r0).] holds
f . g <> 0
; ( f ^ is_right_differentiable_in x0 & Rdiff (f ^ ),x0 = - ((Rdiff f,x0) / ((f . x0) ^2 )) )
set r3 = min r0,r2;
0 <= min r0,r2
by A5, A3, XXREAL_0:15;
then A7:
x0 <= x0 + (min r0,r2)
by A1, XREAL_1:8;
min r0,r2 <= r2
by XXREAL_0:17;
then A8:
x0 + (min r0,r2) <= x0 + r2
by XREAL_1:9;
then
x0 <= x0 + r2
by A7, XXREAL_0:2;
then A9:
x0 in [.x0,(x0 + r2).]
by XXREAL_1:1;
x0 + (min r0,r2) in { g where g is Real : ( x0 <= g & g <= x0 + r2 ) }
by A7, A8;
then
x0 + (min r0,r2) in [.x0,(x0 + r2).]
by RCOMP_1:def 1;
then
[.x0,(x0 + (min r0,r2)).] c= [.x0,(x0 + r2).]
by A9, XXREAL_2:def 12;
then A10:
[.x0,(x0 + (min r0,r2)).] c= dom f
by A4, XBOOLE_1:1;
min r0,r2 <= r0
by XXREAL_0:17;
then A11:
x0 + (min r0,r2) <= x0 + r0
by XREAL_1:9;
then
x0 <= x0 + r0
by A7, XXREAL_0:2;
then A12:
x0 in [.x0,(x0 + r0).]
by XXREAL_1:1;
x0 + (min r0,r2) in { g where g is Real : ( x0 <= g & g <= x0 + r0 ) }
by A7, A11;
then
x0 + (min r0,r2) in [.x0,(x0 + r0).]
by RCOMP_1:def 1;
then A13:
[.x0,(x0 + (min r0,r2)).] c= [.x0,(x0 + r0).]
by A12, XXREAL_2:def 12;
A14:
[.x0,(x0 + (min r0,r2)).] c= dom (f ^ )
A20:
x0 in [.x0,(x0 + (min r0,r2)).]
by A7, XXREAL_1:1;
A21:
for h being convergent_to_0 Real_Sequence
for c being V8() Real_Sequence st rng c = {x0} & rng (h + c) c= dom (f ^ ) & ( for n being Element of NAT holds h . n > 0 ) holds
( (h " ) (#) (((f ^ ) /* (h + c)) - ((f ^ ) /* c)) is convergent & lim ((h " ) (#) (((f ^ ) /* (h + c)) - ((f ^ ) /* c))) = - ((Rdiff f,x0) / ((f . x0) ^2 )) )
proof
let h be
convergent_to_0 Real_Sequence;
for c being V8() Real_Sequence st rng c = {x0} & rng (h + c) c= dom (f ^ ) & ( for n being Element of NAT holds h . n > 0 ) holds
( (h " ) (#) (((f ^ ) /* (h + c)) - ((f ^ ) /* c)) is convergent & lim ((h " ) (#) (((f ^ ) /* (h + c)) - ((f ^ ) /* c))) = - ((Rdiff f,x0) / ((f . x0) ^2 )) )let c be
V8()
Real_Sequence;
( rng c = {x0} & rng (h + c) c= dom (f ^ ) & ( for n being Element of NAT holds h . n > 0 ) implies ( (h " ) (#) (((f ^ ) /* (h + c)) - ((f ^ ) /* c)) is convergent & lim ((h " ) (#) (((f ^ ) /* (h + c)) - ((f ^ ) /* c))) = - ((Rdiff f,x0) / ((f . x0) ^2 )) ) )
assume that A22:
rng c = {x0}
and A23:
rng (h + c) c= dom (f ^ )
and A24:
for
n being
Element of
NAT holds
h . n > 0
;
( (h " ) (#) (((f ^ ) /* (h + c)) - ((f ^ ) /* c)) is convergent & lim ((h " ) (#) (((f ^ ) /* (h + c)) - ((f ^ ) /* c))) = - ((Rdiff f,x0) / ((f . x0) ^2 )) )
A25:
lim h = 0
by FDIFF_1:def 1;
A26:
for
m being
Element of
NAT holds
c . m = x0
A27:
(dom f) \ (f " {0 }) c= dom f
by XBOOLE_1:36;
rng (h + c) c= (dom f) \ (f " {0 })
by A23, RFUNCT_1:def 8;
then A28:
rng (h + c) c= dom f
by A27, XBOOLE_1:1;
then A29:
lim ((h " ) (#) ((f /* (h + c)) - (f /* c))) = Rdiff f,
x0
by A2, A22, A24, Th15;
Rdiff f,
x0 = Rdiff f,
x0
;
then A30:
(h " ) (#) ((f /* (h + c)) - (f /* c)) is
convergent
by A2, A22, A24, A28, Th15;
then A31:
- ((h " ) (#) ((f /* (h + c)) - (f /* c))) is
convergent
by SEQ_2:23;
x0 in dom (f ^ )
by A20, A14;
then A32:
x0 in (dom f) \ (f " {0 })
by RFUNCT_1:def 8;
rng c c= (dom f) \ (f " {0 })
then A33:
rng c c= dom (f ^ )
by RFUNCT_1:def 8;
then A34:
f /* c is
non-zero
by RFUNCT_2:26;
then A36:
h (#) ((h " ) (#) ((f /* (h + c)) - (f /* c))) = (f /* (h + c)) - (f /* c)
by FUNCT_2:113;
A37:
f /* (h + c) is
non-zero
by A23, RFUNCT_2:26;
then A38:
(f /* (h + c)) (#) (f /* c) is
non-zero
by A34, SEQ_1:43;
then A39:
(f /* c) + ((f /* (h + c)) - (f /* c)) = f /* (h + c)
by FUNCT_2:113;
dom (f ^ ) = (dom f) \ (f " {0 })
by RFUNCT_1:def 8;
then A40:
dom (f ^ ) c= dom f
by XBOOLE_1:36;
A41:
for
g being
real number st
0 < g holds
ex
n being
Element of
NAT st
for
m being
Element of
NAT st
n <= m holds
abs (((f /* c) . m) - (f . x0)) < g
then A43:
f /* c is
convergent
by SEQ_2:def 6;
then A44:
lim (f /* c) = f . x0
by A41, SEQ_2:def 7;
A45:
h is
convergent
by FDIFF_1:def 1;
then
h (#) ((h " ) (#) ((f /* (h + c)) - (f /* c))) is
convergent
by A30, SEQ_2:28;
then A46:
f /* (h + c) is
convergent
by A43, A36, A39, SEQ_2:19;
lim (h (#) ((h " ) (#) ((f /* (h + c)) - (f /* c)))) =
(lim h) * (lim ((h " ) (#) ((f /* (h + c)) - (f /* c))))
by A45, A30, SEQ_2:29
.=
0
by A25
;
then
0 = (lim (f /* (h + c))) - (f . x0)
by A43, A44, A36, A46, SEQ_2:26;
then A47:
lim ((f /* (h + c)) (#) (f /* c)) = (f . x0) ^2
by A43, A44, A46, SEQ_2:29;
A48:
lim ((f /* (h + c)) (#) (f /* c)) <> 0
then A50:
(h " ) (#) (((f ^ ) /* (h + c)) - ((f ^ ) /* c)) = (- ((h " ) (#) ((f /* (h + c)) - (f /* c)))) /" ((f /* (h + c)) (#) (f /* c))
by FUNCT_2:113;
A51:
(f /* (h + c)) (#) (f /* c) is
convergent
by A43, A46, SEQ_2:28;
then lim ((h " ) (#) (((f ^ ) /* (h + c)) - ((f ^ ) /* c))) =
(lim (- ((h " ) (#) ((f /* (h + c)) - (f /* c))))) / ((f . x0) ^2 )
by A38, A47, A48, A31, A50, SEQ_2:38
.=
(- (Rdiff f,x0)) / ((f . x0) ^2 )
by A30, A29, SEQ_2:24
.=
- ((Rdiff f,x0) / ((f . x0) ^2 ))
by XCMPLX_1:188
;
hence
(
(h " ) (#) (((f ^ ) /* (h + c)) - ((f ^ ) /* c)) is
convergent &
lim ((h " ) (#) (((f ^ ) /* (h + c)) - ((f ^ ) /* c))) = - ((Rdiff f,x0) / ((f . x0) ^2 )) )
by A38, A51, A48, A31, A50, SEQ_2:37;
verum
end;
0 < min r0,r2
by A5, A3, XXREAL_0:15;
hence
( f ^ is_right_differentiable_in x0 & Rdiff (f ^ ),x0 = - ((Rdiff f,x0) / ((f . x0) ^2 )) )
by A14, A21, Th15; verum