chisq_skydip

mod1 = (exp(-x(1)./sind(el)) - exp(-x(1)./sind(el0)))./(exp(-x(1)./sind(el)) - ...
    exp(-x(1)./sind(el1)));

0001 function delta = chisq_skydip(x, yvals, el, el0, el1)
0002 
0003 %mod1 = (exp(-x(1)./sind(el)) - exp(-x(1)./sind(el0)))./(exp(-x(1)./sind(el)) - ...
0004 %    exp(-x(1)./sind(el1)));
0005 
0006 %mod2 = (exp(-x(2)./sind(el)) - exp(-x(2)./sind(el0)))./(exp(-x(2)./sind(el)) - ...
0007 %    exp(-x(2)./sind(el1)));
0008 
0009 %mod1 = (x(1,1) + x(1,2)./sind(el))./(x(1,1) + x(1,2)./sind(el0));
0010 %mod2 = (x(2,1) + x(2,2)./sind(el))./(x(2,1) + x(2,2)./sind(el0));
0011 
0012 %mod1 = exp(-x(1)./sind(el))./exp(-x(1)./sind(el0)));
0013 %mod2 = (exp(-x(2)./sind(el))./exp(-x(2)./sind(el0)));
0014 
0015 mod1 = (exp(-x(1)./sind(el0)) - ...
0016     exp(-x(1)./sind(el)))./(exp(-x(1)./sind(el1)) - exp(-x(1)./sind(el0)));
0017 
0018 mod2 = (exp(-x(2)./sind(el0)) - ...
0019     exp(-x(2)./sind(el)))./(exp(-x(2)./sind(el1)) - exp(-x(2)./sind(el0)));
0020 
0021 modval(:,1) = mod1;
0022 modval(:,2) = mod2;
0023 
0024 modval = modval(:);
0025 
0026 delta = (modval - yvals);
0027 
0028 
0029 return;
0030     
0031