How to find the sum of values of a function at mid point of every subinterval?
3 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Mubashara Wali
le 20 Mar 2021
Commenté : Mubashara Wali
le 22 Mar 2021
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 commentaires
Réponse acceptée
Matt J
le 20 Mar 2021
Modifié(e) : Matt J
le 20 Mar 2021
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 commentaires
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Calculus dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!