Error in using 'cumtrapz' for integration
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Hello
I need to use 'cumtrapz' for a particular part of my project, but before that i am trying to understand how it works for multiple time cumulative integration. As a test, i was trying the following basic code:
ny = 101;
eta = linspace(-1,1,ny)';
f1 = sin(eta);
% Calculation of first cumulative integral
g1 = cumtrapz(eta,f1);
% Calculation of second cumulative integral
h1 = cumtrapz(eta,g1);
% Calculation of third cumulative intergal
k1 = cumtrapz(eta,h1);
We know that first integration of sinx is -cosx and then second intergration is -sinx.
So the value of g1 at each grid point must match {-cos(grid point) + cos(-1)}, Which it does.
But then the value of h1 at each grid point must match {-sin(grid point) + sin(-1)}, which it doesn't.
I was wondering why this happens.
0 commentaires
Réponse acceptée
Matt J
le 10 Août 2022
Modifié(e) : Matt J
le 10 Août 2022
But next, value of h1 at each grid point must match {-sin(grid point value) - (-sin(-1))},
No, it will be -sin(x)+x*cos(-1) -( -sin(-1)-1*cos(-1))
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Numerical Integration and Differentiation 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!