I need help plotting this, i Cant figure out what is wrong

1 次查看（过去 30 天）
ragheed idrees
ragheed idrees 2019-6-5
评论： ragheed idrees 2019-6-6
fun = @(x) x.^3 - 10*x.^2 + 16*x +80;
a = 0;
b = 6;
for n = [2, 3, 4, 5, 6, 8, 10, 12, 14, 16, 18, 20]
I = A9P1TRAPIdrees(fun,a,b,n);
err = abs((372-I)/372)*100;
fprintf('n = %d value=%f error= %f\n', n, I, err);
end
subplot(3,1,1); %f(x) vs. x
fplot(fun,[0 6])
title({['Ragheed Idrees ',datestr(now)];'';'f(x) vs x)'})
xlabel('x')
ylabel('f(x)')
subplot(3,1,2); %Integral vs number of sugments
plot(n,I)
title('Integral vs number of sugments')
xlabel('number of segments')
ylabel('Integral')
subplot(3,1,3); %Percent error vs number of sugments
loglog(n,err)
title('Percent error vs number of sugments')
xlabel('number of segments')
ylabel('Percent Error')
-----------------------------------------------------------------------------A9P1TRAPIdrees-------------------------------------------------------------------------------------------------------------
function I = A9P1TRAPIdrees(fun,a,b,n)
%this function uses the trapezoidal rule to integrate a function with given
%boundaries and number of segments
%defining the output
h = (b-a)/n;
s=0.5 * (fun(a)+fun(b));
I=0;
for jj = 1 : n-1
s = s +fun(a + jj*h);
end
I = h *s;
it gives me a warning "Function behaves unexpectedly on array inputs. To improve performance, properly vectorize your function to return an output with the same size
and shape as the input arguments. "
and the last 2 subplots wont plot!

请先登录，再回答此问题。

采纳的回答

Walter Roberson
Walter Roberson 2019-6-6
fun = @(x) x^3 - 10*x.^2 + 16*x +80;
should be
fun = @(x) x.^3 - 10*x.^2 + 16*x +80;
  5 个评论
显示 3更早的评论
Walter Roberson
Walter Roberson 2019-6-6
nvals = [2, 3, 4, 5, 6, 8, 10, 12, 14, 16, 18, 20];
numn = length(nvals);
err = zeros(1, numn);
I = zeros(1, numn);
for nidx = 1 : numn
n = nvals(nidx);
I(nidx) = A9P1TRAPIdrees(fun,a,b,n);
err(nidx) = abs((372-I(nidx))/372)*100;
fprintf('n = %d value=%f error= %f\n', n, I(nidx), err(nidx));
end
subplot(3,1,1); %f(x) vs. x
fplot(fun,[0 6])
title({['Ragheed Idrees ',datestr(now)];'';'f(x) vs x)'})
xlabel('x')
ylabel('f(x)')
subplot(3,1,2); %Integral vs number of sugments
plot(nvals, I)
title('Integral vs number of sugments')
xlabel('number of segments')
ylabel('Integral')
subplot(3,1,3); %Percent error vs number of sugments
loglog(nvals, err)
title('Percent error vs number of sugments')
xlabel('number of segments')
ylabel('Percent Error')
-----------------------------------------------------------------------------A9P1TRAPIdrees-------------------------------------------------------------------------------------------------------------

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Line Plots 的更多信息

标签

产品


版本

R2019a

Community Treasure Hunt

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

Start Hunting!

Translated by