function out1 = fun(x)
    out1 = x.^3-2*x-5;