To make better decisions on approximation, we may need to increase reliable and useful information on different aspects of approximation. To enhance information about the quality and certainty of approximating the solution of an Apollonius-type quadratic functional equation, we need to measure both the quality and the certainty of the approximation and the maximum errors. To measure the quality of it, we use fuzzy sets, and to achieve its certainty, we use the probability distribution function. To formulate the above problem, we apply the concept of Z-numbers and introduce a special matrix of the form $\mathrm{diag}(A, B, C)$
diag
(
A
,
B
,
C
)
(named the generalized Z-number) where A is a fuzzy time-stamped set, B is the probability distribution function, and C is a degree of reliability of A that is described as a value of $A\ast B$
A
∗
B
. Using generalized Z-numbers, we define a novel control function to investigate H–U–R stability to approximate the solution of an Apollonius-type quadratic functional equation with quality and certainty of the approximation.