Fiber Recruitment Model Solve & Plot

Every 0.1mm, a new fiber with a stiffness 20 N/mm is recruited starting from x=0. The net force is given by the following equation: . If p=10, solve for F and plot F vs. x.

p=10; %number of fibers
k(1,1:p)=20; %stiffness (N/mm)
dx=0.1; %step size
x=zeros(p*dx); %allocate space for x-values
x0=0;
for i = 1:p*dx
    F=k(i)*(x(i)-(x0+i*dx));
    x=x(i)+i*dx;
end
plot(x,F)
xlabel('x')
ylabel('F')