Numerical solution of parquet equations for the Anderson impurity model