Half Linear Potential - The Answers

Read in the zeros and normalization constants

data={{2.33811,1.4261},{4.08795,1.24516},{5.52056,1.1558},{6.78671,1.09787},{7.94413,1.05559},{9.02265,1.02258},{10.04017,0.995649},{11.00853,0.97301},{11.93601,0.953543},{12.82877,0.93651},{13.69149,0.921402},{14.52783,0.907849},{15.34075,0.895579},{16.13268,0.884383},{16.90563,0.874098},{17.6613,0.864596},{18.40113,0.855773},{19.12638,0.847543},{19.83813,0.839838},{20.53734,0.832597}};

Define (Subscript[ψ, n])(x) .

ψ[n_,x_]:=data[[n,2]]*AiryAi[x-data[[n,1]]];

Graph (Subscript[Ψ, n])(x,0).

Manipulate[
Show[
Plot[
ψ[n,x],{x,0,25},
PlotRange->{{0,25},{-0.55,0.8}},
Ticks->{  {{25,""}},{ {-0.4,""},{-0.2,""},{0.2,""},{0.4,""},{0.6,""},{0.8,""}}},
PlotStyle->{Orange,Thickness[0.005]},
AxesLabel->{x,"ψ(x,0)"},
ImageSize->800,
LabelStyle->({FontSize->16,FontFamily->"Times"})
],
ListLinePlot[
{{data[[n,1]],-0.55},{data[[n,1]],+0.8}},
PlotStyle->{Blue,Thickness[0.001],Dotted}
]
],
{{n,1,"n"},Range[1,20,1],Setter}
]

Graph ((||Ψ(x)||)^2).

Manipulate[
Show[
Plot[
Conjugate[ψ[n,x]]*ψ[n,x],{x,0,25},
PlotRange->{{0,25},{-0.1,0.6}},
Ticks->{  {{25,""}},{ {0.1,""},{0.2,""},{0.3,""},{0.4,""},{0.5,""},{0.6,""}}},
PlotStyle->{Orange,Thickness[0.005]},
AxesLabel->{x,"\[LeftBracketingBar]Ψ(\[RightBracketingBar]^2)"},
ImageSize->800,
LabelStyle->({FontSize->16,FontFamily->"Times"})
],
ListLinePlot[
{{data[[n,1]],0},{data[[n,1]],+0.6}},
PlotStyle->{Blue,Thickness[0.001],Dotted}
]
],
{{n,1,"n"},Range[1,20,1],Setter}
]