l=1000; %mm
a=50; %mm^2
e=100; %KN/mm^2
h=0;X=(0:10:500);Y=[];o=h/l;lp=l/cos(o);ld=(lp+l)/l;while h < 500 t=e.*a.*ld./l; p=2.*t.*h./l; o=h./l; h=h+10; lp=l./cos(o); ld=(lp+l)./l; Y(h)=p; endplot(x,y,'r-')
Here is my code. Am I missing something? Could someone fix it and explain to me, thank you.
Best Answer