/*  10-4-2023   FP bias   */

fx = f0 + x f1

p0 = Integrate[ fx, {x,-h,0} ]
p1 = Integrate[ fx, {x, 0,h} ]

m0 = (1/2 - x/h)
m1 = (1/2 + x/h)

ev$x = m0 n p0 / ( n h ) + m1 n p1 / ( n h )

Expand[ ev$x ]
bias = ev$x - fx

fx = f0 + x f1 + x^2 f2 /2

ibs$0 = Expand[ Integrate[ bias^2, {x, -h/2, h/2 } ] ]




var$x = m0^2/(n h)^2 n p0 (1-p0) +
        m1^2/(n h)^2 n p1 (1-p1) +
	2 m0 m1 / (n h)^2 (-n p0 p1)

var$0 = Integrate[ var$x, {x, -h/2, h/2}]