function z = my_th(x)
I = imread('cameraman.tif');
L = 256;
[m,n]=size(I);
N = m*n;
p=imhist(I);
p=p/N;
ip = p;
for i = 0:L-1
ip(i+1) = i * ip(i+1);
end
uT = sum(ip);
w0 = sum(p(1:x+1));
wt = w0;
w1 = 1 - wt;
ut = sum(ip(1:x+1));
u0 = ut/(wt+eps);
u1 = (uT-ut)/(1-wt+eps);
z = -(w0*w1*(u1-u0)^2);