plotCDM_Part_2.m

%Problem 2a-d
clear;
clc;

%iterates through each number of segments and applied load to fin three
%vectors of values of deflection for each load
N = [6,10,20];
q = [1,10,20,30,50];
data = zeros(3,5);
for n = [1:3]
    for m =[1:5]
        data(n,m) = max(CDM(q(m),N(n),0));
    end
end

%Plotting routine for Part 2
plot(q,data(1,:),q,data(2,:),q,data(3,:));
ylabel('Deflection (meters)')
xlabel('Applied Load q (N/m)')
title('Deflection vs. Applied Load q')