“a state variable approach with numerical implementation”