matlab code -- 2

I am looking for someone to help me develop some MATLAB code to solve a specific problem. The code should be of intermediate complexity and I have particular requirements and constraints that need to be met. The code will help with data analysis, image processing, and optimization. I am looking for someone with expertise and experience in developing MATLAB code, so please only submit proposals if you meet this criterion and the due date for this assignment is tomorrow 11 am. I already have the code partially for question 1 however I am not sure if does the job : % Crater diameter samples between 0.4 and 2.7 cm D = 0.4:0.001:2.7; % Material density values p = [2.7 3.2 4.51 8.0]; % Al SiC Ti Steel density pAl = 2.7; pSiC = 3.2; pTi = 4.51; pSteel = 8.0; % g/cm^3 % Angle of collision in radians and ISS velocity in km/s teta = 0.698; V = 8; % Calculate impactor diameters (dp) for each material dpAl = (D ./ (30.9 * pAl^(0.44) * (V * cosd(teta))^(0.44)).^(1/1.33)); dpSiC = (D ./ (30.9 * pSiC^(0.44) * (V * cosd(teta))^(0.44)).^(1/1.33)); dpTi = (D ./ (30.9 * pTi^(0.44) * (V * cosd(teta))^(0.44)).^(1/1.33)); dpSteel = (D ./ (30.9 * pSteel^(0.44) * (V * cosd(teta))^(0.44)).^(1/1.33)); %% Impactor vs Crater diameter figure; plot(dpAl, D, 'LineWidth', 1.5); hold on; plot(dpSiC, D, 'LineWidth', 1.5); plot(dpTi, D, 'LineWidth', 1.5); plot(dpSteel, D, 'LineWidth', 1.5); hold off; xlabel("Impactor Diameter (cm)"); ylabel("Crater Diameter (cm)"); title("Impactor Diameter vs Crater Diameter"); legend("Al", "SiC", "Ti", "Steel"); %% Crater vs velocity figure; velocity = 3:0.0001:11; % Velocity samples between 3 and 11 km/s DAlv = (30.9 * pAl^(0.44) * dpAl.^(1.33) .* (velocity * cosd(teta)).^(0.44)); DSiCv = (30.9 * pSiC^(0.44) * dpSiC.^(1.33) .* (velocity * cosd(teta)).^(0.44)); DTiv = (30.9 * pTi^(0.44) * dpTi.^(1.33) .* (velocity * cosd(teta)).^(0.44)); DSteelv = (30.9 * pSteel^(0.44) * dpSteel.^(1.33) .* (velocity * cosd(teta)).^(0.44)); plot(velocity, DAlv, 'LineWidth', 1.5); hold on; plot(velocity, DSiCv, 'LineWidth', 1.5); plot(velocity, DTiv, 'LineWidth', 1.5); plot(velocity, DSteelv, 'LineWidth', 1.5); hold off; xlabel("Impactor Velocity (km/s)"); ylabel("Crater Diameter (mm)"); title("Impactor Velocity vs Crater Diameter"); legend("Al", "SiC", "Ti", "Steel"); %% Impactor angle figure; tetavariable = 0:0.0001:pi/2; % Variables between 0 and pi/2 DAlteta = (30.9 * pAl^(0.44) * dpAl.^(1.33) .* (V * cosd(tetavariable)).^(0.44)); DSiCteta = (30.9 * pSiC^0.44 * dpSiC.^(1.33) .* (V * cosd(tetavariable)).^(0.44)); DTiteta = (30.9 * pTi^0.44 * dpTi.^(1.33) .* (V * cosd(tetavariable)).^(0.44)); DSteelteta = (30.9 * pSteel^(0.44) * dpSteel.^(1.33) .* (V * cosd(tetavariable)).^(0.44)); plot(tetavariable, DAlteta, 'LineWidth', 1.5); hold on; plot(tetavariable, DSiCteta, 'LineWidth', 1.5); plot(tetavariable, DTiteta, 'LineWidth', 1.5); plot(tetavariable, DSteelteta, 'LineWidth', 1.5); hold off; xlabel("Impact Angle (radians)"); ylabel("Crater Diameter (mm)"); title("Impact Angle vs Crater Diameter"); legend("Al", "SiC", "Ti", "Steel");