clc
clear
calib=input('enter calib matrix')
meas=input('enter meas matrix')
n=size(calib,2)/2
l=calib'
syms a b cx cy
for k=1:2:2*n
x=a*meas(1,k)+b*meas(1,k+1)+cx
y=-b*meas(1,k)+a*meas(1,k+1)+cy
A(k,:)=jacobian(x,[a b cx cy])
A(k+1,:)=jacobian(y,[a b cx cy])
end
A=eval(A)
xcap=inv(A'*A)*(A'*l)
r=l-(A*xcap)