function J = yst_lstm(p, tdata, pdata)
% Least Squares fit to Logistic Growth
N = length(tdata);
yst = p(1)*exp(p(2)*tdata);
err = pdata - yst;
J = err*err'; % Sum of square errors
end
