(*Задание Хюльтеновского анзаца, энергий связи и масс частиц*) 
m1=938.2723;
m2=939.5656;
eeD=2.225;
(*Eb2=1.296;
Eb3=7.77;*)
kProton=(Sqrt[2*eeD*((m2*m1)/(m2+m1))]);
j=197.327;
radCh = 2.1421;

Solve[(radii^2)/4 + 0.64 -0.125==(radCh^2),radii,Assumptions->{radii>0}]

Psi[r_, q_, k_]:=((E^(-k r))-(E^(-q r)))

(*Нормировка*)

Norm1[q_,k_] := Integrate[((Psi[r,q,k])^2),{r,0,Infinity},Assumptions->{k>0,q>0}]

PsiNorm[r_,q_,k_] := ((Psi[r,q,k])/(Sqrt[Norm1[q,k]]))

Integrate[((PsiNorm[r,q,k])^2),{r,0,Infinity},Assumptions->{k>0,q>0}]

(*Задание численного значения параметра q для Eb=2.224 МэВ. Поиск точки, в окрестности которой будем искать решение уравнения относительно q*)
Plot[Integrate[((PsiNorm[r,(q/j),(kProton/j)])^2)*(r^2),{r,0,Infinity}]-(4.04^2),{q,225,230}]

FindRoot[Integrate[((PsiNorm[r,(q/j),(kProton/j)])^2)*(r^2),{r,0,Infinity},Assumptions->Re[q]>0]==(4.04^2),{q,228,229}]

PsiNorm[r,((228.4649720223125`)/j),(kProton/j)]

PsiD[r_]:=PsiNorm[r,((228.4649720223125`)/j),(kProton/j)]

Plot[PsiD[rad],{rad,0,20}]

Psi2P[p_]:=Simplify[((Sqrt[4 π])/p)Integrate[PsiD[r] Sin[p r],{r,0,Infinity}],Assumptions->p>0]

Print[Psi2P[p]]

(*Проверка нормировки получившейся ВФ*)
Simplify[Integrate[4 π (p^2) ((Psi2P[p])^2) (1/(((2 π))^3)),{p,0,Infinity}]]