In this paper, we investigate joint optimization of secure AirComp and reliable multicasting assisted by a multipleinput multiple-output untrusted two-way relay, where artificial noise is employed at the access point (AP) to interfere the relay for ensuring secure AirComp. We aim to minimize the computation distortion at the AP by jointly designing the transmit variables of all the nodes and the aggregation variables at the AP and relay, under the secure AirComp constraint, the reliable multicasting constraint, and the transmit power constraints of all the nodes. We consider two scenarios that perfect and imperfect channel state information (CSI) are available. For the former, the formulated optimization problem is highly nonconvex, and we propose an effective block coordinate descent (BCD)penalty successive convex approximation (penaltySCA) method to solve the nonconvex problem. For the latter, we model the imperfect CSI by using worst-case criterion and the formulated robust optimization problem is much more challenging than its counterpart with perfect CSI. To solve the robust problem effectively, we first transform it to a deterministic optimization problem by employing some powerful mathematical lemmas, and then apply the proposed BCD-penaltySCA method to solve the reformulated deterministic problem. The proposed methods are shown by simulations to significantly reduce the computation distortion compared with other benchmarks under considering secure AirComp and reliable multicasting.