//Example 23 // force constant and work done
clc;
clear;
close;
m1=1;//assume
m2=19;//assume
mh=1.66*10^-27;//kg
mu=((m1*m2)/(m1+m2))*mh;//kg
w=7.55*10^14;//radians/sec
k=mu*(w)^2;//N/m
disp(k,"force constant is,(N/m)=")
x=0.5;//arngstrom
wh=((1/2)*k*(x*10^-10)^2);//joule
disp(wh,"work done is ,(joule)=")