canitlever_beam_deflection.sci
clc clear P = 1000.0 // N L = 2.0 //m E = 200e9 // Pa b = 0.1 //m h = 0.2 //m I = b*h^3/12.0 x = 0 : 0.1 : L w1 = P/(6*E* I ) * x.^3 - P*L / (2.0* E* I ) * x.^2 P=2*P w2 = P/(6*E* I ) * x.^3 - P*L / (2.0* E* I ) * x.^2 plot(x ,w1, "r") plot(x ,w2, "b")