How to save outputs of a function to different matrixs

function F=fxyz(X)
  x=X(1);
  y=X(2);
z=X(3);
F(1)=x+y+z;
F(2)=3*x+5*y+6*z;
F(3)=x-3*y-6*z-1;
end

[B(1,2),C(1,3),D(1,1)]=fsolve('fxyz',[-1,1,-1])
Assignment has more non-singleton rhs dimensions than non-singleton subscripts

Answer 1

[Stalin Samuel]

The function 'xyz' returns the vector 'F' which has three elements.So follow the below steps

F =fsolve('fxyz',[-1,1,-1]);
B(1,2)=F(1);
C(1,3)=F(2);
D(1,1)=F(3);