For reducing mathematical complexity in the solution of real-life systems, additional assumptions like zero initial condition, integer order system, etc. result in inaccurate modeling of systems, leading to erroneous estimation. This paper presents a complete recursive method for solving irrational systems considering nonzero initial values. Generalized operational integration matrices based on orthogonal HF have been utilized to solve the differential equations describing a fractional order system. The orthogonal HF set is effective in producing piecewise linear solutions of various sorts of control systems when it is used with function samples. A few illustrative instances are provided to confirm the effectiveness of the technique.