Make graph using different x value compared to integral equation
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
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 Kommentare
Antworten (0)
Siehe auch
Kategorien
Mehr zu Symbolic Math Toolbox finden Sie in Help Center und File Exchange
Produkte
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!