function [xpred, Ppred] = predict_nonlin(process_func, x, u, P, Q)

[xpred, df] = feval(process_func, x, u);
Ppred = df*P*df' + Q;


