Untitled

x0= [100; 100; 100; 100];
[x, fval]= fsolve(@sistema_equazioni, x0);
function F= sistema_equazioni(x)
A = [0.37 4.08 3.05 4.51;
    0.04 0.09 0.006 0.067;
    0.05 0.7 0.58 0.74;
    0.001 0.08 0.002 0.14];
b= [476; 10; 80; 10];

%x = A\b;
F= A*x-b; %perchè devo azzerarlo
end 

fprintf('grammi di verdure: %.2f g\n', x(1));
fprintf('grammi di pane: %.2f g\n', x(2));
fprintf('grammi di dessert: %.2f g\n', x(3));
fprintf('grammi di pesce: %.2f g\n', x(4));

ris=A*x; %per verificare che sia giusto