We have investigated the femtosecond soliton propagation in inhomogeneous fiber, which is described by the modified inhomogeneous Hirota equation with variable coefficient (MIH-vc). With the aid of AKNS method, corresponding Lax pair is constructed. By virtue of the Darboux transformation method and symbolic computation, the analytic one- and two-soliton solutions are explicitly obtained. Using obtained solutions, we graphically discuss the features of femtosecond solitons in modified inhomogeneous Hirota system by changing the profile of variable coefficients. We analyze various form of group velocity dispersion, third order dispersion and nonlinearity parameter for periodic amplification system, exponentially distributed system, parabolic solitons, periodic exponentially modulated system, which will be observable in the future experiments. These results are potentially useful in future experiments and soliton control for long-distance optical communication. Finally, the soliton solutions of the MIH-vc equation in double Wronskian form is constructed and further verified using the Wronskian technique by substitute in bilinear equations.