# How can I plot integral of this equation?

1 view (last 30 days)
E kara on 29 Jul 2021
Commented: Star Strider on 29 Jul 2021
I couldn't manage to plot this equation over time.
I wrote the formula for integral, but these are scalar, but i had to wrote '.*' instead of *
I should have a decreasing trend, but there is no plot. Iam new in Matlab, could you please let me know where is the my mistakes?
clc;
figure('color','white');
while(1)
sigma.phi=input('Enter 0.1 in a.u.: ');
mu=input('Enter 2 in mm-1: ');
n=input('number of evicted electrons: ');
for t=0:300
fun=@(x)(n).*(exp(-((mu).*x))).*(exp(-((sigma.phi).*t.*(exp(-(mu).*x)))));
Q=integral(fun,0,1);
end
plot(Q,[0,300]);
xlabel('Time / s ');
title('Signal');
ylabel('Intensity');
hold all;
ch= input('Press 1 to continue and 0 to exit: ')
if ch == 0
break;
end
end

Star Strider on 29 Jul 2021
Try something like this —
% clc;
figure('color','white');
while(1)
sigma.phi=input('Enter 0.1 in a.u.: ');
mu=input('Enter 2 in mm-1: ');
n=input('number of evicted electrons: ');
t=0:300;
fun=@(x,t)(n).*(exp(-((mu).*x))).*(exp(-((sigma.phi).*t.*(exp(-(mu).*x)))));
for k = 1:numel(t)
Q(k)=integral(@(x)fun(x,t(k)),0,1);
end
plot(t,Q);
xlabel('Time / s ');
title('Signal');
ylabel('Intensity');
hold all;
ch= input('Press 1 to continue and 0 to exit: ')
if ch == 0
break;
end
end
.
##### 2 CommentsShowHide 1 older comment
Star Strider on 29 Jul 2021
As always, my pleasure!
.