Make graph using different x value compared to integral equation

1 次查看(过去 30 天)
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?

回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Symbolic Math Toolbox 的更多信息

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by