Problem with a program concerning composite trapezoid rule

4 次查看(过去 30 天)
In this program I'm using the composite trapezoid rule on integral dx/(1+x^2) from 1 to 0. I've used the same program with different functions and haven't had issues till I tried it with this function. It says I have an error with line y=f(x); Here is code:
f = @(x) 1/(1+x^2);
a=0;
b=1;
n=100;
format long
h = (b-a)/n;
S = f(a);
i=1:1:n-1;
x=a+h*i;
y=f(x);
S=S+2*sum(y);
S=S+f(b);
F=h*S/2

采纳的回答

sixwwwwww
sixwwwwww 2013-10-22
Dear Drake, re-write first line as:
f = @(x) 1./(1+x.^2);
then it should work. Good luck!

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by