%% Characterization case thickness
syms ro
pi = 3000;
po = 14.7;
ri = 1.75;
yield = 35000;
hoop_str = ((pi*ri^2-po*ro^2)-(ro^2*(po-pi)))/(ro^2-ri^2) == yield;
sol = double(solve(hoop_str, ro))