IN THE FOLLOWING IS IT POSSIBLE TO MODIFY THE CODE SOTHAT IT WILL GIVE ME OUTPUT FOR DIFFERENT VALUE OF 'b' (LIKE ATTACHING A FOR LOOP IN THAT AS WE DO IN C++.)Clear[k, \[Kappa], a, b, V, \[Sigma], T, g]
\[Sigma] = 1; a = 0.1; b = 26; V = 3;
\[Kappa] = Sqrt[V] - k*k;
T =...
I am giving the following command for calculation of transmission probability through a double barrier but mathematica is not giving any solution to this infact it is taking a lot of time saying 'running'.
I want to know that whether there is something wrong with the command lines(I am sure...
Could anyone tell me where can i find more about the fact that if your Hamiltonian is non-central then the total energy is not dependent only on principle quantum number n but also on l.
i have one problem in numerical approach that how to calculate j(inc) as how shall i know that which part is incident and which part is reflected.
I also tried the transfer matrix approach but it is giving the correct result for one rectangular potential barrier but not for two.
i actually want to calculate the reflection coefficient for any potential.so could you suggest any numerical method to do this as i want to change the parameters (like height ,length of the potential) and get the corresponding T and R.
I have a question
instead of using
y1[a] == 1, y1'[a] == 1
as other two conditions if i will use
y1[-5]==1,y1'[5]==0
as my condition then will that work?
i am taking this condition as my wave func will be a plane wave for a<=0 and i can take a condition like above.
please tell me where did...
Thanks again for the replyI am now learning to use mathematica for my computation and simulation(i used to do these with C++).I want to write a code in mathematica using Numerov method to solve this problem.But i need to know how to write code in mathematica. so is there any quick references for...
Thanks for the reply
1--
could you please tell me why r u putting the condition
y1[a] == 1, y1'[a] == 1
when solving numerically.
2--
again even if u do it numerically the plots of y1 ,y2, y3 look same.
I am trying to solve the potential barrier problem numerically in Mathematica
I am giving the following command but it is showing error which i could not figure out
a = 2;
s = NDSolve[{y1''[x] + k1*y1[x] == 0, y2''[x] + k2*y2[x] == 0,
y3''[x] + k1*y3[x] == 0}, {y1[0] == y2[0], y2[a] ==...