MATLAB code for given equation

[pappu_sam]

MATLAB code

am trying to write matlab code for given equation. and am getting error.
please could anybody tell me how to solve this equation in matlab.
here in this equation bessel function is used.
am getting an error while integrating the bessel function.
and all others are constants as given in the code.
please check this code and please suggest me how to solve it.
thanks in advance.




clc
clear all
close all


ds = 200;
a = 0.025;
L = sqrt(ds * a);
Vw = 100;
Vs = 30000;
t = (L)/(Vw);
u = 55;
den = (7*(10^(-6)));
Uch = 13.8;
e = 1-((0.45 * Uch))/(u);
q0 = ((e*u*a))/(t);
k = 65;
Cp = 500;
E = (1.2*(10^5));
v = 0.3;
h = ((1.5)*(10^(-2)));
y = 300;
H = 10000;


b = ((Vw)*(den)*(Cp))/(k);

x = 1;
x1 = 0.001;
z = {(b/2)*(x-(x1))};
m = z;
p = exp(m);
f = Inlinefunction;
f(x1) = 0.5*{((q0)/((pi)*k))}*{1+((x1)/(a))}*{p}*{besselj(1,z)};
T = diff(fnval(fnint(f),[-a,a]));