Make graph using different x value compared to integral equation
1 view (last 30 days)
Show older comments
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 Comments
Answers (0)
See Also
Categories
Find more on Symbolic Math Toolbox in Help Center and File Exchange
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!