How can I the show the solution/results in table form?

Question

Christoppe 2022-8-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1772630-how-can-i-the-show-the-solution-results-in-table-form

评论： Benjamin Kraus 2022-8-2

Hi! I'm doing a Simpsons 1/3 rule Numerical Methods, How can I show the results each solution in table form before the answer?

Example:

x f(x) (Odd or Even) f(x)*(Odd or Even)

-2.00000 0.05399 1.00000 0.05399

-1.55000 0.12001 4.00000 0.48004

-1.10000 0.21785 2.00000 0.43570

-0.65000 0.32297 4.00000 1.29188

-0.20000 0.39104 2.00000 0.78208

0.25000 0.38667 4.00000 1.54668

0.70000 0.31225 2.00000 0.62450

1.15000 0.20594 4.00000 0.82376

1.60000 0.11092 1.00000 0.11092

summation of [f(xi)* (odd or even)] 6.14955

Answer = Sq. Unit

Here;s my code:

eqn= 0.2+(25*x)-(200*x^2)+(675*x^3)-(900*x^4)+(400*x^5)

a= 0

b= 0.8

n= 8

clear
clc
format short
P=input ('Enter the Equation: ','s');
f=inline(P)
a=input ('Enter the Lower Value Limit: ')
b=input ('Enter the Upper Value Limit: ')
disp('Number of Segments "n" should be Divisible by 2!')
n=input ('Enter the Number of Segments: ')
h=round((b-a)/n,4)
odd=0;
even=0;
for i=1:2:n-1
    x=round(a+i*h,4);
    odd=round(odd+f(x),4);
end
for i=2:2:n-2
    x=round(a+i*h,4);
    even=round(even+f(x),4);
end
interc(i,:)=[x,f(x)]
I=round((h/3)*(f(a)+(4*odd)+(2*even)+f(b)),4);
fprintf('Integrated Value is %0.4f Sq. Unit', I)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Benjamin Kraus 2022-8-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1772630-how-can-i-the-show-the-solution-results-in-table-form#answer_1019770

在 MATLAB Online 中打开

There is still some polish you probably need to make, but here is one approach to displaying your results in a table:

a = 0;
b = 0.8;
n = 8;
f = @(x) 0.2+(25*x)-(200*x^2)+(675*x^3)-(900*x^4)+(400*x^5);
h = round((b-a)/n,4);
odd = 0;
even = 0;
oddOrEven = NaN(n-1,1);
x = NaN(n-1,1);
interc = NaN(n-1,1);
for i = 1:2:n-1
    x(i) = round(a+i*h,4);
    interc(i) = f(x(i));
    odd = round(odd+interc(i),4);
    oddOrEven(i) = odd;
end
for i = 2:2:n-2
    x(i) = round(a+i*h,4);
    interc(i) = f(x(i));
    even = round(even+interc(i),4);
    oddOrEven(i) = even;
end
I = round((h/3)*(f(a)+(4*odd)+(2*even)+f(b)),4);
results = table(x, interc, oddOrEven, interc.*oddOrEven, ...
    'VariableNames',{'x','f(x)','(Odd or Even)','f(x)*(Odd or Even)'});
disp(results)
     x     f(x)     (Odd or Even)    f(x)*(Odd or Even)
    ___    _____    _____________    __________________

    0.1    1.289        1.289              1.6615      
    0.2    1.288        1.288              1.6589      
    0.3    1.607        2.896              4.6539      
    0.4    2.456        3.744              9.1953      
    0.5    3.325        6.221              20.685      
    0.6    3.464        7.208              24.969      
    0.7    2.363        8.584              20.284      
fprintf('Integrated Value is %0.4f Sq. Unit\n', I)
Integrated Value is 1.6395 Sq. Unit

2 个评论
显示无隐藏无

Christoppe 2022-8-2

编辑：Christoppe 2022-8-2

Hello Sir! I appreciate your prompt response. May I ask is it possible to change the value of Odd or Even. Like example First and Last Rows will be multiplied by 1, All Odd rows will be multiplied by 4, and All Even rows will be multiplied by 2...

Benjamin Kraus 2022-8-2