Applied Mechanics-Buckling

import numpy as np

import matplotlib.pyplot as plt

n=1; L=5;

x = np.linspace(0, L, 100)

y = 1.0*np.sin(n*np.pi*x/L)

plt.plot(y,x)

n=2

y = 1.0*np.sin(n*np.pi*x/L)

plt.plot(y,x)

n=3

y = 1.0*np.sin(n*np.pi*x/L)

plt.plot(y,x)

E=10000; n=1; A=10

lamda=np.linspace(0.1, 0.4, 100)

Pcr = n*np.pi**2*E*A/(lamda**2)

plt.plot(lamda,Pcr)