[Scilab-users] Using lsqrsolve with a certain kind of function

[Viktor Goryainov]

Hello,


I have a function that takes tabulated values from several arrays 
(a_w15, b_w15, etc.) and thus returns a vector:


function r = vr3(C, mu)
  a = a_w15 + a_chl15 * C(1) + a_sm15 * C(2) + a_doc15 * C(3);
  bb = b_w15 + b_chl15 * C(1) + b_sm15 * C(2);
  dd = bb ./ a;


  if dd < 0.25 then
   r = 0.319 * dd / mu;
  else
   r = (0.267 * dd + 0.013) / mu;
  end
endfunction

Now I need to use it twice in my task: first to generate a vector for 
some C_inh, then to try and restore C_inh from an initial guess C_init 
using lsqrsolve. How should I modify the vr3 to achieve this?

Best regards,
Viktor