Q4 CODE
function f=question4(x)
N1=x(1);
W2=x(2);
H2=x(3);
r1=x(4);
r2=x(5);
r3=x(6);
Q=x(7);
f(1) = 0.25*N1-r1-r2-r3-0.2*1250;
f(2) = r3-1250+362.5+H2+W2;
f(3) = 0.75*N1-2*r1-3*r2+r3-H2;
f(4) = r2-r3-W2;
f(5) = r1-100;
f(6) = r2-12.5;
f(7) = r1*-24773.8 + r2*-52572.8 + r3*-12753.7 + 303.75*N1 - Q;
end
script file:
x0=10*ones(1,7);
x=fsolve(@question4,x0);
fprintf('Q = %8.2f\n',x(7));
N1=x(1);
r1=x(4);
r2=x(5);
r3=x(6);
Tplot=550:50:700;
Q = r1*-24773.8 + r2*-52572.8 + r3*-12753.7 - ((2.025*N1)*(Tplot-700));
plot (Q,Tplot)
xlabel('Heat transfer rate');
ylabel('Temperature');
title('Q vs T');