In dr1.m you should have (lines 209 to 217):
if ~isreal(jacobia_)
if max(max(abs(imag(jacobia_)))) < 1e-15
jacobia_ = real(jacobia_);
else
info(1) = 6;
info(2) = sum(sum(imag(jacobia_).^2));
return
end
end
If the jacobian matrix is complex dr1 returns info(1)=6. Can you post (or send) your mod file if you do not fall in the trap ?
Stéphane.