Hi I am currently tackling on how to find points of intersection on R studio, including a plot of the functions on the same plot. As you can see I have wrote out the constants and the simultaneous equation that needs to be solved to find theta, in order to find the value for 'k' and struggling on how to attempt this. Any help would be much appreciated!

#Constants given

L<-1.2*10^(-8)
u0<-1.6*10^(-19)

*232*10^(-3)

m_star<-9.109

*10^(-31)*

hbar<-(6.6310^(-34))/(2*pi)

hbar<-(6.63

me<-0.067*9.109*10^(-31)

mh<-0.5*9.109*10^(-31)

#Equation for theta (unknown value is k)

theta<-(k*L)/2
theta0<-sqrt(((m_star)*(u0)

*(L^(2)))/(2*(hbar^(2))))

#Simultaneous equations to solve

y<-theta

y<-theta0*cos(theta)