Why is my MATLAB plot wrong?

Hi! I just started my Matlab course and i have run into a little problem which i hope you can help me with. I am suppose to write a function which integrates a function e^(-t^2) from 0 to some x which is a variable in the function. I did that like this:
function y = simpson_method(x)
y1 = exp(-0.^2);
y2 = exp(-(x).^2);
c = 0;
d = 0;
incr = 0.1;
for i = 1 : (x - 0)/incr - 0.5
y = exp(-(0 + i*incr).^2);
if i/2 == floor(i/2)
d = d + y;
continue
else
c = c + y;
continue
end
end
y = incr/3 * (y1 + 4*c + 2*d + y2);
The function integrates correctly and gives me a desired output for every x. But the problem is when i try to plot it:
x = 0:0.1:15;
y = simpson_method(x);
plot(x,y)
It gives me the wrong graph which i find weird because all the outputs are correct. I get this which is nonsense. Whats happening?
%

5 Kommentare

Adam
Adam am 18 Mai 2018
What is the right graph?
And what are the values of y that you get out?
Edward Lisa
Edward Lisa am 18 Mai 2018
This is the right plot and my function outputs y as very close to 0 when x=0 and reaches 0.884.. when x goes larger which is exactly the plot of the function when integrated on https://www.integral-calculator.com/ with input e^(-x^2) from 0 to your chosen upper limit.
Stephen23
Stephen23 am 18 Mai 2018
@Edward Lisa: can you please confirm the ambiguous syntax:
  1. e^((-t)^2)
  2. e^(-(t^2))
Edward Lisa
Edward Lisa am 18 Mai 2018
Its e^(-(t^2))
Adam
Adam am 21 Mai 2018
If your y values are correct I can see no reason at all why plot would not give the correct answer so your y values cannot be correct if you are literally just calling plot.

Melden Sie sich an, um zu kommentieren.

Antworten (1)

KALYAN ACHARJYA
KALYAN ACHARJYA am 21 Mai 2018
Bearbeitet: KALYAN ACHARJYA am 21 Mai 2018

0 Stimmen

Shown Graph the x ranges are different

Kategorien

Mehr zu Programming finden Sie in Hilfe-Center und File Exchange

Produkte

Version

R2017b

Gefragt:

am 18 Mai 2018

Bearbeitet:

am 21 Mai 2018

Community Treasure Hunt

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

Start Hunting!

Translated by