test1

%[xv yv] = meshgrid(linspace(-5,-3,21),linspace(-4,-1,21));
[xv yv] = meshgrid(linspace(-5,-3,21),linspace(-4,3,21));
for n=1:21
    for m=1:21
        [y q] = fval4([xv(n,m) yv(n,m)],imgpt);
        z1(n,m) = q(1);
        z2(n,m) = q(2);
        z3(n,m) = y;
    end
end
Undefined function or variable 'imgpt'.

Error in ==> test1 at 7
        [y q] = fval4([xv(n,m) yv(n,m)],imgpt);
mesh(xv,yv,z1)
xlabel('x-axis');
ylabel('y-axis');
grid
mesh(xv,yv,z2)
xlabel('x-axis');
ylabel('y-axis');
grid
mesh(xv,yv,log10(z3))
xlabel('x-axis');
ylabel('y-axis');
grid
v = [-3.7 -2];
[vf fmin exitflag options] = fminsearch(@fval4,v,[],imgpt)
vfit = [0 vf(1) 0 vf(2)];
[pto Tn pn] = back_projection(vfit,imgpt2);
draw5(pto);
%axis([-60 60 -120 0])