Make graph using different x value compared to integral equation

syms x f

a = 1.27;

b = 0.19 ;

c = 0.44 ;

d = 0.2 ;

e = 300 ;

w = i*e*log(2);

f = 5*10^21 ;

g = 2*10^23 ;

h1= linspace (0.2, 1.7);

k1 = linspace (0, 1.5);

i = 8.61733*10^(-5);

j = 2*10^21;

B1 = 0.085 ;

A1 = j.*(f/g).^((i.*e)/(4*B1)).*((2*B1.^2)./((2*B1)-(i.*e))).*exp(-(1/(2*B1)).*(a-k1-((b.^2)./(4*B1))));

P1 = (1./sqrt(2*b.^2*pi)).*exp(-((x+w)+((b.^2)./(2*B1))-a).^2./(2*b.^2));

F1 = (2*exp((h1-(x+w))./(i.*e)))./(1+2*exp((h1-(x+w))./(i.*e))+exp(((2*h1)-(2*(x+w))-d)./(i.*e)));

D1 = A1.*((2./F1).^((i.*e)./(2*B1))).*P1 ;

P11 = (1./sqrt(2*b.^2*pi)).*exp(-((x-d-w)+((b.^2)./(2*B1))-a).^2./(2*b.^2));

F11 = (2*exp((h1-(x-d-w))./(i.*e)))./(1+2*exp((h1-(x-d-w))./(i.*e))+exp(((2*h1)-(2*(x-d-w))-d)./(i.*e)));

D11 = A1.*((2./F11).^((i.*e)./(2*B1))).*P11;

f = (A1.*((2./F1).^((i.*e)./(2*B1))).*P1)+(A1.*((2./F11).^((i.*e)./(2*B1))).*P11);

F(x) = int(f,x,[0 1.7]);

plot(F,[0 1.7],'k'); xlabel('x'), ylabel ('y'); xlim([0,1.7]); ylim([10^15,10^20]); set(gca, 'YScale','log');set(gca, 'XScale','linear');

I'd like to obtain different graph compared to integral's x: x into h1-k1.

How can I solve this coding problem?