% Material properties
materials = {'Mild Steel', 'AHSS', 'Aluminum', 'Magnesium', 'CFRP', ...
             'Glass Fiber Comp.', 'Natural Fiber Comp.', 'Polycarbonate', 'Titanium'};
density = [7.85, 7.85, 2.70, 1.80, 1.60, 1.90, 1.30, 1.20, 4.50];
yield_strength = [250, 600, 280, 200, 500, 250, 120, 70, 900];

% Different marker shapes for clarity
markers = {'o', 's', 'd', '^', 'v', '*', 'x', 'p', 'h'};
figure; hold on;

for i = 1:length(materials)
    plot(density(i), yield_strength(i), markers{i}, 'MarkerSize', 10, 'LineWidth', 2, ...
        'DisplayName', materials{i});
    % Add label next to each point
    text(density(i)+0.07, yield_strength(i), materials{i}, 'FontSize', 9, ...
        'VerticalAlignment', 'bottom', 'HorizontalAlignment', 'left');
end

xlabel('Density (g/cm^3)');
ylabel('Yield Strength (MPa)');
title('Ashby Diagram: Density vs. Yield Strength');
grid on;
legend('Location', 'northeastoutside');
hold off;