Problem with FORM analysis

Hello everyone,

I tried to analyze the reliability of a limit state function using the Monte Carlo and FORM methods. The solution of the Monte Carlo method is Pf=1.2850e-03. But the FORM method has a zero probability of failure, which clearly indicates a calculation error.
I know this is not a UQLab software issue. But I wanted to ask experts in this field to help me interpret this strange computational difference.

Why does the FORM method not work for the function I have defined ?!

I have enclosed the codes, here:

uq_Example_Reliability_07_RubOver.m (2.3 KB)

uq_OvertoppingFailure2.m (364 Bytes)

Thanks in advance.