How to find the sum of values of a function at mid point of every subinterval?
3 ビュー (過去 30 日間)
古いコメントを表示
I am trying to run this code but don't know how to write code for finding value of function F at mid point of subinterval say, if my interval is [ti, ti+1] for i=0,1,2, n. and need to find sum of functional value at every (ti +ti+1)/2. I need to do it for Rsum2 in this code.
clc; clear all; format long
%%%%%%%%%%%%%%%%%%%%%%%%%%
alpha=0.8;% fractional index
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
TN=1 % time
N=10
T0=0
tau=TN/N
T=[T0:tau:TN]
X0=0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
X(1)=X0
F =@(X,T)=X+T
for n=0:N-1
Rsum=0;
for j=1:1:n+2
Rsum=Rsum+2*F (j)
end
Rsum2=0;
for j=1:1:n+1
Rsum2=Rsum2+4*F (j+1/2)
end
X(n+2)=X(1) +(tau^alpha)*Rsum + (tau^alpha)*Rsum2
end
0 件のコメント
採用された回答
Matt J
2021 年 3 月 20 日
編集済み: Matt J
2021 年 3 月 20 日
If your function is vectorized, like in the following example, it is quite simple:
fun=@(tm) tm.^2+ sqrt(tm); %A vectorized function
n=8;
t=sort(rand(1,n)) %interval end points t(i)
tmid = t(1:end-1)/2 +t(2:end)/2 %interval mid-points
result = sum(fun(tmid))
3 件のコメント
その他の回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Calculus についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!