changed small bug in calc_F.m (length(rL))