bar_tension.sci
clc clear P = 1000.0 // N L = 2.0 //m E = 200e9 //GPa b = 0.1 //m h = 0.2 //m A = b*h q = 10.0 // N/m x = 0:0.1:L u = 1.0/(E*A)*(-0.5*q*x.^2+(P+q*L)*x) plot(x,u,'r')