how to create this function \min _{\rho, \kappa, \theta, \sigma} \sum_{i=1}^n\left[b_{i t}-\sum_{k=1}^{\lambda_i} p_k\left(\rho_t, \kappa_t, \theta_t, \sigma_t, \tau_{i k}\right) c_{i k}\right]^2, in matlab

main.m
% Define constants
b = [...]; % vector of size n x 1
p = @(rho, kappa, theta, sigma, tau) [...]; % function handle, returns a vector of size lambda_i x 1
c = [...]; % matrix of size n x lambda_i
lambda = [...]; % vector of size n x 1
tau = [...]; % matrix of size n x lambda_i

% Define function to minimize as anonymous function
fun = @(x) sum((b - sum(p(x(1), x(2), x(3), x(4), tau(:,1:lambda(i))) .* c, 2)).^2);

% Use optimization function (e.g. fminsearch) to minimize the function
x0 = [...]; % initial guess for values of rho, kappa, theta, sigma
x = fminsearch(fun, x0);
579 chars
14 lines