Help with an implicit solution of a system.

Hi, guys.

I have the system below:

V(:,:,n)=(D+0.5*deltat(1,n)*(-K))\(D*m1(:,:,n)-(-K)*m(:,:,n)+Fload);

which I have to solve implicitly. "K" and "D" are constant "81x81" matrices. "Fload" and "m" are constant "81x1" matrices and my time step is an adaptive scalar which changes each round. So the solution should be a "81x1". The problem is for some time steps my solution gets so close to singularity that MATLAB warns me. What can I do about it? is there any other solvers besides "\" that have more accuracy?