%% Stress field

for i = 1:nelm
    if sub_D(i) == 1 || sub_D(i) == 4 || sub_D(i) == 5   %Aluminium
        D = hooke(2, E_al, v_al);
    elseif sub_D(i) == 2       %Copper
        D = hooke(2, E_cu, v_cu);
    elseif sub_D(i) == 3       %Electricity Core
        D = hooke(2, E_core, v_core);
    elseif sub_D(i) == 6       %Steel
        D = hooke(2, E_steel, v_steel);
    end
    [es, et] = plants(Ex(i,:), Ey(i,:), ep,D, ed(i,:));
    Es(i,:) = es; %Element stresses
    Es(i,1:3) = Es(i,1:3)-Sigma_0(i,1:3); %Subtracting thermal strain
end