<?xml version="1.0" encoding="UTF-8"?>
<Worksheet>
<Version major="12" minor="2"/>
<Label-Scheme value="2" prefix=""/>
<View-Properties presentation="false"><Zoom percentage="150"/></View-Properties>
<MapleNet-Properties warnlevel="3" longdelim="true" plotoptions="" echo="1" errorcursor="false" elisiontermsthreshold="10000" elisiontermsbefore="100" errorbreak="1" prettyprint="3" prompt="> " latexwidth="6.0" useclientjvm="true" ansi="false" screenheight="25" labelwidth="20" plotdevice="inline" typesetting="standard" displayprecision="-1" indentamount="4" elisiondigitsthreshold="10000" quiet="false" plotoutput="terminal" rtablesize="10" preplot="" showassumed="1" plotdriver="openviz" elisiontermsafter="100" labelling="true" postplot="" screenwidth="79" elisiondigitsbefore="100" elisiondigitsafter="100" verboseproc="1" imaginaryunit="I" ShowLabels="true"/>
<Styles><Font name="Heading 4" background="[255,255,255]" bold="false" executable="false" family="Serif" foreground="[0,0,0]" italic="true" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Ordered List 5" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Ordered List 1" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Text Output" background="[255,255,255]" bold="false" executable="false" family="Monospaced" foreground="[0,0,255]" italic="false" opaque="false" readonly="true" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Bullet Item" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Equation Label" background="[255,255,255]" bold="true" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Author" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Text" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Diagnostic" background="[255,255,255]" bold="false" executable="false" family="Monospaced" foreground="[40,120,40]" italic="false" opaque="false" readonly="true" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Heading 2" background="[255,255,255]" bold="true" executable="false" family="Serif" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="16" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Ordered List 3" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Maple Plot" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Code" background="[255,255,255]" bold="false" executable="false" family="Monospaced" foreground="[255,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Header and Footer" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="10" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Normal" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Maple Input" background="[255,255,255]" bold="true" executable="true" family="Monospaced" foreground="[255,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Annotation Title" background="[255,255,255]" bold="true" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="18" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="2D Output" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,255]" italic="false" opaque="false" readonly="true" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="2D Inert Output" background="[255,255,255]" bold="false" executable="true" family="Times New Roman" foreground="[144,144,144]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Dictionary Hyperlink" background="[255,255,255]" bold="false" executable="false" family="Serif" foreground="[147,0,15]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="true" placeholder="false"/>
<Font name="Dash Item" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Heading 3" background="[255,255,255]" bold="true" executable="false" family="Serif" foreground="[0,0,0]" italic="true" opaque="false" readonly="false" size="14" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Ordered List 4" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Maple Output" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="List Item" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Line Printed Output" background="[255,255,255]" bold="false" executable="false" family="Monospaced" foreground="[0,0,255]" italic="false" opaque="false" readonly="true" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="2D Input" background="[255,255,255]" bold="false" executable="true" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Maple Input Placeholder" background="[255,255,255]" bold="true" executable="true" family="Monospaced" foreground="[200,0,200]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="true"/>
<Font name="2D Math" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Error" background="[255,255,255]" bold="false" executable="false" family="Monospaced" foreground="[255,0,255]" italic="false" opaque="false" readonly="true" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Warning" background="[255,255,255]" bold="false" executable="false" family="Monospaced" foreground="[0,0,255]" italic="false" opaque="false" readonly="true" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Heading 1" background="[255,255,255]" bold="true" executable="false" family="Serif" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="18" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Annotation Text" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Title" background="[255,255,255]" bold="true" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="18" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Ordered List 2" background="[255,255,255]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="false" placeholder="false"/>
<Font name="Hyperlink" background="[255,255,255]" bold="false" executable="false" family="Serif" foreground="[0,128,128]" italic="false" opaque="false" readonly="false" size="12" subscript="false" superscript="false" underline="true" placeholder="false"/>
<Layout name="Ordered List 5" alignment="left" bullet="ROMAN" firstindent="0" leftmargin="144" rightmargin="0" linespacing="0.0" spaceabove="3" spacebelow="3" linebreak="space" pagebreak-before="false" initial="-1" bulletsuffix=""/>
<Layout name="Heading 4" alignment="left" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="0" spacebelow="0" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Text Output" alignment="left" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="0" spacebelow="0" linebreak="newline" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Ordered List 1" alignment="left" bullet="numeric" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="3" spacebelow="3" linebreak="space" pagebreak-before="false" initial="-1" bulletsuffix=""/>
<Layout name="Annotation Title" alignment="centred" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="12" spacebelow="12" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Bullet Item" alignment="left" bullet="dot" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="3" spacebelow="3" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Author" alignment="centred" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="8" spacebelow="8" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Dash Item" alignment="left" bullet="dash" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="3" spacebelow="3" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Heading 3" alignment="left" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="0" spacebelow="0" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Diagnostic" alignment="left" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="0" spacebelow="0" linebreak="any" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Ordered List 4" alignment="left" bullet="ALPHABETIC" firstindent="0" leftmargin="108" rightmargin="0" linespacing="0.0" spaceabove="3" spacebelow="3" linebreak="space" pagebreak-before="false" initial="-1" bulletsuffix=""/>
<Layout name="Maple Output" alignment="centred" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.3" spaceabove="0" spacebelow="0" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Line Printed Output" alignment="left" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="0" spacebelow="0" linebreak="any" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="List Item" alignment="left" bullet="indent" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="3" spacebelow="3" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Heading 2" alignment="left" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="8" spacebelow="2" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Ordered List 3" alignment="left" bullet="roman" firstindent="0" leftmargin="72" rightmargin="0" linespacing="0.0" spaceabove="3" spacebelow="3" linebreak="space" pagebreak-before="false" initial="-1" bulletsuffix=""/>
<Layout name="Error" alignment="left" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="0" spacebelow="0" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Maple Plot" alignment="centred" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="0" spacebelow="0" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Warning" alignment="left" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="0" spacebelow="0" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Heading 1" alignment="left" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="8" spacebelow="4" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Title" alignment="centred" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="12" spacebelow="12" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Layout name="Ordered List 2" alignment="left" bullet="alphabetic" firstindent="0" leftmargin="36" rightmargin="0" linespacing="0.0" spaceabove="3" spacebelow="3" linebreak="space" pagebreak-before="false" initial="-1" bulletsuffix=""/>
<Layout name="Normal" alignment="left" bullet="none" firstindent="0" leftmargin="0" rightmargin="0" linespacing="0.0" spaceabove="0" spacebelow="0" linebreak="space" pagebreak-before="false" initial="0" bulletsuffix=""/>
<Pencil-style name="Pencil 3" pen-color="[0,0,0]" pen-height="3.0" pen-width="3.0" pen-opacity="1.0"/>
<Pencil-style name="Pencil 1" pen-color="[0,0,0]" pen-height="1.0" pen-width="1.0" pen-opacity="1.0"/>
<Pencil-style name="Pencil 5" pen-color="[255,0,0]" pen-height="5.0" pen-width="5.0" pen-opacity="1.0"/>
<Pencil-style name="Pencil 2" pen-color="[0,0,255]" pen-height="1.0" pen-width="1.0" pen-opacity="1.0"/>
<Pencil-style name="Pencil 4" pen-color="[0,0,255]" pen-height="3.0" pen-width="3.0" pen-opacity="1.0"/>
<Highlighter-style name="Highlighter 2" pen-color="[255,204,0]" pen-height="14.0" pen-width="14.0" pen-opacity="0.8"/>
<Highlighter-style name="Highlighter 4" pen-color="[0,255,255]" pen-height="32.0" pen-width="32.0" pen-opacity="0.8"/>
<Highlighter-style name="Highlighter 3" pen-color="[51,255,0]" pen-height="24.0" pen-width="24.0" pen-opacity="0.8"/>
<Highlighter-style name="Highlighter 1" pen-color="[255,153,255]" pen-height="12.0" pen-width="8.0" pen-opacity="0.8"/>
<Highlighter-style name="Highlighter 5" pen-color="[255,255,0]" pen-height="48.0" pen-width="48.0" pen-opacity="0.8"/>
</Styles>
<Task-table>
<Task-category name="<default>">
</Task-category>
</Task-table>
<Task>
</Task>
<Group labelreference="L376" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">restart; Digits:=14: with(LinearAlgebra):</Text-field>
</Input>
</Group>
<Group labelreference="L324" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">with(orthopoly):</Text-field>
</Input>
</Group>
<Group labelreference="L384" drawlabel="true">
<Input>
<Text-field style="Text" layout="Normal">Laguerre basis set for central potentials.</Text-field>
<Text-field style="Text" layout="Normal"></Text-field>
<Text-field style="Text" layout="Normal"></Text-field>
</Input>
</Group>
<Group labelreference="L383" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">L(2,1,x);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbW5HRiQ2JFEiM0YnLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGJy1JI21vR0YkNi1RKCZtaW51cztGJ0YvLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y4LyUpc3RyZXRjaHlHRjgvJSpzeW1tZXRyaWNHRjgvJShsYXJnZW9wR0Y4LyUubW92YWJsZWxpbWl0c0dGOC8lJ2FjY2VudEdGOC8lJ2xzcGFjZUdRLDAuMjIyMjIyMmVtRicvJSdyc3BhY2VHRkctRiM2JkYrLUYzNi1RMSZJbnZpc2libGVUaW1lcztGJ0YvRjZGOUY7Rj1GP0ZBRkMvRkZRJjAuMGVtRicvRklGUC1JI21pR0YkNiVRInhGJy8lJ2l0YWxpY0dRJXRydWVGJy9GMFEnaXRhbGljRidGLy1GMzYtUSIrRidGL0Y2RjlGO0Y9Rj9GQUZDRkVGSC1GIzYoLUkmbWZyYWNHRiQ2KC1GLDYkUSIxRidGLy1GLDYkUSIyRidGLy8lLmxpbmV0aGlja25lc3NHUSIxRicvJStkZW5vbWFsaWduR1EnY2VudGVyRicvJSludW1hbGlnbkdGaG8vJSliZXZlbGxlZEdGOEZMLUYjNiQtSSVtc3VwR0YkNiVGUkZgby8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRidGLy8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGOEYvRi8=">LCgiIiQiIiJJInhHNiIhIiQqJEYlIiIjI0YkRik=</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L317" drawlabel="true">
<Input>
<Text-field style="Normal" layout="Normal">Radial orbitals (multiplied by r, leaving out Y_Lm): look almost like hydrogenic orbitals, but form an enumerably infinite set.</Text-field>
</Input>
</Group>
<Group labelreference="L347" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">uL:=(n,l,beta)->r^(l+1)*exp(-beta*r/2)*L(n-1,2*l+2,beta*r);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEjdUxGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSM6PUYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZMLUYjNigtSShtZmVuY2VkR0YkNiQtRiM2Ki1GLDYlUSJuRidGL0YyLUY2Ni1RIixGJ0Y5RjsvRj9GMUZARkJGREZGRkgvRktRJjAuMGVtRicvRk5RLDAuMzMzMzMzM2VtRictRiw2JVEibEYnRi9GMkZZLUYsNiVRJyYjOTQ2O0YnL0YwRj1GOS8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGPUY5RjktRjY2LVEoJiM4NTk0O0YnRjlGO0Y+RkBGQkZERkZGSEZnbi9GTkZobi1GIzYoLUklbXN1cEdGJDYlLUYsNiVRInJGJ0YvRjItRiM2JkZbby1GNjYtUSIrRidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSwwLjIyMjIyMjJlbUYnL0ZORmlwLUkjbW5HRiQ2JFEiMUYnRjlGOS8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRictRjY2LVExJkludmlzaWJsZVRpbWVzO0YnRjlGO0Y+RkBGQkZERkZGSEZnbkZqby1GXnA2JS1GNjYtUS8mRXhwb25lbnRpYWxFO0YnRjlGO0Y+RkBGQkZERkZGSEZnbi9GTlEsMC4xMTExMTExZW1GJy1GIzYnLUY2Ni1RKiZ1bWludXMwO0YnRjlGO0Y+RkBGQkZERkZGSEZocEZqcC1GIzYoLUkmbWZyYWNHRiQ2KEZbcS1GXHE2JFEiMkYnRjkvJS5saW5ldGhpY2tuZXNzR1EiMUYnLyUrZGVub21hbGlnbkdRJ2NlbnRlckYnLyUpbnVtYWxpZ25HRl5zLyUpYmV2ZWxsZWRHRj1GYnEtRiM2JkZeb0ZicUZgcEY5RmJvRmVvRjlGYm9GZW9GOUZfcUZicS1GIzYoLUYsNiVRIkxGJ0YvRjItRjY2LVEwJkFwcGx5RnVuY3Rpb247RidGOUY7Rj5GQEZCRkRGRkZIRmduRmpvLUZSNiQtRiM2Ki1GIzYmRlYtRjY2LVEoJm1pbnVzO0YnRjlGO0Y+RkBGQkZERkZGSEZocEZqcEZbcUY5RlktRiM2Ji1GIzYmRmZyRmJxRltvRjlGZXBGZnJGOUZZRmNzRmJvRmVvRjlGOUZib0Zlb0Y5RjlGYm9GZW9GOUZib0Zlb0Y5">Zio2JUkibkc2IkkibEdGJUklYmV0YUdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKigpSSJyR0YlLCY5JSIiIkYwRjBGMC1JJGV4cEdGJTYjLCQqJjkmRjBGLUYwIyEiIiIiI0YwLUkiTEdGJTYlLCY5JEYwRjhGMCwmRi9GOUY5RjBGNUYwRiVGJUYl</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L356" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">plot([seq(uL(n,0,1.5),n=1..4)],r=0..20,color=[red,blue,green,black]);</Text-field>
<Text-field style="Maple Plot" layout="Maple Plot"></Text-field>
</Input>
<Output>
<Text-field style="Maple Plot" layout="Maple Plot"><Plot height="384" type="two-dimensional" width="716" plot-scale="1.0" plot-xtrans="0.0" plot-ytrans="0.0" gridlinevisibility="1" legendvisibility="false">NiktJSdDVVJWRVNHNiQ3ZW83JCQiIiEhIiIkIiIhISIiNyQkIjJLTEwzRldZcyMhIz0kIjJiYEAjXDtgcEUhIz03JCQiMmttbTthKUdcYSEjPSQiMjkncFhKdzFKXyEjPTckJCIsRCJHJFI8KSEjNyQiMVgtM2pZKXlvKCEjPDckJCIyTExMJDN4JikqMyIhIzwkIjJ0eGgoM3JKLzUhIzw3JCQiLERjJ3lNOyEjNiQiMiZmXXgocFpoVyIhIzw3JCQiMm1tbW1UOih6QCEjPCQiMmt5U2MiMyk0Jj0hIzw3JCQiMUxMJDNGV1lzIyEjOyQiMmR3LSkpNHI1QSMhIzw3JCQiK0RKZHBLISM1JCIxVVMieiQ+YmVEISM7NyQkIjFubTt6Pl05USEjOyQiMXIpPkhkVGEnRyEjOzckJCIxTExMTDNWZlYhIzskIjEyT1QlKilbTzkkISM7NyQkIjFtbSJ6PjV4SSYhIzskIjJ4Jz1eIlIzWmMkISM8NyQkIixEYyopZkQnISM2JCIyWSEpUXkpKlJKIlIhIzw3JCQiMUxMM0Yqb1U/KCEjOyQiMWJyRSM9SnA+JSEjOzckJCIxbm1tIkhbRDopISM7JCIwZCYqeTE6TFUlISM6NyQkIi12JHBVJkc1ISM2JCIyRVM5QjsiZmJaISM8NyQkIjJMTEwkZTAkPUMiISM7JCIxMFR4JkdtSCpbISM7NyQkIjJMTCQzeGZdJkgiISM7JCIwKGVIc3UvLlwhIzo3JCQiMk1MTGVSIj1cOCEjOyQiMllIZi1ZOlohXCEjPDckJCIyT0wkZTlvJkdTIiEjOyQiMVFZU3UlPScpKlshIzs3JCQiMk5MTExCS2xYIiEjOyQiMU1IZCRvbWApWyEjOzckJCIyT0xMMzIkKVFjIiEjOyQiMTxsWFs0bVJbISM7NyQkIjJPTEwkM1JCcjshIzskIjEvSjlVPXpyWiEjOzckJCIyLStdUDk6XCk9ISM7JCIyMnVyJT4yJ1xlJSEjPDckJCIybW1tInpqZik0IyEjOyQiMUJMP3cnRylbViEjOzckJCIsdkIxbkgjISM1JCIyYTk/RiNIQi1UISM8NyQkIjJNTExlNDtbXCMhIzskIjI4QiNccCNHMyVRISM8NyQkIjFubTt6dCUqKnAjISM6JCIyOVNzKGVcIVJjJCEjPDckJCIsRG15XSFIISM1JCIydyczbTAqZXlHJCEjPDckJCIxbm07SGRBPEohIzokIjJOdT4pZVUrNEkhIzw3JCQiMlBMTGV6cyRITCEjOyQiMUMqSExhWDV1IyEjOzckJCIybm1tVF1SM2EkISM7JCIyKVJlKEdILXdbIyEjPDckJCIsREAxQnYkISM1JCIyKWU3PUcyYVxBISM8NyQkIjFNTGU5ZCMpcFIhIzokIjEsNyssJilwQD8hIzs3JCQiMW5tbTtfTSg9JSEjOiQiMSUpZTlcIXk5Ij0hIzs3JCQiMU1MTDN5X3FYISM6JCIyViVHJDQzc0xbIiEjPDckJCIyMCsrK2wrPismISM7JCIyRWY8d2llWTwiISM8NyQkIip2V11WJiEiKSQiMDhfJXoiKVxCIyohIzs3JCQiKk5mQyZlISIpJCIxOio9PTZOQUUoISM8NyQkIjFNTCRlejY6QichIzokIjFcI1sxVWMiPmUhIzw3JCQiMW5tbTs9QyNvJyEjOiQiMVAhNCdmLTxdVyEjPDckJCIxbm1tbSNwUzEoISM6JCIyOi0kPj1BKEhgJCEjPTckJCIsRE9EIzN2ISM1JCIyOG1HaUZON3AjISM9NyQkIjFtbW1tKHk4IXohIzokIjFSR11aWCIqM0AhIzw3JCQiMSwrXWkudEskKSEjOiQiMjE1P2cqSEs0OyEjPTckJCIsdjN6TXUpISM1JCIyOzY8PycpUjRDIiEjPTckJCIxb21tIkhfPzwqISM6JCIxUFlEXWNEUiUqISM9NyQkIjFubTt6aWhsJiohIzokIjExWzsjUjgiR3QhIz03JCQiMUxMTDMjRywqKiohIzokIjFuSWBpWVdtYiEjPTckJCIyTUwkZXp3NVY1ISM6JCIxWzpKUEdedlQhIz03JCQiLXYkUSNcIjMiISM1JCIxJjMtcTI5aUMkISM9NyQkIjJOTExlIipbSDciISM6JCIxRyM+NyUqNCpwQyEjPTckJCIqZHhkOyIhIigkIjJqbnJKK20nZj0hIz43JCQiLEQweHc/IiEiKiQiMmlfcUMrLXFTIiEjPjckJCIyLCtdaSZwQFs3ISM6JCIyKSpSNWVxb0gyIiEjPjckJCIsdmdIS0giISIqJCIxVmMlMyI0YEp6ISM+NyQkIjJsbW1tWnZPTCIhIzokIjAwaF8lKmYkUmchIz03JCQiMiwrK10yZ29QIiEjOiQiMSVwJnkpKlsjKjRYISM+NyQkIjJLTCRlUjwqZlQiISM6JCIyYjguJmYpUiVlTSEjPzckJCIyLSsrXSlIeGU5ISM6JCIyajZBaVUmKlxlIyEjPzckJCIybW07SCFvLSpcIiEjOiQiMigpUmlYZzdUJz4hIz83JCQiMiwrXTdrLjZhIiEjOiQiMjhQJj5BSHhzOSEjPzckJCIybW1tO1dUQWUiISM6JCIyWkN1SipIbTU2ISM/NyQkIi1EMSozYGkiISM1JCIxMHZUIil6c2YjKSEjPzckJCIyTkxMTCp6eW07ISM6JCIxUXUlPTkiKmU/JyEjPzckJCIyT0xMM04xIzQ8ISM6JCIxPiVHOSkqKnBIWSEjPzckJCIxbm0iSFl0N3YiISM5JCIxKClRVSFRPSxZJCEjPzckJCIyLSsrK3hHKip5IiEjOiQiMnlCMz0lUlhZRSEjQDckJCIxbm07OUBCTT0hIzkkIjFuQSYzKCpmXyU+ISM/NyQkIjJPTExMYmRRKD0hIzokIjImZTQnR3ZkaloiISNANyQkIi1ET2w1Oz4hIzUkIjItMlhATm0nKjQiISNANyQkIjItK11QP1dsJj4hIzokIjExMWYlNCQ+IkgpISNANyQkIiQrIyEiIiQiMTtsLjVrLz1oISNALSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSdDVVJWRVNHNiQ3am83JCQiIiEhIiIkIiIhISIiNyQkIjJLTEwzRldZcyMhIz0kIjE2biZmMyNcKip5ISM8NyQkIjJrbW07YSlHXGEhIz0kIjJlIip6Uyo9Y0U6ISM8NyQkIixEIkckUjwpISM3JCIxdFknKT5dNTdBISM7NyQkIjJMTEwkM3gmKSozIiEjPCQiMjoyJypbInB3W0chIzw3JCQiLERjJ3lNOyEjNiQiMjtCbyhwOyNRKVIhIzw3JCQiMm1tbW1UOih6QCEjPCQiMFJeUnRdeCVcISM6NyQkIjFMTCQzRldZcyMhIzskIjFIJFEoKioqb2F2JiEjOzckJCIrREpkcEshIzUkIjFqU1o6KlwzVSchIzs3JCQiMW5tO3o+XTlRISM7JCIwQSc0TykqeWNwISM6NyQkIjFMTExMM1ZmViEjOyQiMSRvPWBwb19QKCEjOzckJCIxbW0iej41eEkmISM7JCIxLzs/NyRmZyZ5ISM7NyQkIixEYyopZkQnISM2JCIxdSNRYV9OdDEpISM7NyQkIjFtO3pXI0gsdCchIzskIjF3NGhuJT5tMykhIzs3JCQiMUxMM0Yqb1U/KCEjOyQiMTNWKDRUP2EwKSEjOzckJCIudiQ0JzMleXchIzgkIjF4c2RKWmh5eiEjOzckJCIxbm1tIkhbRDopISM7JCIxJFJYJ1s8dmd5ISM7NyQkIi12JHBVJkc1ISM2JCIxJlFfLW4heUhwISM7NyQkIjJMTEwkZTAkPUMiISM7JCIxbShvYUZZWGMmISM7NyQkIjJNTExlUiI9XDghIzskIjFQNF8vLDcpeSUhIzs3JCQiMk5MTExCS2xYIiEjOyQiMSNvLjgpKmVEKVIhIzs3JCQiMk9MTDMyJClRYyIhIzskIjFSSWQxcilmOyQhIzs3JCQiMk9MTCQzUkJyOyEjOyQiMmpTI0dKXT9gQiEjPDckJCIyb21UZ191IXk8ISM7JCIxRF9PeVYpKWY6ISM7NyQkIjItK11QOTpcKT0hIzskIjEsWG96SipbInohIzw3JCQiMk1MZTl3YjwqPiEjOyQiMVYmKT01TyYqSGIhIz03JCQiMm1tbSJ6amYpNCMhIzskITFPc007OW9KayEjPDckJCIsdkIxbkgjISM1JCEyTFF0UDlQZCM9ISM8NyQkIjJNTExlNDtbXCMhIzskITFtWWh0YXZdRyEjOzckJCIxbm07enQlKipwIyEjOiQhMVp4Qi4qPT11JCEjOzckJCIsRG15XSFIISM1JCEydlBPN1RjT1klISM8NyQkIjFubTtIZEE8SiEjOiQhMVUjKW9sYWdVXSEjOzckJCIyUExMZXpzJEhMISM7JCExLVg7bXEhZVkmISM7NyQkIjJubW1UXVIzYSQhIzskITE0IXlfI1FcXGQhIzs3JCQiLERAMUJ2JCEjNSQhMSJwJiopKkdFRyJmISM7NyQkIjFuO2FqZjFoUSEjOiQhMU0/JlFbb2MmZiEjOzckJCIxTUxlOWQjKXBSISM6JCEwNzU7YSFmdGYhIzo3JCQiLkRjWSZleVMhIzckITFDeGlpIT0pb2YhIzs3JCQiMW5tbTtfTSg9JSEjOiQhMXVRJSlvPFxWZiEjOzckJCIxLCtdN2wkKnlWISM6JCExSnchKiopZTFiZSEjOzckJCIxTUxMM3lfcVghIzokITFUWkB3UGQ+ZCEjOzckJCIyMCsrK2wrPismISM7JCEwL2dtdTckKkcmISM6NyQkIip2V11WJiEiKSQhMjBDJXBcJHBDdiUhIzw3JCQiKk5mQyZlISIpJCEyPHdqUCQpPm0+JSEjPDckJCIxTUwkZXo2OkInISM6JCExX2VPYlZkJHAkISM7NyQkIjFubW07PUMjbychIzokITFvVDcjKmZeREohIzs3JCQiMW5tbW0jcFMxKCEjOiQhMmp3P1hTI28kbyMhIzw3JCQiLERPRCMzdiEjNSQhMTZtdWcmKmVCQSEjOzckJCIxbW1tbSh5OCF6ISM6JCExYUBxSGIjbyc9ISM7NyQkIjEsK11pLnRLJCkhIzokITElKlsqKipHNihHOiEjOzckJCIsdjN6TXUpISM1JCEyJUhPRCx5QmI3ISM8NyQkIjFvbW0iSF8/PCohIzokITJcYXFuXyNbOjUhIzw3JCQiMW5tO3ppaGwmKiEjOiQhMSg+LGo0YWlKKSEjPDckJCIxTExMMyNHLCoqKiEjOiQhMSh6TXgoNFxybSEjPDckJCIyTUwkZXp3NVY1ISM6JCEyL0MvITQvaCFHJiEjPTckJCItdiRRI1wiMyIhIzUkITFoIVx0WnBBSCUhIzw3JCQiMk5MTGUiKltINyIhIzokITEpUSJvYDVTPk0hIzw3JCQiKmR4ZDsiISIoJCEyMjhLLmROU3AjISM9NyQkIixEMHh3PyIhIiokITJra21nVzBuNyMhIz03JCQiMiwrXWkmcEBbNyEjOiQhMic9VFYnKWYwKG8iISM9NyQkIix2Z0hLSCIhIiokITEmbyFvandrKzghIzw3JCQiMmxtbW1adk9MIiEjOiQhMiM0ZXFSNitGNSEjPTckJCIyLCsrXTJnb1AiISM6JCEwJlIlPiRlS2h6ISM8NyQkIjJLTCRlUjwqZlQiISM6JCExTUJdXjk6M2ohIz03JCQiMi0rK10pSHhlOSEjOiQhMVVhSjhPKTMpWyEjPTckJCIybW07SCFvLSpcIiEjOiQhMGA0ZlBfciNRISM8NyQkIjIsK103ay42YSIhIzokITFCZi4kcDZGJ0ghIz03JCQiMm1tbTtXVEFlIiEjOiQhMllaT2BpMUdJIyEjPjckJCItRDEqM2BpIiEjNSQhMmRdbiN5JioqZXciISM+NyQkIjJOTExMKnp5bTshIzokITFGJD5BKiozYU8iISM9NyQkIjJPTEwzTjEjNDwhIzokITIvcmpAKmUyWzUhIz43JCQiMW5tIkhZdDd2IiEjOSQhMWRFOkxaUV4hKSEjPjckJCIyLSsrK3hHKip5IiEjOiQhMU83JTMqKTQ6SichIz43JCQiMW5tOzlAQk09ISM5JCEyTGBhQFo0Jm9aISM/NyQkIjJPTExMYmRRKD0hIzokITE5RzJrKT1vcSQhIz43JCQiLURPbDU7PiEjNSQhMUNiIj4jcHJJRyEjPjckJCIyLStdUD9XbCY+ISM6JCEyS0ZyZTJ4WD0jISM/NyQkIiQrIyEiIiQhMiNmKTQyYHM9bCIhIz8tJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIjNSEiIi0lJ0NVUlZFU0c2JDdlcDckJCIiISEiIiQiIiEhIiI3JCQiMktMTDNGV1lzIyEjPSQiMmkyOlp6KyRlOiEjPDckJCIya21tO2EpR1xhISM9JCIxJmUqcD1YTnBIISM7NyQkIixEIkckUjwpISM3JCIxKFFdbSd5WVRVISM7NyQkIjJMTEwkM3gmKSozIiEjPCQiMD0qW0lgZSNRJiEjOjckJCIybG07YThBQk8iISM8JCIxUDBtRUZJK2shIzs3JCQiLERjJ3lNOyEjNiQiMSUpSHhuLSk9SSghIzs3JCQiMk1MJGUqKTREMj4hIzwkIjE4PiVHPmJVNCkhIzs3JCQiMm1tbW1UOih6QCEjPCQiMCg0KCpwUjAleSkhIzo3JCQiK0RKZHBLISM1JCIyYnJjMiEqeVIxIiEjOzckJCIxTExMTDNWZlYhIzskIjI9JiozbC5JNjgiISM7NyQkIjJrO3pXbitsZiUhIzwkIjF0T1h2WzNMNiEjOjckJCIuRGNecU4kWyEjOCQiMiMpPVBLYlI2OCIhIzs3JCQiMUszeGMua3FdISM7JCIyKVEpUTt1YGI3IiEjOzckJCIxbW0iej41eEkmISM7JCIxVERQQmZkOzYhIzo3JCQiMUwkMy0pKVw9eSYhIzskIjJHZ2kyZSVRKjMiISM7NyQkIixEYyopZkQnISM2JCIyR0hePThXODAiISM7NyQkIjFMTDNGKm9VPyghIzskIjE0PiQ+OFsxXCohIzs3JCQiMW5tbSJIW0Q6KSEjOyQiMT12eiwqPjBAKSEjOzckJCIxTUxlOXcpKj0jKiEjOyQiMTwrSnNdbXlsISM7NyQkIi12JHBVJkc1ISM2JCIxMTMhUlokUVhbISM7NyQkIjJsbVRnaSc9TjYhIzskIjI7Qm9CYC9VNCQhIzw3JCQiMkxMTCRlMCQ9QyIhIzskIjF0eiFSYlsjKlEiISM7NyQkIjJNTExlUiI9XDghIzskITJ0KipIaGBbIT1CISM9NyQkIjJOTExMQktsWCIhIzskITI6PWx3cic+QTwhIzw3JCQiMk9MTDMyJClRYyIhIzskITJgIm8zJikpW3owJCEjPDckJCIyT0xMJDNSQnI7ISM7JCEyJmVNOztBQ0NVISM8NyQkIjItK11QOTpcKT0hIzskITEzR1hOYWo8ZyEjOzckJCIybW1tInpqZik0IyEjOyQhMS4kcCU9MSYpPXIhIzs3JCQiMk9MTCQzOGwoPiMhIzskITEwSyEzKHB0NHUhIzs3JCQiLHZCMW5IIyEjNSQhMV4pMypHZSFHZCghIzs3JCQiMk5MJDMtUEJZQiEjOyQhMSlbZ2E1WiQ0dyEjOzckJCIybG1tbTtoZFIjISM7JCExeHJyZGxhPHchIzs3JCQiLURKJylHWEMhIzYkITE9PEx6OG8pZighIzs3JCQiMk1MTGU0O1tcIyEjOyQhMXczajIiW1NiKCEjOzckJCIxbm07enQlKipwIyEjOiQhMSI+PS80cEw3KCEjOzckJCIsRG15XSFIISM1JCExYDsjKVFZWmxqISM7NyQkIjFubTtIZEE8SiEjOiQhMVIucEZhIjNMJiEjOzckJCIyUExMZXpzJEhMISM7JCEwKFtRbFIheTclISM6NyQkIjJubW1UXVIzYSQhIzskITJFd3lMJj1oT0chIzw3JCQiLERAMUJ2JCEjNSQhMm9UYTdxRGpeIiEjPDckJCIxTUxlOWQjKXBSISM6JCExZlBxLDVSND0hIzw3JCQiMW5tbTtfTSg9JSEjOiQiMjsxJ2Z2QmsqMyIhIzw3JCQiMSwrXTdsJCp5ViEjOiQiMD1pbjhYLjgjISM6NyQkIjFNTEwzeV9xWCEjOiQiMWlsMSZvNUAzJCEjOzckJCIxbm07SFVAJ3klISM6JCIwbG1pIikqXE5TISM6NyQkIjIwKysrbCs+KyYhIzskIjJtaC02JUc+ZFshIzw3JCQiKnZXXVYmISIpJCIxdi1IIWY5eDUnISM7NyQkIiowX1BrJiEiKSQiMSYqZlo1WTdHbCEjOzckJCIqTmZDJmUhIikkIjEoUSg+KFFEJlJvISM7NyQkIjFubSJIZCYpPi8nISM6JCIxQyM9PmpwXy4oISM7NyQkIjFNTCRlejY6QichIzokIjEjXHY7eC9jOighIzs3JCQiMW47LyxWPldqISM6JCIwVyRvRGBZJT4oISM6NyQkIi1EMW8ob1gnISM2JCIwRCdbOiEpKjNAKCEjOjckJCIxTSRlOUpmJnBsISM6JCIxOjA9KFIiWzFzISM7NyQkIjFubW07PUMjbychIzokIjE8RW5FWngjPSghIzs3JCQiMW9tbVRiOnRvISM6JCIxOicqUW5VJ0c1KCEjOzckJCIxbm1tbSNwUzEoISM6JCIxbzVVX3A4enAhIzs3JCQiLERPRCMzdiEjNSQiMUYsYWUlXChlbCEjOzckJCIxbW1tbSh5OCF6ISM6JCIxeVNyJXBxJXpnISM7NyQkIjEsK11pLnRLJCkhIzokIjEqXFA9Y3gwXCYhIzs3JCQiLHYzek11KSEjNSQiMS4pKVE7TzcyXCEjOzckJCIxb21tIkhfPzwqISM6JCIyTnVIZigqW19JJSEjPDckJCIxbm07emlobCYqISM6JCIyJD5yVlU+RnhQISM8NyQkIjFMTEwzI0csKioqISM6JCIxTTpEXiQ0dEMkISM7NyQkIjJNTCRlenc1VjUhIzokIjJZRnohZScpUltGISM8NyQkIi12JFEjXCIzIiEjNSQiMTFXMiIpPnlmQiEjOzckJCIyTkxMZSIqW0g3IiEjOiQiMF8lKXBoZ3opPiEjOjckJCIqZHhkOyIhIigkIjIjNGYoeXohM2E7ISM8NyQkIixEMHh3PyIhIiokIjIpenQjPmopXHQ4ISM8NyQkIjIsK11pJnBAWzchIzokIjJyIj42KlErOjkiISM8NyQkIix2Z0hLSCIhIiokIjFNdDYlUXlZQyohIzw3JCQiMmxtbW1adk9MIiEjOiQiMUJVYzVOYjl3ISM8NyQkIjIsKytdMmdvUCIhIzokIjJhYyUzSU1DamghIz03JCQiMktMJGVSPCpmVCIhIzokIjImM3l2KkguLjImISM9NyQkIjItKytdKUh4ZTkhIzokIjJPNy5RLiEzIjMlISM9NyQkIjJtbTtIIW8tKlwiISM6JCIxQV5xNEtdO0whIzw3JCQiMiwrXTdrLjZhIiEjOiQiMiYqeiJcO1toaEUhIz03JCQiMm1tbTtXVEFlIiEjOiQiMiQ+QSRbWE4uOSMhIz03JCQiLUQxKjNgaSIhIzUkIjFwcGBaWnQpcCIhIzw3JCQiMk5MTEwqenltOyEjOiQiMkoheSozU0NpTiIhIz03JCQiMk9MTDNOMSM0PCEjOiQiMlUyb18mKnBYMiIhIz03JCQiMW5tIkhZdDd2IiEjOSQiMWJDPWVWUTUmKSEjPTckJCIyLSsrK3hHKip5IiEjOiQiMSR5eipwMEhibyEjPTckJCIxbm07OUBCTT0hIzkkIjFDKilIUyczJ1FgISM9NyQkIjJPTExMYmRRKD0hIzokIjF5KVxKInlvZ1UhIz03JCQiLURPbDU7PiEjNSQiMlohZV1yQHpWTCEjPjckJCIyLStdUD9XbCY+ISM6JCIxJTRySDomM1pFISM9NyQkIiQrIyEiIiQiMiRwQXgkZmpjMCMhIz4tJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiM1ISIiJCIiISEiIi0lJ0NVUlZFU0c2JDdecjckJCIiISEiIiQiIiEhIiI3JCQiMm1tO2E4QUJPIiEjPSQiMnglb2hWJGU1SyIhIzw3JCQiMktMTDNGV1lzIyEjPSQiMVJbYSpcUzpjIyEjOzckJCItRDFrJ3AzJSEjOCQiMTxCOSgqXEpDUCEjOzckJCIya21tO2EpR1xhISM9JCIxKip6VHFPPDdbISM7NyQkIixEIkckUjwpISM3JCIxO2gmKUd3JFJ4JyEjOzckJCIyTExMJDN4JikqMyIhIzwkIjFXTiYqeXNwbiUpISM7NyQkIjJsbTthOEFCTyIhIzwkIjEoKilIZ1EnPjgqKiEjOzckJCIsRGMneU07ISM2JCIyRyE0OGgyIkg2IiEjOzckJCIyTUwkZSopNEQyPiEjPCQiMi01dSFlITNMQCIhIzs3JCQiMm1tbW1UOih6QCEjPCQiMkE6UGUhKXlUSCIhIzs3JCQiMUxMJDNGV1lzIyEjOyQiMSlmNCpmeGEuOSEjOjckJCIrREpkcEshIzUkIjJadGg/Rm1EWCIhIzs3JCQiMk9MJDNfdi5VTiEjPCQiMmQoKWVqOCZwZDkhIzs3JCQiMW5tO3o+XTlRISM7JCIyJHBPNGNnW145ISM7NyQkIi1EMWsncDMlISM3JCIyayZ6QlRzLk45ISM7NyQkIjFMTExMM1ZmViEjOyQiMU08RSwsUTQ5ISM6NyQkIjFtbSJ6PjV4SSYhIzskIjIkXDJbdmhhaDchIzs3JCQiLERjKilmRCchIzYkIjIjW106djtqWzUhIzs3JCQiMW07elcjSCx0JyEjOyQiMWB3QXgiUXZFKiEjOzckJCIxTEwzRipvVT8oISM7JCIxJnlxOS1FYyl6ISM7NyQkIi52JDQnMyV5dyEjOCQiMUVOQWg2IWVtJyEjOzckJCIxbm1tIkhbRDopISM7JCIxaip5JFJXL0lgISM7NyQkIjEsXTdgencmbykhIzskIjEvOyJ5KT1SS1EhIzs3JCQiMU1MZTl3KSo9IyohIzskIjItTVVQPlc7TyMhIzw3JCQiMW47L3dzP18oKiEjOyQiMVItJ3AybS1QKiEjPDckJCItdiRwVSZHNSEjNiQhMk5NLDp0WWxEJSEjPTckJCIySyRlKilmWSc9MyIhIzskITFyM14wb2o4PCEjOzckJCIybG1UZ2knPU42ISM7JCExWEkjKmV1Jm8iSCEjOzckJCIvdj0jZjMmKT0iISM4JCEyTyRlImZMVXctJSEjPDckJCIyTExMJGUwJD1DIiEjOyQhMTY5X0cyVVNdISM7NyQkIjJNTExlUiI9XDghIzskITFkYzlMWm1vbiEjOzckJCIyTkxMTEJLbFgiISM7JCExaU4iUjg/VjIpISM7NyQkIjJPTEwzMiQpUWMiISM7JCExeCVwMydIUmwqKSEjOzckJCIyT0xMJDNSQnI7ISM7JCEwJj00bE1LayUqISM6NyQkIjE8L3dpUyV6cCIhIzokITFGNUV3TjpKJiohIzs3JCQiMixdKD08VWxDPCEjOyQhMU0jeSJbIWVpZCohIzs3JCQiMk5lOTtQazh2IiEjOyQhMCpwPkA4QisnKiEjOjckJCIyb21UZ191IXk8ISM7JCExQWtMNmxvLicqISM7NyQkIjItdm8vbyV5Lz0hIzskITAkZT5NR0QoZSohIzo3JCQiMk4kZSpbJFtcSj0hIzskITEmZiY0Ij10OmIqISM7NyQkIjE8SEsqKVw/ZT0hIzokITFmNWpTPUkoXCohIzs3JCQiMi0rXVA5OlwpPSEjOyQhMXQnKTNmQjVEJSohIzs3JCQiMk1MZTl3YjwqPiEjOyQhMSRROClwITQwKCopISM7NyQkIjJtbW0iempmKTQjISM7JCExJlFUXSZvSCRHKSEjOzckJCIsdkIxbkgjISM1JCExJFxTVykpMyRSbCEjOzckJCIyTUxMZTQ7W1wjISM7JCEyJylSeCJcWDctVyEjPDckJCIsdnQiUShmIyEjNSQhMU5PI1x6cXRAJCEjOzckJCIxbm07enQlKipwIyEjOiQhMic9KHArIz11Oj8hIzw3JCQiMk5MTDMtOEQhRyEjOyQhMSJma1U7LExAKSEjPDckJCIsRG15XSFIISM1JCIxTiMzTyoqKnBbTSEjPDckJCIxbm07SGRBPEohIzokIjEsNUhGQGUkZiMhIzs3JCQiMlBMTGV6cyRITCEjOyQiMmNrSU1HKnBGWCEjPDckJCIybm1tVF1SM2EkISM7JCIxQWQuSVdBcGchIzs3JCQiLERAMUJ2JCEjNSQiMURlNGVVYCE+KCEjOzckJCIxbjthamYxaFEhIzokIjE9UkkleXY4ZyghIzs3JCQiMU1MZTlkIylwUiEjOiQiMGQoUixYdSx6ISM6NyQkIjFuVDUhZjBVLSUhIzokIjE3SytOJz06LCkhIzs3JCQiLkRjWSZleVMhIzckIjE0dCdSJFEvJjQpISM7NyQkIjFNZTlUYCdIOCUhIzokIjEvaTZhIj5IOikhIzs3JCQiMW5tbTtfTSg9JSEjOiQiMCJ6XnQsIWU9KSEjOjckJCIxLF1pU0lDTlUhIzokIjEoNEd1PT5ZPikhIzs3JCQiMU1MZWszOSRHJSEjOiQiMWNGWTNfNSY9KSEjOzckJCIxbzthKW9RNUwlISM6JCIxWl1NLEZ6ZCIpISM7NyQkIjEsK103bCQqeVYhIzokIjFkNz4oeU9LNikhIzs3JCQiMW9tVGdAdHVXISM6JCIxWTEkbzItWih6ISM7NyQkIjFNTEwzeV9xWCEjOiQiMCdmYXBPRHV4ISM6NyQkIjFubTtIVUAneSUhIzokIjEiPS1gZScpPjcoISM7NyQkIjIwKysrbCs+KyYhIzskIjFUaWRWKm8jUWkhIzs3JCQiKUZaPV8hIigkIjF5Sl9DTjl2XiEjOzckJCIqdlddViYhIikkIjFwXlAleTMsKlIhIzs3JCQiKjBfUGsmISIpJCIyVzBKVm5ZInpGISM8NyQkIipOZkMmZSEiKSQiMmBRejNgai5hIiEjPDckJCIxbm0iSGQmKT4vJyEjOiQiMVdcPWAnZSMqPiUhIzw3JCQiMU1MJGV6NjpCJyEjOiQhMVVnVyVIMyVSbiEjPDckJCItRDFvKG9YJyEjNiQhMiMqPj9dOiVvOj4hIzw3JCQiMW5tbTs9QyNvJyEjOiQhMk03KlJua05ySSEjPDckJCIxb21tVGI6dG8hIzokITFqaCcpKSplOnFSISM7NyQkIjFubW1tI3BTMSghIzokITBbPDxMZHd5JSEjOjckJCIsRE9EIzN2ISM1JCExVnomeXBcUE4nISM7NyQkIjFtbW1tKHk4IXohIzokITEmKVJpLSFvTk0oISM7NyQkIjFNTGVrWDA8IikhIzokITEnSHhEJ1sqXHQoISM7NyQkIjEsK11pLnRLJCkhIzokITBmNlYiPTdFISkhIzo3JCQiK0RaNVEmKSEiKiQhMSg0Jio9REZsQCkhIzs3JCQiLHYzek11KSEjNSQhMUR4RGAzMUgkKSEjOzckJCIxTkxlKm9seCYqKSEjOiQhMTFsbl41OXIkKSEjOzckJCIxb21tIkhfPzwqISM6JCExKT5SPHNFV00pISM7NyQkIjFvbVQmR00pbyQqISM6JCExQGtgLmpPbSMpISM7NyQkIjFubTt6aWhsJiohIzokITEteiUpeiYpXFYiKSEjOzckJCIxTExMMyNHLCoqKiEjOiQhMShSd0E0SlJ2KCEjOzckJCIyTUwkZXp3NVY1ISM6JCExKilHJVI5bnRAKCEjOzckJCItdiRRI1wiMyIhIzUkITFIQFJZMSg+bychIzs3JCQiMk5MTGUiKltINyIhIzokITFCZk4lPSdSbmchIzs3JCQiKmR4ZDsiISIoJCExJ3lONydbc0FhISM7NyQkIixEMHh3PyIhIiokITJjNCl5SElLMFshIzw3JCQiMiwrXWkmcEBbNyEjOiQhMiU9JGUsK2ZkQiUhIzw3JCQiLHZnSEtIIiEiKiQhMUdraShlTXNrJCEjOzckJCIybG1tbVp2T0wiISM6JCExJnBWTys0UzskISM7NyQkIjIsKytdMmdvUCIhIzokITFqZCI0dz8oKXAjISM7NyQkIjJLTCRlUjwqZlQiISM6JCExb2cmKVFWZEFCISM7NyQkIjItKytdKUh4ZTkhIzokITJ3XS1hPmEkZj4hIzw3JCQiMm1tO0ghby0qXCIhIzokITIvInBQYVQpM20iISM8NyQkIjIsK103ay42YSIhIzokITIqKSpcKkcnZU8hUiIhIzw3JCQiMm1tbTtXVEFlIiEjOiQhMnhtJ1FMIVtKOyIhIzw3JCQiLUQxKjNgaSIhIzUkITF2Kj1ZKWVuMCcqISM8NyQkIjJOTExMKnp5bTshIzokITF2Kj5TZzloJnohIzw3JCQiMk9MTDNOMSM0PCEjOiQhMTwvKT1oTmlgJyEjPDckJCIxbm0iSFl0N3YiISM5JCExLyVvaCpwKillYCEjPDckJCIyLSsrK3hHKip5IiEjOiQhMSI0JjNXJz46WCUhIzw3JCQiMW5tOzlAQk09ISM5JCExKlsuM0Jvb2UkISM8NyQkIjJPTExMYmRRKD0hIzokITIlXHJGemlQW0ghIz03JCQiLURPbDU7PiEjNSQhMm1PKT5fKnBiUSMhIz03JCQiMi0rXVA/V2wmPiEjOiQhMmNFTShmJyl5VT4hIz03JCQiJCsjISIiJCExdiNcIil5JClSYiIhIzwtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lJVZJRVdHNiQ7JCIiISEiIiQiJCsjISIiJShERUZBVUxURy0lK0FYRVNMQUJFTFNHNictSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USJyNiIvJSdmYW1pbHlHUSE2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiJRITYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM1JStIT1JJWk9OVEFMRyUrSE9SSVpPTlRBTEctJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiRxJiEiIi0lKUJPVU5EU19ZRzYjJCIjXSEiIi0lLUJPVU5EU19XSURUSEc2IyQiJXFrISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJT9QISIiLSUpQ0hJTERSRU5HNiI=</Plot></Text-field>
</Output>
</Group>
<Group labelreference="L377" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">int(uL(1,0,1.5)^2,r=0..infinity);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">JCIvZiNmI2YjZiNmISM5</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L330" drawlabel="true">
<Input>
<Text-field style="Normal" layout="Normal">This basis is orthogonal, but not yet normalized. Orthogonality is demonstrated below.</Text-field>
</Input>
</Group>
<Group labelreference="L326" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">int(uL(1,0,1.5)*uL(2,0,1.5),r=0..infinity);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">JCIiIUYj</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L341" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># Define a basis with fixed parameter beta:</Text-field>
</Input>
</Group>
<Group labelreference="L358" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">phi:=n->uL(n,0,3/2);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEnJiM5NjY7RicvJSdpdGFsaWNHUSZmYWxzZUYnLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGJy1JI21vR0YkNi1RIzo9RidGMi8lJmZlbmNlR0YxLyUqc2VwYXJhdG9yR0YxLyUpc3RyZXRjaHlHRjEvJSpzeW1tZXRyaWNHRjEvJShsYXJnZW9wR0YxLyUubW92YWJsZWxpbWl0c0dGMS8lJ2FjY2VudEdGMS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkktRiM2KC1GLDYlUSJuRicvRjBRJXRydWVGJy9GM1EnaXRhbGljRictRjY2LVEoJiM4NTk0O0YnRjJGOUY7Rj1GP0ZBRkNGRS9GSFEmMC4wZW1GJy9GS0ZZLUYjNigtRiw2JVEjdUxGJ0ZRRlMtRjY2LVEwJkFwcGx5RnVuY3Rpb247RidGMkY5RjtGPUY/RkFGQ0ZFRlhGWi1JKG1mZW5jZWRHRiQ2JC1GIzYqRk4tRjY2LVEiLEYnRjJGOS9GPEZSRj1GP0ZBRkNGRUZYL0ZLUSwwLjMzMzMzMzNlbUYnLUkjbW5HRiQ2JFEiMEYnRjJGYm8tSSZtZnJhY0dGJDYoLUZpbzYkUSIzRidGMi1GaW82JFEiMkYnRjIvJS5saW5ldGhpY2tuZXNzR1EiMUYnLyUrZGVub21hbGlnbkdRJ2NlbnRlckYnLyUpbnVtYWxpZ25HRmpwLyUpYmV2ZWxsZWRHRjEvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRjFGMkYyRl9xRmJxRjJGX3FGYnFGMkZfcUZicUYy">Zio2I0kibkc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkjdUxHRiU2JTkkIiIhIyIiJCIiI0YlRiVGJQ==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L337" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">OL:=(m,n)->int(phi(n)*phi(m),r=0..infinity);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEjT0xGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSM6PUYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZMLUYjNigtSShtZmVuY2VkR0YkNiQtRiM2KC1GLDYlUSJtRidGL0YyLUY2Ni1RIixGJ0Y5RjsvRj9GMUZARkJGREZGRkgvRktRJjAuMGVtRicvRk5RLDAuMzMzMzMzM2VtRictRiw2JVEibkYnRi9GMi8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGPUY5RjktRjY2LVEoJiM4NTk0O0YnRjlGO0Y+RkBGQkZERkZGSEZnbi9GTkZobi1GIzYqLUkobXN1YnN1cEdGJDYnLUY2Ni1RKCYjODc0NztGJ0Y5RjtGPkZARkJGREZGRkhGZ25GZm8tSSNtbkdGJDYkUSIwRidGOS1GLDYlUSgmIzg3MzQ7RidGL0YyLyUxc3VwZXJzY3JpcHRzaGlmdEdRIjJGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnLUYjNiYtRiM2KC1GLDYlUScmIzk2NjtGJy9GMEY9RjktRjY2LVEwJkFwcGx5RnVuY3Rpb247RidGOUY7Rj5GQEZCRkRGRkZIRmduRmZvLUZSNiQtRiM2JkZbb0Zeb0Zhb0Y5RjlGXm9GYW9GOS1GNjYtUTEmSW52aXNpYmxlVGltZXM7RidGOUY7Rj5GQEZCRkRGRkZIRmduRmZvLUYjNihGYHFGZHEtRlI2JC1GIzYmRlZGXm9GYW9GOUY5Rl5vRmFvRjlGOS1JJ21zcGFjZUdGJDYmLyUnaGVpZ2h0R1EmMC4wZXhGJy8lJndpZHRoR1EmMC4zZW1GJy8lJmRlcHRoR0Zpci8lKmxpbmVicmVha0dRJWF1dG9GJy1GNjYtUTAmRGlmZmVyZW50aWFsRDtGJ0Y5RjtGPkZARkJGREZGRkhGZ25GZm8tRiw2JVEickYnRi9GMkZeb0Zhb0Y5Rl5vRmFvRjlGXm9GYW9GOQ==">Zio2JEkibUc2IkkibkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkkaW50RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiQqJi1JJHBoaUdGJTYjOSUiIiItRjI2IzkkRjUvSSJyR0YlOyIiIUkpaW5maW5pdHlHRi1GJUYlRiU=</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L351" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">OL(2,2);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">IyIjOyIiKg==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L423" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># a few orthogonality checks:</Text-field>
</Input>
</Group>
<Group labelreference="L349" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">OL(1,2);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">IiIh</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L344" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">OL(2,4);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">IiIh</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L370" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># Now define the exponential potential and its matrix elements:</Text-field>
</Input>
</Group>
<Group labelreference="L375" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">V0:=-2;</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEjVjBGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSM6PUYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZMLUYjNiUtRjY2LVEqJnVtaW51czA7RidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSwwLjIyMjIyMjJlbUYnL0ZORlUtSSNtbkdGJDYkUSIyRidGOUY5LyUrZm9yZWdyb3VuZEdRKFswLDAsMF1GJy8lKXJlYWRvbmx5R0Y9Rjk=">ISIj</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L388" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">r0:=4.0; # to be consistent with the potential used in PH5000L5.mw where phase shifts are calculated.</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEjcjBGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSM6PUYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZMLUkjbW5HRiQ2JFEkNC4wRidGOS8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGPUY5">JCIjUyEiIg==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L348" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">V:=V0*exp(-r/r0);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEiVkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RIzo9RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y9LyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkwtRiM2Jy1GNjYtUSomdW1pbnVzMDtGJ0Y5RjtGPkZARkJGREZGRkgvRktRLDAuMjIyMjIyMmVtRicvRk5GVS1GIzYmLUkjbW5HRiQ2JFEiMkYnRjktRjY2LVExJkludmlzaWJsZVRpbWVzO0YnRjlGO0Y+RkBGQkZERkZGSC9GS1EmMC4wZW1GJy9GTkZbby1JJW1zdXBHRiQ2JS1GNjYtUS8mRXhwb25lbnRpYWxFO0YnRjlGO0Y+RkBGQkZERkZGSEZqbi9GTlEsMC4xMTExMTExZW1GJy1GIzYnRlEtRiM2Ji1GWjYkUTEwLjI1MDAwMDAwMDAwMDAwRidGOUZnbi1GLDYlUSJyRidGL0YyRjkvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRj1GOS8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRidGOUZfcEZicEY5Rl9wRmJwRjk=">LCQtSSRleHBHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHNiI2IywkSSJyR0YoJCEvKysrKysrRCEjOSEiIw==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L352" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">VME:=(m,n)->int(phi(n)*V*phi(m),r=0..infinity);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEkVk1FRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LVEjOj1GJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRj0vJSlzdHJldGNoeUdGPS8lKnN5bW1ldHJpY0dGPS8lKGxhcmdlb3BHRj0vJS5tb3ZhYmxlbGltaXRzR0Y9LyUnYWNjZW50R0Y9LyUnbHNwYWNlR1EsMC4yNzc3Nzc4ZW1GJy8lJ3JzcGFjZUdGTC1GIzYoLUkobWZlbmNlZEdGJDYkLUYjNigtRiw2JVEibUYnRi9GMi1GNjYtUSIsRidGOUY7L0Y/RjFGQEZCRkRGRkZIL0ZLUSYwLjBlbUYnL0ZOUSwwLjMzMzMzMzNlbUYnLUYsNiVRIm5GJ0YvRjIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRj1GOUY5LUY2Ni1RKCYjODU5NDtGJ0Y5RjtGPkZARkJGREZGRkhGZ24vRk5GaG4tRiM2Ki1JKG1zdWJzdXBHRiQ2Jy1GNjYtUSgmIzg3NDc7RidGOUY7Rj5GQEZCRkRGRkZIRmduRmZvLUkjbW5HRiQ2JFEiMEYnRjktRiw2JVEoJiM4NzM0O0YnRi9GMi8lMXN1cGVyc2NyaXB0c2hpZnRHUSIyRicvJS9zdWJzY3JpcHRzaGlmdEdRIjBGJy1GIzYoLUYjNigtRiw2JVEnJiM5NjY7RicvRjBGPUY5LUY2Ni1RMCZBcHBseUZ1bmN0aW9uO0YnRjlGO0Y+RkBGQkZERkZGSEZnbkZmby1GUjYkLUYjNiZGW29GXm9GYW9GOUY5Rl5vRmFvRjktRjY2LVExJkludmlzaWJsZVRpbWVzO0YnRjlGO0Y+RkBGQkZERkZGSEZnbkZmby1GLDYlUSJWRidGL0YyRltyLUYjNihGYHFGZHEtRlI2JC1GIzYmRlZGXm9GYW9GOUY5Rl5vRmFvRjlGOS1JJ21zcGFjZUdGJDYmLyUnaGVpZ2h0R1EmMC4wZXhGJy8lJndpZHRoR1EmMC4zZW1GJy8lJmRlcHRoR0Zccy8lKmxpbmVicmVha0dRJWF1dG9GJy1GNjYtUTAmRGlmZmVyZW50aWFsRDtGJ0Y5RjtGPkZARkJGREZGRkhGZ25GZm8tRiw2JVEickYnRi9GMkZeb0Zhb0Y5Rl5vRmFvRjlGXm9GYW9GOQ==">Zio2JEkibUc2IkkibkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkkaW50RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiQqKC1JJHBoaUdGJTYjOSUiIiJJIlZHRiVGNS1GMjYjOSRGNS9JInJHRiU7IiIhSSlpbmZpbml0eUdGLUYlRiVGJQ==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L345" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">VME(1,1);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">JCEvP0peb2JqdSEjOQ==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L366" drawlabel="true">
<Input>
<Text-field style="Normal" layout="Normal">Had we defined the basis with a non-float (using 15/10 instead of 1.5) the numbers would be computed symbolically exactly.</Text-field>
</Input>
</Group>
<Group labelreference="L329" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">VME(2,1);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">JCEvMyo+QW4nKT4kISM5</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L333" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">VME(3,1);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">JCEvIilvMCNcIVIiKiEjOg==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L355" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">VME(1,6);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">JCEvdidSVGdiSyohIzw=</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L381" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">TME:=(m,n)->-1/2*int(phi(n)*diff(phi(m),r$2),r=0..infinity);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEkVE1FRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LVEjOj1GJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRj0vJSlzdHJldGNoeUdGPS8lKnN5bW1ldHJpY0dGPS8lKGxhcmdlb3BHRj0vJS5tb3ZhYmxlbGltaXRzR0Y9LyUnYWNjZW50R0Y9LyUnbHNwYWNlR1EsMC4yNzc3Nzc4ZW1GJy8lJ3JzcGFjZUdGTC1GIzYoLUkobWZlbmNlZEdGJDYkLUYjNigtRiw2JVEibUYnRi9GMi1GNjYtUSIsRidGOUY7L0Y/RjFGQEZCRkRGRkZIL0ZLUSYwLjBlbUYnL0ZOUSwwLjMzMzMzMzNlbUYnLUYsNiVRIm5GJ0YvRjIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRj1GOUY5LUY2Ni1RKCYjODU5NDtGJ0Y5RjtGPkZARkJGREZGRkhGZ24vRk5GaG4tRiM2Jy1GNjYtUSomdW1pbnVzMDtGJ0Y5RjtGPkZARkJGREZGRkgvRktRLDAuMjIyMjIyMmVtRicvRk5GXXAtRiM2KC1JJm1mcmFjR0YkNigtSSNtbkdGJDYkUSIxRidGOS1GZXA2JFEiMkYnRjkvJS5saW5ldGhpY2tuZXNzR1EiMUYnLyUrZGVub21hbGlnbkdRJ2NlbnRlckYnLyUpbnVtYWxpZ25HRmBxLyUpYmV2ZWxsZWRHRj0tRjY2LVExJkludmlzaWJsZVRpbWVzO0YnRjlGO0Y+RkBGQkZERkZGSEZnbkZmby1GIzYkLUYjNiotSShtc3Vic3VwR0YkNictRjY2LVEoJiM4NzQ3O0YnRjlGO0Y+RkBGQkZERkZGSEZnbkZmby1GZXA2JFEiMEYnRjktRiw2JVEoJiM4NzM0O0YnRi9GMi8lMXN1cGVyc2NyaXB0c2hpZnRHUSIyRicvJS9zdWJzY3JpcHRzaGlmdEdRIjBGJy1GIzYmLUYjNigtRiw2JVEnJiM5NjY7RicvRjBGPUY5LUY2Ni1RMCZBcHBseUZ1bmN0aW9uO0YnRjlGO0Y+RkBGQkZERkZGSEZnbkZmby1GUjYkLUYjNiZGW29GXm9GYW9GOUY5Rl5vRmFvRjlGZXEtRlI2JC1GIzYoLUZicDYoLUYjNiYtSSVtc3VwR0YkNiUtRjY2LVErJlBhcnRpYWxEO0YnRjlGO0Y+RkBGQkZERkZGSEZnbkZmb0ZocC9GaXJGXXNGXm9GYW9GOS1GIzYmLUYjNidGaHQtRmZ0NiUtRiw2JVEickYnRi9GMkZocEZbdUZeb0Zhb0Y5Rl5vRmFvRjlGW3FGXnFGYXFGY3FGZXEtRiM2KEZic0Zmcy1GUjYkLUYjNiZGVkZeb0Zhb0Y5RjlGXm9GYW9GOUZeb0Zhb0Y5RjlGOS1JJ21zcGFjZUdGJDYmLyUnaGVpZ2h0R1EmMC4wZXhGJy8lJndpZHRoR1EmMC4zZW1GJy8lJmRlcHRoR0Zgdi8lKmxpbmVicmVha0dRJWF1dG9GJy1GNjYtUTAmRGlmZmVyZW50aWFsRDtGJ0Y5RjtGPkZARkJGREZGRkhGZ25GZm9GYnVGXm9GYW9GOUY5Rl5vRmFvRjlGXm9GYW9GOUZeb0Zhb0Y5Rl5vRmFvRjk=">Zio2JEkibUc2IkkibkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCQtSSRpbnRHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiU2JComLUkkcGhpR0YlNiM5JSIiIi1JJWRpZmZHRi42JC1GMzYjOSQtSSIkR0YuNiRJInJHRiUiIiNGNi9GQDsiIiFJKWluZmluaXR5R0YuIyEiIkZBRiVGJUYl</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L332" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">TME(1,1);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">IyIiIiIiJw==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L342" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">TME(1,3);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">IyIiIiIiJA==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L321" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">TME(1,2),TME(2,1);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">NiQjIiIiIiIkRiM=</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L364" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">TME(10,10);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">IyIkOigiIic=</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L338" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">TME(10,11);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">IyIkJlEiIiQ=</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L325" drawlabel="true">
<Input>
<Text-field style="Normal" layout="Normal">The kinetic energy is not diagonal, but symmetric!</Text-field>
<Text-field style="Normal" layout="Normal"></Text-field>
</Input>
</Group>
<Group labelreference="L365" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Nb:=10;</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEjTmJGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSM6PUYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZMLUkjbW5HRiQ2JFEjMTBGJ0Y5LyUrZm9yZWdyb3VuZEdRKFswLDAsMF1GJy8lKXJlYWRvbmx5R0Y9Rjk=">IiM1</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L369" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">VM:=Matrix(Nb,Nb,shape=symmetric):</Text-field>
</Input>
</Group>
<Group labelreference="L319" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">HM:=Matrix(Nb,Nb,shape=symmetric):</Text-field>
</Input>
</Group>
<Group labelreference="L340" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">OM:=Matrix(Nb,Nb,shape=symmetric):</Text-field>
<Text-field prompt="> " style="Maple Input" layout="Normal">for m from 1 to Nb do: for n from 1 to m do: VM[m,n]:=VME(m,n): OM[m,n]:=OL(m,n): od: od:</Text-field>
<Text-field prompt="> " style="Maple Input" layout="Normal">#print(map(evalf,VM));</Text-field>
<Text-field prompt="> " style="Maple Input" layout="Normal">for m from 1 to Nb do: for n from 1 to m do: HM[m,n]:=evalf(VM[m,n])+evalf(TME(m,n)): od: od:</Text-field>
</Input>
</Group>
<Group labelreference="L382" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"></Text-field>
</Input>
</Group>
<Group labelreference="L350" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">E10:=(convert(sort(Re(Eigenvalues(HM,OM))),list));</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEkRTEwRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LVEjOj1GJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRj0vJSlzdHJldGNoeUdGPS8lKnN5bW1ldHJpY0dGPS8lKGxhcmdlb3BHRj0vJS5tb3ZhYmxlbGltaXRzR0Y9LyUnYWNjZW50R0Y9LyUnbHNwYWNlR1EsMC4yNzc3Nzc4ZW1GJy8lJ3JzcGFjZUdGTC1JKG1mZW5jZWRHRiQ2Ji1GIzY4LUYjNiUtRjY2LVEqJnVtaW51czA7RidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSwwLjIyMjIyMjJlbUYnL0ZORlotSSNtbkdGJDYkUTQxLjAwOTg3NzcyNzUwNDA5OTk5RidGOUY5LUY2Ni1RIixGJ0Y5RjsvRj9GMUZARkJGREZGRkgvRktRJjAuMGVtRicvRk5RLDAuMzMzMzMzM2VtRictRiM2JUZWLUZnbjYkUTUwLjQ5NjA5MjMzNjMzNzczMTk3NkYnRjlGOUZqbi1GIzYlRlYtRmduNiRRNTAuMjEzNTM3MzI1MTQ0Mjg1MDAxRidGOUY5RmpuLUYjNiVGVi1GZ242JFE2MC4wNjMwNDg5OTA5Njg1NjkwOTkxRidGOUY5RmpuLUYjNiVGVi1GZ242JFE4MC4wMDA4ODY5NjE1MjkyNzIxMTMwMjlGJ0Y5RjlGam4tRmduNiRRNjAuMDQ1MjkyNDYxNzgzNTQ0OTk4OUYnRjlGam4tRmduNiRRNTAuMTM1MTMyMzQwMTI3Mzc3OTg3RidGOUZqbi1GZ242JFE1MC40MjQzMTc5Mjc3MjExMTA5ODhGJ0Y5RmpuLUZnbjYkUTQxLjkyNjI4MTM0NDk3NDE4MDAyRidGOUZqbi1GZ242JFE0MTEuOTc4OTE3ODQ4OTcwOTk5M0YnRjkvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRj1GOUY5LyUlb3BlbkdRIltGJy8lJmNsb3NlR1EiXUYnRmVxRmhxRjk=">NywkITMqKio0L3ZzeCk0NSEjPCQhM3c+dFBqTCM0J1whIz0kITMsXUdXXktQTkBGKCQhMyIqNHAmbzQqKltJJyEjPiQhM0hJNnMjSDoncCkpISNAJCIzKikqXGEkeWhDSFhGLSQiMygpelBGLE1LXjhGKCQiMykpNDZAeCN6SkMlRigkIjMtIT11XE0iR0U+RiUkIjMkKio0KCpbeSIqeT4iISM7</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L449" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Nb:=20;</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEjTmJGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSM6PUYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZMLUkjbW5HRiQ2JFEjMjBGJ0Y5LyUrZm9yZWdyb3VuZEdRKFswLDAsMF1GJy8lKXJlYWRvbmx5R0Y9Rjk=">IiM/</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L450" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">VM:=Matrix(Nb,Nb,shape=symmetric):</Text-field>
</Input>
</Group>
<Group labelreference="L447" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">HM:=Matrix(Nb,Nb,shape=symmetric):</Text-field>
</Input>
</Group>
<Group labelreference="L446" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">OM:=Matrix(Nb,Nb,shape=symmetric):</Text-field>
<Text-field prompt="> " style="Maple Input" layout="Normal">for m from 1 to Nb do: for n from 1 to m do: VM[m,n]:=VME(m,n): OM[m,n]:=OL(m,n): od: od:</Text-field>
<Text-field prompt="> " style="Maple Input" layout="Normal">#print(map(evalf,VM));</Text-field>
<Text-field prompt="> " style="Maple Input" layout="Normal">for m from 1 to Nb do: for n from 1 to m do: HM[m,n]:=evalf(VM[m,n])+evalf(TME(m,n)): od: od:</Text-field>
</Input>
</Group>
<Group labelreference="L448" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"></Text-field>
</Input>
</Group>
<Group labelreference="L451" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">E20:=(convert(sort(Re(Eigenvalues(HM,OM))),list));</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEkRTIwRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LVEjOj1GJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRj0vJSlzdHJldGNoeUdGPS8lKnN5bW1ldHJpY0dGPS8lKGxhcmdlb3BHRj0vJS5tb3ZhYmxlbGltaXRzR0Y9LyUnYWNjZW50R0Y9LyUnbHNwYWNlR1EsMC4yNzc3Nzc4ZW1GJy8lJ3JzcGFjZUdGTC1JKG1mZW5jZWRHRiQ2Ji1GIzZMLUYjNiUtRjY2LVEqJnVtaW51czA7RidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSwwLjIyMjIyMjJlbUYnL0ZORlotSSNtbkdGJDYkUTQxLjAwOTg3ODQ4Mzk5MTAzOTk3RidGOUY5LUY2Ni1RIixGJ0Y5RjsvRj9GMUZARkJGREZGRkgvRktRJjAuMGVtRicvRk5RLDAuMzMzMzMzM2VtRictRiM2JUZWLUZnbjYkUTUwLjQ5NjEwNDc3NTkxMzYzNTk5OEYnRjlGOUZqbi1GIzYlRlYtRmduNiRRNTAuMjEzNTcxNzgzNzM0MjIzMDAzRidGOUY5RmpuLUYjNiVGVi1GZ242JFE2MC4wNjMwODcwMzg4MDk5NjI0MDA4RidGOUY5RmpuLUYjNiVGVi1GZ242JFE3MC4wMDM3MDM5NDUwNDk3NDY2NTAxMEYnRjlGOUZqbi1GZ242JFE3MC4wMDQ5Nzg3NDAxNDg4MzQ4NTk5NkYnRjlGam4tRmduNiRRNjAuMDE2NzcyNjAzNTQ5NDcyNTk5OUYnRjlGam4tRmduNiRRNjAuMDM0OTY0NzE5NjAwNDQ4NzAyOUYnRjlGam4tRmduNiRRNjAuMDYwNjkxNjc0MDY0NjMyMDk5M0YnRjlGam4tRmduNiRRNjAuMDk2MDAyMjQyOTQ1OTAwMjA2OUYnRjlGam4tRmduNiRRNTAuMTQ0MTgxODU4NjU3NDkxMDAyRidGOUZqbi1GZ242JFE1MC4yMTA2MTE0NjU0MDg0MDAwMDZGJ0Y5RmpuLUZnbjYkUTUwLjMwNDI3MjMwOTk4MzI2MDk5MUYnRjlGam4tRmduNiRRNTAuNDQxODI4MjkxODQ1ODg2MDEyRidGOUZqbi1GZ242JFE1MC42NTc5OTYyNTYyODU2NjMwNTBGJ0Y5RmpuLUZnbjYkUTQxLjA0OTcwMTIwNTQ3NTM5MDAzRidGOUZqbi1GZ242JFE0MS45MDUxMjA3NDU1MDE1Njk5NkYnRjlGam4tRmduNiRRNDQuMDY2MjE1NjEzMDQzODgwMjRGJ0Y5RmpuLUZnbjYkUTQxMC44MDQxMjAxMDUzMjk1MDA4RidGOUZqbi1GZ242JFE0NDguMzM3ODA3MzYwNzc5NTAwMkYnRjkvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRj1GOUY5LyUlb3BlbkdRIltGJy8lJmNsb3NlR1EiXUYnRmNzRmZzRjk=">NzYkITMoKlI1KlJbeSk0NSEjPCQhMykqZmo4ZngvaFwhIz0kITMuSUFNUHlyTkBGKCQhMzNTaSo0KVFxM2ohIz4kITM1XW11XF0lUnEkISM/JCIzJypmWyQpWyx1eVxGMCQiMyoqZnMlXE5nc24iRi0kIjNIcVsvZz5aJ1wkRi0kIjMkKjRLWTF1O3BnRi0kIjNwPytmJUhDLWcqRi0kIjMtNVxkJ2U9PVciRigkIjMxK1MzYVk2MUBGKCQiMyIqNEUkKSo0QkYvJEYoJCIzN2cpZSU9SEc9V0YoJCIzXUltJkdjaSp6bEYoJCIzLiFSdmE/LChcNUYlJCIzJypwOl1YMjcwPkYlJCIzQyEpUS84Y0BtU0YlJCIzM11IYDU/VCEzIiEjOyQiMy1dejJPMnlMW0ZN</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L444" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># We see that the first 4 negative-E eigenvalues are confirmed in their lower significant digits.</Text-field>
</Input>
</Group>
<Group labelreference="L445" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># Problem 10.8 in Quantum Mechanics by Ghatak and Lokanathan, QC174.12 G485 2004 asks to transform the SE for the s-states of an exp. potential to the Bessel equation.</Text-field>
</Input>
</Group>
<Group labelreference="L362" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># The exact s-state eigenvalues are obtained from:</Text-field>
</Input>
</Group>
<Group labelreference="L389" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">g:=sqrt(r0^2*abs(8*V0));</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEiZ0YnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RIzo9RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y9LyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkwtSSNtbkdGJDYkUTAxNi4wMDAwMDAwMDAwMDBGJ0Y5LyUrZm9yZWdyb3VuZEdRKFswLDAsMF1GJy8lKXJlYWRvbmx5R0Y9Rjk=">JCIvKysrKysrOyEjNw==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L431" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">plot(BesselJ(nu,g),nu=0..30);</Text-field>
</Input>
<Output>
<Text-field style="Maple Plot" layout="Maple Plot"><Plot height="400" type="two-dimensional" width="980" plot-scale="1.0" plot-xtrans="0.0" plot-ytrans="0.0" gridlinevisibility="1" legendvisibility="false">NiYtJSdDVVJWRVNHNiQ3XXM3JCQiIiEhIiIkIS9qJClSMioqWzwhIzk3JCQiLERjJ3lNOyEjNiQhLyczXS8iW1w5ISM5NyQkIitESmRwSyEjNSQhLy9CbkRRZTUhIzk3JCQiLUQxaydwMyUhIzckIS9OYm5HbW0kKSEjOjckJCIsdm9mViFcISM2JCEvRiI+VSlbQmchIzo3JCQiLXZvSHZAZCEjNyQhLyR6XyYzYSNmJCEjOjckJCIqRFkiUmwhIiokIS8uYTc5I0c2IiEjOjckJCIudm9IbDoneiEjOCQiLzdqJm96RD8kISM6NyQkIi12VlYpUlEqISM3JCIvQTx5eEpZdCEjOjckJCIvRDFSLmshMyIhIzgkIi9ELTh5Mzg2ISM5NyQkIi12VkEpR0EiISM2JCIvdGlkWC9SOSEjOTckJCIvdj1VIltHUSIhIzgkIi88NWBjI2ZzIiEjOTckJCIuRDEvOUdhIiEjNyQiLiNSJj5sWSI+ISM4NyQkIi9RJSkqKXB6QTshIzgkIi8vYGo+S3A+ISM5NyQkIi9EMVIqekZxIiEjOCQiL1wqUUc7bio+ISM5NyQkIi89bjg5eFU8ISM4JCIvKFxDZFMsKyMhIzk3JCQiLzdHKSlHdyN5IiEjOCQiLyIqcHZzdCcqPiEjOTckJCIvMSpHT2FGIz0hIzgkIi89TTI3YicpPiEjOTckJCIsdiRldWk9ISM1JCIvWnRTKlwncD4hIzk3JCQiLXYkNHNQLSMhIzYkIi8nM1A0KHpOPSEjOTckJCIqTil6JT0jISIpJCIvIiopPmFGVWciISM5NyQkIi1EMVkjZU0jISM2JCItZiU0JTQqRyIhIzc3JCQiLEQnMyZvXSMhIzUkIi8pNC9xI3AlMyohIzo3JCQiL0QxKnk2cm0jISM4JCIvLSEpKipbMWBbISM6NyQkIi5EY3JzdCNHISM3JCIvW3JRNXUwUyEjOzckJCIvdj1VT2ooKUghIzgkIS9sY3VYKXovJSEjOjckJCItdm9YKnk5JCEjNiQhLyQzdSNIaHYjKSEjOjckJCIsRCdwWidIJCEjNSQhL1R4VkNTIz0iISM5NyQkIi1EYyRmXVckISM2JCEvPWE4M0soWyIhIzk3JCQiKnZUT2YkISIpJCEvYmlkbSg0dCIhIzk3JCQiLXZWVEFVUCEjNiQhLyU0Rz5bWyE+ISM5NyQkIit2JFsiPlEhIiokIS9uXCplOWInPiEjOTckJCItRDFFMicqUSEjNiQhL0RKIWY+YisjISM5NyQkIix2JG8qSChSISM1JCEvbzheR2xDPyEjOTckJCItdm81IypcUyEjNiQhL0tPYlgjSC0jISM5NyQkIi1ESiZwUD8lISM2JCEvTndpOip6Jj4hIzk3JCQiLXYkKnpoZFYhIzYkIS8lPikqKT0mXCI9ISM5NyQkIi12JEhHbl4lISM2JCEuV2o5WURmIiEjODckJCItdiRmUWVuJSEjNiQhL0hfOnY3MDghIzk3JCQiLXYkKilbXCRbISM2JCEvWCcpcEssZScqISM6NyQkIi12JD5mUypcISM2JCEvVSE9RGFJKmUhIzo3JCQiK0QjZkU6JiEiKiQhLz4xKVJfaCM+ISM6NyQkIi1EYyNmN0omISM2JCIvPngkZlZqNCMhIzo3JCQiLHZHZilwYSEjNSQiLiNbeEJKP2chIzk3JCQiLXY9JGYlR2MhIzYkIi93TSlSODdxKiEjOjckJCIvREpYKilmInomISM4JCIvQFhMUXI0OCEjOTckJCIudj1kUVomZiEjNyQiL2ZlanB4KGYiISM5NyQkIi92Vik+eXk2JyEjOCQiL3ApUWcnb0M9ISM5NyQkIitEeSwiRychIiokIi8kbzNOKVwkKT4hIzk3JCQiL3YkNGNrR04nISM4JCIvUCcqUW1tST8hIzk3JCQiLnZvSDZaVSchIzckIi9GPVhBZ2o/ISM5NyQkIi9EIkcuZWxcJyEjOCQiL2whZU9FQTMjISM5NyQkIi12b1pTb2whIzYkIi9APUU8YyczIyEjOTckJCIuPm44R1ZnJyEjNyQiLz0hNEgwTTMjISM5NyQkIi92by86RFNtISM4JCIuQG9xSG4yIyEjODckJCIuY0UoWzx3bSEjNyQiLyNmWilcY20/ISM5NyQkIi5EMUMpNDduISM3JCIvendIbyVIMCMhIzk3JCQiL0Rjd1wlUnknISM4JCIvVk5lMV86PyEjOTckJCIsRHIiemJvISM1JCIvWSFlO1hbJz4hIzk3JCQiLkQiR2xiPHEhIzckIi90S0t2JFshPSEjOTckJCItdlY4S3pyISM2JCIvOWwlcGxwZSIhIzk3JCQiLnYkZmgzVHQhIzckIi9iJzNgJWU+OCEjOTckJCIrdjQmR10oISIqJCIvaztxVUU3NSEjOTckJCIsRCwhR2x3ISM1JCIvKCoqUlIzKFJuISM6NyQkIiowNHgjeSEiKSQiL3J4J2ZxJm9KISM6NyQkIix2M1EsKnohIzUkIS8iUXgib2MvWyEjOzckJCIrRHJjXyIpISIqJCEvQnM2diIqKjQlISM6NyQkIil3NDQkKSEiKCQhL0VEaGUvbHUhIzo3JCQiK3YhR2NZKSEiKiQhL1lZcyQ+RDEiISM5NyQkIipiZUBpKSEiKSQhL20rdDRrXTghIzk3JCQiK0QhKm95KCkhIiokIS9wKHlaVVhnIiEjOTckJCIudiRmJHlIMSohIzckIS8oKjNwND1pPiEjOTckJCItdiRwbnNNKiEjNiQhLzJRRyk0LDwjISM5NyQkIi5jRWR6PFYqISM3JCEvKFFjIltqLEEhIzk3JCQiL0RjXjlIOyYqISM4JCEvbnpybkA+QSEjOTckJCIucC9MLjNnKiEjNyQhL1RqXiRwSEEjISM5NyQkIi52JDRfSiZvKiEjNyQhLz0iPTwkNDhBISM5NyQkIi92PW4qUVYmKSohIzgkIS4iR0A4bWBAISM4NyQkIixERk9CKyIhIiokIS94J1xjbFMvIyEjOTckJCItREpMKDQuIiEjNSQhL0RONDAwYzwhIzk3JCQiKlI1J2Y1ISIoJCEvJFI7Nmk8TyIhIzk3JCQiL1FmVmlFdzUhIzckIS5AQG47XjQiISM4NyQkIi92PSg0QUg0IiEjNyQhLyVSWiQpPWEzKSEjOjckJCIvN3ldemQ0NiEjNyQhLik9SzNEIjMmISM5NyQkIi52ViFRQkU2ISM2JCEvSkdvQzApKj4hIzo3JCQiLzcuZHEoNDkiISM3JCIvTCwpejddXighIzs3JCQiL3ZvNC5zYjYhIzckIi8pKVtXdHAiWyQhIzo3JCQiL1FNaU5ZcTYhIzckIi95cmJlXGVoISM6NyQkIis6bz8mPSIhIikkIi8kNFI0QjV2KSEjOjckJCIvUSUpKlwjUSw3ISM3JCIvJVEjeikqZlk2ISM5NyQkIi92byU9ZXZAIiEjNyQiL2FYWlNgLDkhIzk3JCQiLzdgcFF0TDchIzckIi9LPkx0SlA7ISM5NyQkIi52VmI0Klw3ISM2JCIvdV9oIikqPSY9ISM5NyQkIi12M2RyIUciISM1JCIvSFxELSRvPiMhIzk3JCQiLkRKJz1fNjghIzYkIi87PS44MGFDISM5NyQkIi92VmBbbVY4ISM3JCIvKj0wKHkvSEUhIzk3JCQiLXZWeSFlUCIhIzUkIi88J3AiKXlTciMhIzk3JCQiL2Neb3I9JFEiISM3JCIvQCE9VUw9cyMhIzk3JCQiLzdHJFxtMFIiISM3JCIvN0skKkhaREYhIzk3JCQiL28vPWUlelIiISM3JCIvUSJlKj05REYhIzk3JCQiL0QiRzlEYFMiISM3JCItOlkqKik0cyMhIzc3JCQiL1FNI3okMz85ISM3JCIvT3NfXSU9cSMhIzk3JCQiLnY9V1VbViIhIzYkIi9iKXkpXEpwRSEjOTckJCIvRGNPM29tOSEjNyQiL2A8KHl1KGZEISM5NyQkIi1ESiM+JilcIiEjNSQiL1dxdydIdFMjISM5NyQkIi52JD46bWs6ISM2JCIuayUpZl52KyMhIzg3JCQiLkRjZFFBaSIhIzYkIi9cPTh0NEc7ISM5NyQkIi12dExVJW8iISM1JCIvRCQpKUdjJVI3ISM5NyQkIitiam1bPCEiKSQiL281SlJBUyopISM6NyQkIi12eWJeNj0hIzUkIi9ASHQoWzpDJyEjOjckJCIudlZWREIoPSEjNiQiLzdtI3A7M0UlISM6NyQkIi1ENlclKVI+ISM1JCIuUl1NdycqbyMhIzk3JCQiKzpLXis/ISIpJCIvPT9xJCk9RTwhIzo3JCQiLEQ2IUhsPyEiKiQiL0pURnAhUi8iISM6NyQkIi52JDR3KVI3IyEjNiQiL3pEJHoxbVgnISM7NyQkIix2WmYiKT0jISIqJCIvb2swRmtBUCEjOzckJCIudlY/UyZbQSEjNiQiLypIZTxMbzsjISM7NyQkIi52PVliO0ojISM2JCIvLExQSTAuNyEjOzckJCIsRDtpTFAjISIqJCIvMS5oKnBTaSchIzw3JCQiLnYkZkwnelYjISM2JCIvRVpvYD9xTSEjPDckJCIqKj49K0QhIigkIi80bSlcKWZDPSEjPDckJCItREUmNFFjIyEjNSQiLyJcKUhFU3EjKiEjPTckJCIudlY+NXBpIyEjNiQiL19iLDVoWlkhIz03JCQiK2JKKltvIyEiKSQiL3cjMydcI1xVIyEjPTckJCItRHIiWzh2IyEjNSQiL2ZucT9pSDYhIz03JCQiKkwneTVHISIoJCIvaDV0QFU4YyEjPjckJCIudlYhKWZUKEchIzYkIi9FajVdXz9FISM+NyQkIi5EY0k7WyRIISM2JCIvInp3eDNcQyIhIz43JCQiJCskISIiJCIvd0lrJ1FfXSYhIz8tJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lJVZJRVdHNiQ7JCIiISEiIiQiJCskISIiJShERUZBVUxURy0lK0FYRVNMQUJFTFNHNictSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1UScmIzk1Nzs2Ii8lJ2ZhbWlseUdRITYiLyUlc2l6ZUdRIzEwNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJmZhbHNlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1Enbm9ybWFsNiJRITYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM1JStIT1JJWk9OVEFMRyUrSE9SSVpPTlRBTEctJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiQhZiEiIi0lKUJPVU5EU19ZRzYjJCIjXSEiIi0lLUJPVU5EU19XSURUSEc2IyQiJSE0KiEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUhKVEhIiItJSlDSElMRFJFTkc2Ig==</Plot></Text-field>
</Output>
</Group>
<Group labelreference="L464" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># Five bound states exist in this case: the roots wrt nu are to be found:</Text-field>
</Input>
</Group>
<Group labelreference="L385" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">nu0:=fsolve(BesselJ(nu,g),nu=4..5);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEoJiM5NTc7MEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RIzo9RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y9LyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkwtRiM2KC1GLDYlUSdmc29sdmVGJ0YvRjItRjY2LVEwJkFwcGx5RnVuY3Rpb247RidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSYwLjBlbUYnL0ZORlgtSShtZmVuY2VkR0YkNiQtRiM2Ki1GIzYoLUYsNiVRKEJlc3NlbEpGJy9GMEY9RjlGVC1GZW42JC1GIzYoLUYsNiVRJyYjOTU3O0YnRl5vRjktRjY2LVEiLEYnRjlGOy9GP0YxRkBGQkZERkZGSEZXL0ZOUSwwLjMzMzMzMzNlbUYnLUkjbW5HRiQ2JFEwMTYuMDAwMDAwMDAwMDAwRidGOS8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGPUY5RjlGYHBGY3BGOUZmb0Zjb0Zmby1GIzYoLUZdcDYkUSI0RidGOS1GNjYtUSMuLkYnRjlGO0Y+RkBGQkZERkZGSC9GS1EsMC4yMjIyMjIyZW1GJ0ZZLUZdcDYkUSI1RidGOUZgcEZjcEY5RmBwRmNwRjlGOUZgcEZjcEY5RmBwRmNwRjk=">LUknZnNvbHZlRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiUtSShCZXNzZWxKR0YkNiRJI251R0YnJCIvKysrKysrOyEjN0YsOyIiJSIiJg==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L386" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">E0:=-nu0^2/(8.)/r0^2; # this is the energy for given nu.</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEjRTBGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSM6PUYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZMLUYjNictRjY2LVEqJnVtaW51czA7RidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSwwLjIyMjIyMjJlbUYnL0ZORlUtRiM2Ji1JI21uR0YkNiRRMzAuMDA3ODEyNTAwMDAwMDAwMEYnRjktRjY2LVExJkludmlzaWJsZVRpbWVzO0YnRjlGO0Y+RkBGQkZERkZGSC9GS1EmMC4wZW1GJy9GTkZbby1JJW1zdXBHRiQ2JS1GIzYoLUYsNiVRJ2Zzb2x2ZUYnRi9GMi1GNjYtUTAmQXBwbHlGdW5jdGlvbjtGJ0Y5RjtGPkZARkJGREZGRkhGam5GXG8tSShtZmVuY2VkR0YkNiQtRiM2Ki1GIzYoLUYsNiVRKEJlc3NlbEpGJy9GMEY9RjlGZW8tRmlvNiQtRiM2KC1GLDYlUScmIzk1NztGJ0ZicEY5LUY2Ni1RIixGJ0Y5RjsvRj9GMUZARkJGREZGRkhGam4vRk5RLDAuMzMzMzMzM2VtRictRlo2JFEwMTYuMDAwMDAwMDAwMDAwRidGOS8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGPUY5RjlGY3FGZnFGOUZqcEZncEZqcC1GIzYoLUZaNiRRIjRGJ0Y5LUY2Ni1RIy4uRidGOUY7Rj5GQEZCRkRGRkZIRlRGXG8tRlo2JFEiNUYnRjlGY3FGZnFGOUZjcUZmcUY5RjlGY3FGZnFGOS1GWjYkUSIyRidGOS8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRidGOUZjcUZmcUY5RmNxRmZxRjk=">LCQqJC1JJ2Zzb2x2ZUc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYlLUkoQmVzc2VsSkdGJjYkSSNudUdGKSQiLysrKysrKzshIzdGLjsiIiUiIiYiIiMkIS8rKysrXTd5ISM7</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L438" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">nu1:=fsolve(BesselJ(nu,g),nu=1.5..2);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEoJiM5NTc7MUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RIzo9RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y9LyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkwtRiM2KC1GLDYlUSdmc29sdmVGJ0YvRjItRjY2LVEwJkFwcGx5RnVuY3Rpb247RidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSYwLjBlbUYnL0ZORlgtSShtZmVuY2VkR0YkNiQtRiM2Ki1GIzYoLUYsNiVRKEJlc3NlbEpGJy9GMEY9RjlGVC1GZW42JC1GIzYoLUYsNiVRJyYjOTU3O0YnRl5vRjktRjY2LVEiLEYnRjlGOy9GP0YxRkBGQkZERkZGSEZXL0ZOUSwwLjMzMzMzMzNlbUYnLUkjbW5HRiQ2JFEwMTYuMDAwMDAwMDAwMDAwRidGOS8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGPUY5RjlGYHBGY3BGOUZmb0Zjb0Zmby1GIzYoLUZdcDYkUSQxLjVGJ0Y5LUY2Ni1RIy4uRidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSwwLjIyMjIyMjJlbUYnRlktRl1wNiRRIjJGJ0Y5RmBwRmNwRjlGYHBGY3BGOUY5RmBwRmNwRjlGYHBGY3BGOQ==">LUknZnNvbHZlRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiUtSShCZXNzZWxKR0YkNiRJI251R0YnJCIvKysrKysrOyEjN0YsOyQiIzohIiIiIiM=</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L437" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">E1:=-evalf(nu1^2/(8.)/r0^2);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEjRTFGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSM6PUYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZMLUYjNictRjY2LVEqJnVtaW51czA7RidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSwwLjIyMjIyMjJlbUYnL0ZORlUtRiM2Ji1JI21uR0YkNiRRMzAuMDA3ODEyNTAwMDAwMDAwMEYnRjktRjY2LVExJkludmlzaWJsZVRpbWVzO0YnRjlGO0Y+RkBGQkZERkZGSC9GS1EmMC4wZW1GJy9GTkZbby1JJW1zdXBHRiQ2JS1GIzYoLUYsNiVRJ2Zzb2x2ZUYnRi9GMi1GNjYtUTAmQXBwbHlGdW5jdGlvbjtGJ0Y5RjtGPkZARkJGREZGRkhGam5GXG8tSShtZmVuY2VkR0YkNiQtRiM2Ki1GIzYoLUYsNiVRKEJlc3NlbEpGJy9GMEY9RjlGZW8tRmlvNiQtRiM2KC1GLDYlUScmIzk1NztGJ0ZicEY5LUY2Ni1RIixGJ0Y5RjsvRj9GMUZARkJGREZGRkhGam4vRk5RLDAuMzMzMzMzM2VtRictRlo2JFEwMTYuMDAwMDAwMDAwMDAwRidGOS8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGPUY5RjlGY3FGZnFGOUZqcEZncEZqcC1GIzYoLUZaNiRRJDEuNUYnRjktRjY2LVEjLi5GJ0Y5RjtGPkZARkJGREZGRkhGVEZcby1GWjYkUSIyRidGOUZjcUZmcUY5RmNxRmZxRjlGOUZjcUZmcUY5RmByLyUxc3VwZXJzY3JpcHRzaGlmdEdRIjBGJ0Y5RmNxRmZxRjlGY3FGZnFGOQ==">LCQqJC1JJ2Zzb2x2ZUc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYlLUkoQmVzc2VsSkdGJjYkSSNudUdGKSQiLysrKysrKzshIzdGLjskIiM6ISIiIiIjRjYkIS8rKysrXTd5ISM7</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L462" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">nu2:=fsolve(BesselJ(nu,g),nu=5..6);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEoJiM5NTc7MkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RIzo9RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y9LyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkwtSSNtbkdGJDYkUTA1LjIyODQ5NzcxNDM1MTdGJ0Y5LyUrZm9yZWdyb3VuZEdRKFswLDAsMF1GJy8lKXJlYWRvbmx5R0Y9Rjk=">JCIvPE45eFxHXyEjOA==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L463" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">E2:=-evalf(nu2^2/(8.)/r0^2);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEjRTJGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSM6PUYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZMLUYjNiUtRjY2LVEqJnVtaW51czA7RidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSwwLjIyMjIyMjJlbUYnL0ZORlUtSSNtbkdGJDYkUTEwLjIxMzU3MTc4Mzk3NjQxRidGOUY5LyUrZm9yZWdyb3VuZEdRKFswLDAsMF1GJy8lKXJlYWRvbmx5R0Y9Rjk=">JCEvVHdSeXJOQCEjOQ==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L466" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">nu3:=fsolve(BesselJ(nu,g),nu=2..4);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEoJiM5NTc7M0YnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RIzo9RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y9LyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkwtSSNtbkdGJDYkUTAyLjg0MTY3OTI1OTE1MjRGJ0Y5LyUrZm9yZWdyb3VuZEdRKFswLDAsMF1GJy8lKXJlYWRvbmx5R0Y9Rjk=">JCIvQzpmI3o7JUchIzg=</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L465" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">E3:=-evalf(nu3^2/(8.)/r0^2);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEjRTNGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSM6PUYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZMLUYjNiUtRjY2LVEqJnVtaW51czA7RidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSwwLjIyMjIyMjJlbUYnL0ZORlUtSSNtbkdGJDYkUTIwLjA2MzA4NzAzOTE1NTQ0NUYnRjlGOS8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGPUY5">JCEvWGE6UnEzaiEjOg==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L468" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">nu4:=fsolve(BesselJ(nu,g),nu=0.5..1.5);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEoJiM5NTc7NEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RIzo9RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y9LyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkwtSSNtbkdGJDYkUTEwLjY5MDQxNzU1NjUxMDIwRidGOS8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGPUY5">JCIvPzVsYjwvcCEjOQ==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L467" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">E4:=-evalf(nu4^2/(8.)/r0^2);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEjRTRGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSM6PUYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZMLUYjNiUtRjY2LVEqJnVtaW51czA7RidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSwwLjIyMjIyMjJlbUYnL0ZORlUtSSNtbkdGJDYkUTMwLjAwMzcyNDAzNDM5MzI2MTlGJ0Y5RjkvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRj1GOQ==">JCEvPkUkUk1TcyQhIzs=</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L459" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"></Text-field>
</Input>
</Group>
<Group labelreference="L461" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># We can study systematically for each of the 5 eigenvalues how well the numerical calculation is doing.</Text-field>
</Input>
</Group>
<Group labelreference="L460" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"></Text-field>
</Input>
</Group>
<Group labelreference="L390" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># We optimize the calculation wrt the basis parameter.</Text-field>
</Input>
</Group>
<Group labelreference="L391" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">phi:=n->uL(n,0,1.0);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEnJiM5NjY7RicvJSdpdGFsaWNHUSZmYWxzZUYnLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGJy1JI21vR0YkNi1RIzo9RidGMi8lJmZlbmNlR0YxLyUqc2VwYXJhdG9yR0YxLyUpc3RyZXRjaHlHRjEvJSpzeW1tZXRyaWNHRjEvJShsYXJnZW9wR0YxLyUubW92YWJsZWxpbWl0c0dGMS8lJ2FjY2VudEdGMS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkktRiM2KC1GLDYlUSJuRicvRjBRJXRydWVGJy9GM1EnaXRhbGljRictRjY2LVEoJiM4NTk0O0YnRjJGOUY7Rj1GP0ZBRkNGRS9GSFEmMC4wZW1GJy9GS0ZZLUYjNigtRiw2JVEjdUxGJ0ZRRlMtRjY2LVEwJkFwcGx5RnVuY3Rpb247RidGMkY5RjtGPUY/RkFGQ0ZFRlhGWi1JKG1mZW5jZWRHRiQ2JC1GIzYqRk4tRjY2LVEiLEYnRjJGOS9GPEZSRj1GP0ZBRkNGRUZYL0ZLUSwwLjMzMzMzMzNlbUYnLUkjbW5HRiQ2JFEiMEYnRjJGYm8tRmlvNiRRJDEuMEYnRjIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRjFGMkYyRl9wRmJwRjJGX3BGYnBGMkZfcEZicEYy">Zio2I0kibkc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkjdUxHRiU2JTkkIiIhJCIjNSEiIkYlRiVGJQ==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L404" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag:=proc(Nb,beta,N_E) local phi,VME,TME,OL,VM,HM,OM,m,n,res; phi:=n->uL(n,0,beta);</Text-field>
</Input>
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">VME:=(m,n)->int(phi(n)*V*phi(m),r=0..infinity); TME:=(m,n)->-1/2*int(phi(n)*diff(phi(m),r$2),r=0..infinity); OL:=(m,n)->int(phi(n)*phi(m),r=0..infinity);</Text-field>
</Input>
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">VM:=Matrix(Nb,Nb,shape=symmetric):</Text-field>
</Input>
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">HM:=Matrix(Nb,Nb,shape=symmetric):</Text-field>
</Input>
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">OM:=Matrix(Nb,Nb,shape=symmetric):</Text-field>
<Text-field prompt="> " style="Maple Input" layout="Normal">for m from 1 to Nb do: for n from 1 to m do: VM[m,n]:=VME(m,n): OM[m,n]:=OL(m,n): od: od:</Text-field>
<Text-field prompt="> " style="Maple Input" layout="Normal">for m from 1 to Nb do: for n from 1 to m do: HM[m,n]:=evalf(VM[m,n])+evalf(TME(m,n)): od: od:</Text-field>
</Input>
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">res:=(convert(sort(Re(Eigenvalues(HM,OM))),list)); [seq(res[j],j=1..N_E)];</Text-field>
</Input>
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">end:</Text-field>
</Input>
</Group>
<Group labelreference="L402" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"></Text-field>
</Input>
</Group>
<Group labelreference="L395" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag(5,1,5);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">NyckITMhKipHKFItP08vNSEjPCQhMz5nd0JxYlpMWiEjPSQhMyYqcC4kZiswIVw8RigkITNHIVtGJD5ONmhQISM+JCIzOCFSUyhRMHdeXEYo</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L396" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag(5,1/2,5);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">NyckITNxUlZZIj1WaCEqKiEjPSQhMyIpSE82LiZlWkUlRiUkITMxU3M7ckopPkEiRiUkITMpKnB1YzMoUkZUIiEjPiQiMyUpcCJ5QCYpPSEzRyEjPw==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L469" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># A smaller beta-value results in worse (higher) energies for all states, except the last one (which is still coming out as unbound). Note this 5-by-5 calculation has only 5 eigenvalues!</Text-field>
</Input>
</Group>
<Group labelreference="L470" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">j:=4; # pick the state for which the analysis is carried out:</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output" display="LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEiakYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RIzo9RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y9LyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkwtSSNtbkdGJDYkUSI0RidGOS8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGPUY5">IiIl</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L397" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">LL5:=[]: for i from 1 to 20 do: beta:=3/10*i; LL5:=[op(LL5),[beta,Diag(5,beta,5)[j]]]: od:</Text-field>
</Input>
</Group>
<Group labelreference="L405" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">PL1:=plot(LL5): beta:='beta': PL0:=plot(E||(j-1),beta=0..2,color=green): plots[display](PL0,PL1,view=[0..2,-0.1..0.2]);</Text-field>
</Input>
<Output>
<Text-field style="Maple Plot" layout="Maple Plot"><Plot height="400" type="two-dimensional" width="756" plot-scale="1.0" plot-xtrans="0.0" plot-ytrans="0.0" gridlinevisibility="1" legendvisibility="false">NictJSdDVVJWRVNHNiQ3UzckJCIiISEiIiQhL1hhOlJxM2ohIzo3JCQiMUxMTEwzVmZWISM8JCEvWGE6UnEzaiEjOjckJCIxbm1tIkhbRDopISM8JCEvWGE6UnEzaiEjOjckJCIyTExMJGUwJD1DIiEjPCQhL1hhOlJxM2ohIzo3JCQiMkxMTCQzUkJyOyEjPCQhL1hhOlJxM2ohIzo3JCQiMmxtbSJ6amYpNCMhIzwkIS9YYTpScTNqISM6NyQkIjJLTExlNDtbXCMhIzwkIS9YYTpScTNqISM6NyQkIjIoKioqKlxpJ3ldIUghIzwkIS9YYTpScTNqISM6NyQkIjFMTCRlenMkSEwhIzskIS9YYTpScTNqISM6NyQkIjImKioqKlw3aUlfUCEjPCQhL1hhOlJxM2ohIzo3JCQiMW5tbTtfTSg9JSEjOyQhL1hhOlJxM2ohIzo3JCQiMk1MTCQzeV9xWCEjPCQhL1hhOlJxM2ohIzo3JCQiKmwrPismISIqJCEvWGE6UnEzaiEjOjckJCIqdlddViYhIiokIS9YYTpScTNqISM6NyQkIipOZkMmZSEiKiQhL1hhOlJxM2ohIzo3JCQiMUxMJGV6NjpCJyEjOyQhL1hhOlJxM2ohIzo3JCQiMW1tbTs9QyNvJyEjOyQhL1hhOlJxM2ohIzo3JCQiMW1tbW0jcFMxKCEjOyQhL1hhOlJxM2ohIzo3JCQiLERPRCMzdiEjNiQhL1hhOlJxM2ohIzo3JCQiMW1tbW0oeTgheiEjOyQhL1hhOlJxM2ohIzo3JCQiLERPSUZMKSEjNiQhL1hhOlJxM2ohIzo3JCQiLHYzek11KSEjNiQhL1hhOlJxM2ohIzo3JCQiMW5tbSJIXz88KiEjOyQhL1hhOlJxM2ohIzo3JCQiMW5tO3ppaGwmKiEjOyQhL1hhOlJxM2ohIzo3JCQiMUxMTDMjRywqKiohIzskIS9YYTpScTNqISM6NyQkIjJLTCRlenc1VjUhIzskIS9YYTpScTNqISM6NyQkIi12JFEjXCIzIiEjNiQhL1hhOlJxM2ohIzo3JCQiMktMTGUiKltINyIhIzskIS9YYTpScTNqISM6NyQkIipkeGQ7IiEiKSQhL1hhOlJxM2ohIzo3JCQiMikqKioqXF9xbjI3ISM7JCEvWGE6UnEzaiEjOjckJCIyKSoqKlxpJnBAWzchIzskIS9YYTpScTNqISM6NyQkIjIpKioqKlwyJ0hLSCIhIzskIS9YYTpScTNqISM6NyQkIjJsbW1tWnZPTCIhIzskIS9YYTpScTNqISM6NyQkIit2KydvUCIhIiokIS9YYTpScTNqISM6NyQkIjJLTCRlUjwqZlQiISM7JCEvWGE6UnEzaiEjOjckJCIrJilIeGU5ISIqJCEvWGE6UnEzaiEjOjckJCIybG07SCFvLSpcIiEjOyQhL1hhOlJxM2ohIzo3JCQiMioqKipcN2suNmEiISM7JCEvWGE6UnEzaiEjOjckJCIybW1tO1dUQWUiISM7JCEvWGE6UnEzaiEjOjckJCItRDEqM2BpIiEjNiQhL1hhOlJxM2ohIzo3JCQiMk1MTEwqenltOyEjOyQhL1hhOlJxM2ohIzo3JCQiMkxMTDNOMSM0PCEjOyQhL1hhOlJxM2ohIzo3JCQiMm1tO0hZdDd2IiEjOyQhL1hhOlJxM2ohIzo3JCQiKnhHKip5IiEiKSQhL1hhOlJxM2ohIzo3JCQiMm1tbVQ2S1UkPSEjOyQhL1hhOlJxM2ohIzo3JCQiMkxMTExiZFEoPSEjOyQhL1hhOlJxM2ohIzo3JCQiLURPbDU7PiEjNiQhL1hhOlJxM2ohIzo3JCQiLXYuVWFjPiEjNiQhL1hhOlJxM2ohIzo3JCQiIz8hIiIkIS9YYTpScTNqISM6LSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIjNSEiIiQiIiEhIiItJSdDVVJWRVNHNiQ3NjckJCIiJCEiIiQiMEpzW0dFdV0nISM9NyQkIiInISIiJCEwdkJBR1knKlIjISM7NyQkIiIqISIiJCEwUCxvMiNHSkghIzs3JCQiIzchIiIkIS8memFZcyRIYSEjOjckJCIjOiEiIiQhMEs0WUBpcSVSISM7NyQkIiM9ISIiJCIwLDdjUjdQPiIhIzo3JCQiI0AhIiIkIjAtNys5KmVvVSEjOjckJCIjQyEiIiQiMCY0aFUmKnpJJSkhIzo3JCQiI0YhIiIkIjBMMyVbJz4mXDghIzk3JCQiI0khIiIkIjBUc0lsU3ckPiEjOTckJCIjTCEiIiQiMFZHPW85SGcjISM5NyQkIiNPISIiJCIwTDdGTG5FTSQhIzk3JCQiI1IhIiIkIjAlZShScSo9YlQhIzk3JCQiI1UhIiIkIi8lKip5YkIkUl0hIzg3JCQiI1ghIiIkIjB1VjxgWlUqZiEjOTckJCIjWyEiIiQiMFJSUi5lJD5xISM5NyQkIiNeISIiJCIvQSdRaCk+OSIpISM4NyQkIiNhISIiJCIweWoxKmVUeSMqISM5NyQkIiNkISIiJCIwWlNzM3Q2MCIhIzg3JCQiI2chIiIkIjB5M0AnPlIiPSIhIzgtJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lJVZJRVdHNiQ7JCIiISEiIiQiIz8hIiI7JCEiIiEiIiQiIiMhIiItJStBWEVTTEFCRUxTRzYnLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEnJiM5NDY7NiIvJSdmYW1pbHlHUShERUZBVUxUNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EmZmFsc2U2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdub3JtYWw2IlEhNiItJSVGT05URzYlJShERUZBVUxURyUoREVGQVVMVEciIzUlK0hPUklaT05UQUxHJStIT1JJWk9OVEFMRy0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJHEmISIiLSUpQk9VTkRTX1lHNiMkIiQ/IiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJTVwISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJV1QISIiLSUpQ0hJTERSRU5HNiI=</Plot></Text-field>
</Output>
</Group>
<Group labelreference="L406" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">j;# </Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">IiIl</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L407" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">#LL5;</Text-field>
</Input>
</Group>
<Group labelreference="L408" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># Exercise: find the minimum value from the 5-by-5 calculation by optimizing wrt beta; how big is the relative error?</Text-field>
</Input>
</Group>
<Group labelreference="L471" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># Try 10-by-10 calculations:</Text-field>
</Input>
</Group>
<Group labelreference="L414" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">LL10:=[]: for i from 1 to 5 do: beta:=1+1/10*i; LL10:=[op(LL10),[beta,Diag(10,beta,5)[j]]]: od:</Text-field>
</Input>
</Group>
<Group labelreference="L412" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">PL2:=plot(LL10,color=blue,style=point,symbol=cross): plots[display](PL0,PL1,PL2,view=[0..2,-0.065..-0.035]);</Text-field>
</Input>
<Output>
<Text-field style="Maple Plot" layout="Maple Plot"><Plot height="400" type="two-dimensional" width="944" plot-scale="1.0" plot-xtrans="0.0" plot-ytrans="0.0" gridlinevisibility="1" legendvisibility="false">NigtJSdDVVJWRVNHNiQ3UzckJCIiISEiIiQhL1hhOlJxM2ohIzo3JCQiMUxMTEwzVmZWISM8JCEvWGE6UnEzaiEjOjckJCIxbm1tIkhbRDopISM8JCEvWGE6UnEzaiEjOjckJCIyTExMJGUwJD1DIiEjPCQhL1hhOlJxM2ohIzo3JCQiMkxMTCQzUkJyOyEjPCQhL1hhOlJxM2ohIzo3JCQiMmxtbSJ6amYpNCMhIzwkIS9YYTpScTNqISM6NyQkIjJLTExlNDtbXCMhIzwkIS9YYTpScTNqISM6NyQkIjIoKioqKlxpJ3ldIUghIzwkIS9YYTpScTNqISM6NyQkIjFMTCRlenMkSEwhIzskIS9YYTpScTNqISM6NyQkIjImKioqKlw3aUlfUCEjPCQhL1hhOlJxM2ohIzo3JCQiMW5tbTtfTSg9JSEjOyQhL1hhOlJxM2ohIzo3JCQiMk1MTCQzeV9xWCEjPCQhL1hhOlJxM2ohIzo3JCQiKmwrPismISIqJCEvWGE6UnEzaiEjOjckJCIqdlddViYhIiokIS9YYTpScTNqISM6NyQkIipOZkMmZSEiKiQhL1hhOlJxM2ohIzo3JCQiMUxMJGV6NjpCJyEjOyQhL1hhOlJxM2ohIzo3JCQiMW1tbTs9QyNvJyEjOyQhL1hhOlJxM2ohIzo3JCQiMW1tbW0jcFMxKCEjOyQhL1hhOlJxM2ohIzo3JCQiLERPRCMzdiEjNiQhL1hhOlJxM2ohIzo3JCQiMW1tbW0oeTgheiEjOyQhL1hhOlJxM2ohIzo3JCQiLERPSUZMKSEjNiQhL1hhOlJxM2ohIzo3JCQiLHYzek11KSEjNiQhL1hhOlJxM2ohIzo3JCQiMW5tbSJIXz88KiEjOyQhL1hhOlJxM2ohIzo3JCQiMW5tO3ppaGwmKiEjOyQhL1hhOlJxM2ohIzo3JCQiMUxMTDMjRywqKiohIzskIS9YYTpScTNqISM6NyQkIjJLTCRlenc1VjUhIzskIS9YYTpScTNqISM6NyQkIi12JFEjXCIzIiEjNiQhL1hhOlJxM2ohIzo3JCQiMktMTGUiKltINyIhIzskIS9YYTpScTNqISM6NyQkIipkeGQ7IiEiKSQhL1hhOlJxM2ohIzo3JCQiMikqKioqXF9xbjI3ISM7JCEvWGE6UnEzaiEjOjckJCIyKSoqKlxpJnBAWzchIzskIS9YYTpScTNqISM6NyQkIjIpKioqKlwyJ0hLSCIhIzskIS9YYTpScTNqISM6NyQkIjJsbW1tWnZPTCIhIzskIS9YYTpScTNqISM6NyQkIit2KydvUCIhIiokIS9YYTpScTNqISM6NyQkIjJLTCRlUjwqZlQiISM7JCEvWGE6UnEzaiEjOjckJCIrJilIeGU5ISIqJCEvWGE6UnEzaiEjOjckJCIybG07SCFvLSpcIiEjOyQhL1hhOlJxM2ohIzo3JCQiMioqKipcN2suNmEiISM7JCEvWGE6UnEzaiEjOjckJCIybW1tO1dUQWUiISM7JCEvWGE6UnEzaiEjOjckJCItRDEqM2BpIiEjNiQhL1hhOlJxM2ohIzo3JCQiMk1MTEwqenltOyEjOyQhL1hhOlJxM2ohIzo3JCQiMkxMTDNOMSM0PCEjOyQhL1hhOlJxM2ohIzo3JCQiMm1tO0hZdDd2IiEjOyQhL1hhOlJxM2ohIzo3JCQiKnhHKip5IiEiKSQhL1hhOlJxM2ohIzo3JCQiMm1tbVQ2S1UkPSEjOyQhL1hhOlJxM2ohIzo3JCQiMkxMTExiZFEoPSEjOyQhL1hhOlJxM2ohIzo3JCQiLURPbDU7PiEjNiQhL1hhOlJxM2ohIzo3JCQiLXYuVWFjPiEjNiQhL1hhOlJxM2ohIzo3JCQiIz8hIiIkIS9YYTpScTNqISM6LSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIjNSEiIiQiIiEhIiItJSdDVVJWRVNHNiQ3NjckJCIiJCEiIiQiMEpzW0dFdV0nISM9NyQkIiInISIiJCEwdkJBR1knKlIjISM7NyQkIiIqISIiJCEwUCxvMiNHSkghIzs3JCQiIzchIiIkIS8memFZcyRIYSEjOjckJCIjOiEiIiQhMEs0WUBpcSVSISM7NyQkIiM9ISIiJCIwLDdjUjdQPiIhIzo3JCQiI0AhIiIkIjAtNys5KmVvVSEjOjckJCIjQyEiIiQiMCY0aFUmKnpJJSkhIzo3JCQiI0YhIiIkIjBMMyVbJz4mXDghIzk3JCQiI0khIiIkIjBUc0lsU3ckPiEjOTckJCIjTCEiIiQiMFZHPW85SGcjISM5NyQkIiNPISIiJCIwTDdGTG5FTSQhIzk3JCQiI1IhIiIkIjAlZShScSo9YlQhIzk3JCQiI1UhIiIkIi8lKip5YkIkUl0hIzg3JCQiI1ghIiIkIjB1VjxgWlUqZiEjOTckJCIjWyEiIiQiMFJSUi5lJD5xISM5NyQkIiNeISIiJCIvQSdRaCk+OSIpISM4NyQkIiNhISIiJCIweWoxKmVUeSMqISM5NyQkIiNkISIiJCIwWlNzM3Q2MCIhIzg3JCQiI2chIiIkIjB5M0AnPlIiPSIhIzgtJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lJ0NVUlZFU0c2JjcnNyQkIiM2ISIiJCEvdy1PJXklemkhIzo3JCQiIzchIiIkITBMel5lYkFJJyEjOzckJCIjOCEiIiQhMEpaMCFRbixqISM7NyQkIiM5ISIiJCEwS1ZaZCkpPkknISM7NyQkIiM6ISIiJCEwInAmbzQqKltJJyEjOy0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiM1ISIiLSUnU1lNQk9MRzYjJSZDUk9TU0ctJSZTVFlMRUc2IyUmUE9JTlRHLSUlVklFV0c2JDskIiIhISIiJCIjPyEiIjskISNsISIkJCEjTiEiJC0lK0FYRVNMQUJFTFNHNictSSZtZnJhY0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjYqLUkjbW5HNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEiMzYiLyUnZmFtaWx5R1EoREVGQVVMVDYiLyUlc2l6ZUdRIzEwNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJmZhbHNlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1Enbm9ybWFsNiItSSNtbkc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USIyNiIvJSdmYW1pbHlHUShERUZBVUxUNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EmZmFsc2U2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdub3JtYWw2Ii8lLmxpbmV0aGlja25lc3NHUSIxNiIvJStkZW5vbWFsaWduR1EnY2VudGVyNiIvJSludW1hbGlnbkdRJ2NlbnRlcjYiLyUpYmV2ZWxsZWRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiJRITYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM1JStIT1JJWk9OVEFMRyUrSE9SSVpPTlRBTEctJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiRJKCEiIi0lKUJPVU5EU19ZRzYjJCIkPychIiItJS1CT1VORFNfV0lEVEhHNiMkIiVJJykhIiItJS5CT1VORFNfSEVJR0hURzYjJCIlXUshIiItJSlDSElMRFJFTkc2Ig==</Plot></Text-field>
</Output>
</Group>
<Group labelreference="L413" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"></Text-field>
</Input>
</Group>
<Group labelreference="L409" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># we see that the 10-by-10 diagonalization is less sensitive to the basis parameter beta.</Text-field>
</Input>
</Group>
<Group labelreference="L410" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># Exercise: track the relative error (compared to the exact answer) as a function of beta.</Text-field>
</Input>
</Group>
<Group labelreference="L427" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"></Text-field>
</Input>
</Group>
<Group labelreference="L411" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># Now ask the questions: is there a bound state in the L=1 sector (p-wave)?</Text-field>
</Input>
</Group>
<Group labelreference="L415" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"></Text-field>
</Input>
</Group>
<Group labelreference="L422" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag:=proc(Nb,beta,L,N_E) local phi,VME,TME,OL,VM,HM,OM,m,n,res; phi:=n->uL(n,L,beta);</Text-field>
</Input>
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">VME:=(m,n)->int(phi(n)*(V+L*(L+1)/2/r^2)*phi(m),r=0..infinity); TME:=(m,n)->-1/2*int(phi(n)*diff(phi(m),r$2),r=0..infinity); OL:=(m,n)->int(phi(n)*phi(m),r=0..infinity);</Text-field>
</Input>
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">VM:=Matrix(Nb,Nb,shape=symmetric):</Text-field>
</Input>
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">HM:=Matrix(Nb,Nb,shape=symmetric):</Text-field>
</Input>
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">OM:=Matrix(Nb,Nb,shape=symmetric):</Text-field>
<Text-field prompt="> " style="Maple Input" layout="Normal">for m from 1 to Nb do: for n from 1 to m do: VM[m,n]:=VME(m,n): OM[m,n]:=OL(m,n): od: od:</Text-field>
<Text-field prompt="> " style="Maple Input" layout="Normal">for m from 1 to Nb do: for n from 1 to m do: HM[m,n]:=evalf(VM[m,n])+evalf(TME(m,n)): od: od:</Text-field>
</Input>
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">res:=(convert(sort(Re(Eigenvalues(HM,OM))),list)); [seq(res[j],j=1..N_E)];</Text-field>
</Input>
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">end:</Text-field>
</Input>
</Group>
<Group labelreference="L416" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"></Text-field>
</Input>
</Group>
<Group labelreference="L417" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag(5,1,1,1);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">NyMkITNsUkhZKHB6ISlmJyEjPQ==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L418" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag(5,2,1,1);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">JCEzKikqXGlIVlhObCchIz0=</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L419" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag(5,3,1,1);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">JCEzJSpSb104LS9hbSEjPQ==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L420" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag(5,4,1,3);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">NyUkITNUU0FwaUU0Vm0hIz0kITM1NXg6YTlkVj5GJSQiM19IaSp6OyxALSlGJQ==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L452" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag(10,4,1,3);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">NyUkITNNZ0w8RUYvYW0hIz0kITNGU1ApKXpqNSY0JEYlJCEzTj5meTFLWiZlKSEjPg==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L453" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag(20,4,1,3);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">NyUkITNwPnZgdUYvYW0hIz0kITN5PmY3LF8zKDQkRiUkITMrK3RaQjM7NTZGJQ==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L454" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag(30,4,1,4);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">NyYkITNwUkFkdUYvYW0hIz0kITMkKkgibzg/JjMoNCRGJSQhMzAhKTRGPTlANTZGJSQhMy4hWyVvUDhUIXAiISM+</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L455" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag(30,5,1,5);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">NyckITNvejg1dUYvYW0hIz0kITNDPz8+LF8zKDQkRiUkITMwUyNlJmUlNC02IkYlJCEzKikqZTMqNCdIJiplIiEjPiQiM0tdWi1fNTJCTkYs</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L456" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag(30,3,1,5);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">NyckITNJSTRrdUYvYW0hIz0kITMyPyNlNz8mMyg0JEYlJCEzJSpmXCM9VTYtNiJGJSQhMykqKkg9ckVtLXEiISM+JCIzRCFlIW88XXJRbiEjPw==</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L457" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal">Diag(30,2,1,5);</Text-field>
</Input>
<Output>
<Text-field style="2D Output" layout="Maple Output"><Equation executable="false" style="2D Output">NyckITNCNU5ddUYvYW0hIz0kITMjKSoqPkosXzMoNCRGJSQhMyw/dyU9VTYtNiJGJSQhMzs1cXB6Q1crPCEjPiQiMy0/aUAoSCpcN0QhIz8=</Equation></Text-field>
</Output>
</Group>
<Group labelreference="L458" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"></Text-field>
</Input>
</Group>
<Group labelreference="L472" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"># For L>0 we do not have an exact result to compare against. Thus, we need to be diligent when determining the best results from matrix diagonalization.</Text-field>
</Input>
</Group>
<Group labelreference="L473" drawlabel="true">
<Input>
<Text-field prompt="> " style="Maple Input" layout="Normal"></Text-field>
</Input>
</Group>
</Worksheet>