Yes!
This works:
% MAPPING: d(1) = d_out, d(2) = d_in
pipe = @(d,Iz,A) [pi*(d(1)^4-d(2)^4)/64 - Iz; pi*(d(1)^2-d(2)^2)/4 - A];
Iz = 100; % Create Data
A = 25; % Create Data
D = fsolve(@(d) pipe(d,Iz,A), [1; 1]);
You need the Optimization Toolbox to use the fsolve function.
