# How to do Interpolation to this function after given b this values

Mathematica Asked by user69941 on August 29, 2020

  [Chi][b_] := (Z/k)*NIntegrate[Subscript[f, pp1][q]*Subscript[F, pp1][q]*BesselJ[0, q*b]*q, {q, 0, b}];


bList = Table[b, {b, 0, 5.1 , 0.05}]

where:
Z = 6;
k=1.55246;

 Subscript[f, pp1][q_]=(6.254736279890945*(1 - 1.4511668475476842))/E^(0.2115*q^2)

Subscript[F, pp1][q_] := ((4*Pi)/q)*NIntegrate[Subscript[[Rho], p][r]*Sin[q*r]*r, {r, 0, Infinity}];

Subscript[[Rho], p][r_]=0.013132593248303927/(1 + (E^(1.7543859649122808*(-2.380427976610103 - r)) + E^(1.7543859649122808*(-2.380427976610103 + r)))*
(0.5 + 0.08823886490842314*r^2)^1.5)


Clear["Global*"]

Subscript[f, pp1][q_] = (6.254736279890945*(1 - 1.4511668475476842))/
E^(0.2115*q^2);

Subscript[F, pp1][q_?NumericQ] := ((4*Pi)/q)*
NIntegrate[Subscript[ρ, p][r]*Sin[q*r]*r, {r, 0, Infinity}];

Subscript[ρ, p][r_] =
0.013132593248303927/(1 +
(E^(1.7543859649122808*(-2.380427976610103 - r)) +
E^(1.7543859649122808*(-2.380427976610103 + r)))*(0.5 +
0.08823886490842314*r^2)^1.5);

χ[b_?NumericQ] := (Z/k)*NIntegrate[
Subscript[f, pp1][q]*Subscript[F, pp1][q]*BesselJ[0, q*b]*q, {q, 0, b}];

Z = 6; k = 1.55246;

data = {#, χ[#]} & /@ Range[0, 5.1, 0.05] // Quiet;

f = Interpolation[data];

Plot[f[b], {b, 0, 5.1},
Epilog -> {Red, AbsolutePointSize, Point[data]}]
` Answered by Bob Hanlon on August 29, 2020

## Related Questions

### Fit, NIntegrate in FindRoot results in non-numerical values

0  Asked on August 2, 2021

### How do you calculate the area of a surface with edges in mathematica

1  Asked on August 2, 2021

### Generate 3 dimensional random vectors

2  Asked on August 2, 2021

### Reduce : This system cannot be solved with the methods available to Reduce

0  Asked on August 2, 2021

### Sorting Eigensystem According to Complicated Rule

1  Asked on August 2, 2021

### How to find the regions for which the derivative of my function changes sign?

0  Asked on August 1, 2021

### Map Function with Table on each row

1  Asked on August 1, 2021 by smarter-math

### Removal of radicals to form algebraic expression

2  Asked on August 1, 2021

### LinearOptimization duality

1  Asked on August 1, 2021

### Rotate 120-Cell Animation

0  Asked on August 1, 2021

### Determinant of Gauss (Area of a polygon of n vertices) how to implement

1  Asked on August 1, 2021

### Making boundary condition a variable in ParametricNDSolve

1  Asked on July 31, 2021

### Prevent Full-Screen on Presentation Start

0  Asked on July 31, 2021

### WolframScript not simplifying output

1  Asked on July 31, 2021 by charlottec

### How can I check if one expression implies another?

2  Asked on July 31, 2021 by user12404

### sorting out the coordinate and respective solution values for plotting

1  Asked on July 31, 2021

### How to understand the usage of Inner and Outer figuratively?

5  Asked on July 30, 2021 by xyz

### Nested lamdas: MapThreading a MapAt lambda

3  Asked on July 30, 2021

### Form sub-equations from a parent equation

0  Asked on July 30, 2021 by user75147

### Expressing a polynomial as a sum of squares

2  Asked on July 30, 2021 by pizzazz