Re: shooting method

• To: mathgroup at smc.vnet.net
• Subject: [mg92105] Re: shooting method
• From: "Kevin J. McCann" <Kevin.McCann at umbc.edu>
• Date: Sat, 20 Sep 2008 04:58:38 -0400 (EDT)
• Organization: University System of Maryland
• References: <gapfd4\$o7q\$1@smc.vnet.net>

```I think you mean that A is approximately zero not <<1.

To figure out what I should actually use for A and B, I first Plot your
potential. If I am looking for bound states, then it looks as though I
can use A=0.9 and B=4. Further refinements could be necessary. Here is
some code.

A = 0.9;
B = 4.0;
nde[e_] = {-y''[x] + 400 (1/x^12 - 1/x^6) y[x] - e y[x] == 0,
y[A] == 0, y'[A] == 1};
tbl = Table[{e,
y[B] /. NDSolve[nde[e],
y, {x, A, B}][]}, {e, -10, -0.01, .01}];
ListPlot[tbl, PlotRange -> {-100, 100}]

What this does is to start at A and shoot to B for a given trial
eigenenergy e, then evaluate y[B] and plot the results for a lot of e's.
As you can see if you run the code, there is a zero crossing for e ~=
-0.85 (the only zero crossing).

Now you can write a function and use FindRoot to get the eigenvalue,
then you can generate the eigenfunction and normalize it.

Here is the code to find the eigenvalue and plot it:

Clear[f];
f[e_?NumberQ] := y[B] /. NDSolve[nde[e], y, {x, A, B}][]
e0 = e /. FindRoot[f[e] == 0, {e, -9, -8}]
Y = y /. NDSolve[nde[e0], y, {x, A, B}][]

Plot[Y[x], {x, A, B}]

Note that there is another eigenvalue/eigenfunction at about e=-55 as well.

Luca Petrone wrote:
> Dear All,
>
> Using a shooting method, I would like to find the smallest eigenvalue of
>
> -y''[x] + 400((1/x)^12 - (1/x)^6) y[x] ==  E0*y[x]
>
> (the quantum-mechanical steady state energy of an anharmonic
oscillator with Lennard-Jones potential)
> with the boundary conditions :
> y[A]==0
> y[B]==0
> with A<<1 and B >> 1
> Can anyone help about how to implement it ?
>
> Thank you
>
> Luca
>

```

• Prev by Date: Re: NDSolve and Piecewise
• Next by Date: Re: Debracketing array symbols
• Previous by thread: Re: shooting method
• Next by thread: imagesize to full width of current screen
{"ft_lang_label":"__label__en","ft_lang_prob":0.8215586,"math_prob":0.9546239,"size":1879,"snap":"2019-35-2019-39","text_gpt3_token_len":642,"char_repetition_ratio":0.08693333,"word_repetition_ratio":0.011594203,"special_character_ratio":0.38265035,"punctuation_ratio":0.18262807,"nsfw_num_words":0,"has_unicode_error":false,"math_prob_llama3":0.9976231,"pos_list":[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14],"im_url_duplicate_count":[null,null,null,null,null,null,null,null,null,null,null,null,null,null,null],"WARC_HEADER":"{\"WARC-Type\":\"response\",\"WARC-Date\":\"2019-09-17T11:46:44Z\",\"WARC-Record-ID\":\"<urn:uuid:f35be768-06f0-4b6c-8c4d-285480c7804f>\",\"Content-Length\":\"43663\",\"Content-Type\":\"application/http; msgtype=response\",\"WARC-Warcinfo-ID\":\"<urn:uuid:8bdc400c-9c08-4d95-8851-7f33bf342a20>\",\"WARC-Concurrent-To\":\"<urn:uuid:b5c77ca9-8f99-4840-9fbe-dcaf1ae43c35>\",\"WARC-IP-Address\":\"140.177.205.73\",\"WARC-Target-URI\":\"http://forums.wolfram.com/mathgroup/archive/2008/Sep/msg00459.html\",\"WARC-Payload-Digest\":\"sha1:PE4MYCUKSN7H2FSZ2ZQ65ZBMGIKNUYPQ\",\"WARC-Block-Digest\":\"sha1:OE57FPZZSRSJWZW4KWFUBZLHI6XIOABT\",\"WARC-Identified-Payload-Type\":\"text/html\",\"warc_filename\":\"/cc_download/warc_2019/CC-MAIN-2019-39/CC-MAIN-2019-39_segments_1568514573070.36_warc_CC-MAIN-20190917101137-20190917123137-00344.warc.gz\"}"}