function KriegF

     globals;
     
     KriegD;
     
     F = u - v + (mue/h^2)*(D*u);

end