Untitled

T0=273.15;
T4=373.15;
k=[3, 1.5, 5, 2];
z=[0, 0.1, 0.2, 0.4, 0.45];


fun= @(x) [k(1)*(x(1)-T0)/(z(2)-z(1)) -x(4); %eq1
    k(2)*(x(2)-x(1))/(z(3)-z(2)) -x(4); %eq2
    k(3)*(x(3)-x(2))/(z(4)-z(3)) -x(4); %eq3
    k(4)*(T4-x(3))/(z(5)-z(4)) -x(4)]; %eq4


x0=[280 300 320 10000];

sol = fsolve(fun, x0);

T1 = sol(1);
T2 = sol(2);
T3 = sol(3);
q = sol(4);

fprintf('T1= %.4f K\n', T1);
fprintf('T2= %.4f K\n', T2);
fprintf('T3= %.4f K\n', T3);
fprintf('q= %.4f W/m^2\n', q);