# Numerical integration using Trapz

8 visualizzazioni (ultimi 30 giorni)
Waseem Akhtar il 6 Mag 2021
Commentato: Star Strider il 7 Mag 2021
Hi all,
I am trying to calculate a double integral over a rectangular grid using trapz.
The grid is defined as:
[x,z] = meshgrid(xp,zp);
xp and zp are the 1d arrays of x and z coordinates of the grid.
I trying to calculate the function values (F) at the grid nodes as:
F(i,j)=trapz(zp,trapz(xp,g.*f,2));
here,
g is a function of x and z i.e. g(x,z) whereas, f values are already calculated values at each node of the grid i.e. not a function of x and z.
Can someone suggest how may I do it, please.
##### 0 CommentiMostra -2 commenti meno recentiNascondi -2 commenti meno recenti

Accedi per commentare.

### Risposta accettata

Star Strider il 6 Mag 2021
That should produce one scalar value as ‘F’, so the subscripts are not necessary.
I would do something like this —
F = trapz(zp,trapz(xp,g.*f));
The result of the first trapz call will be a vector, and of the second will be a scalar. So long as ‘g’ and ‘f’ have compatible dimensions, and they match ‘xp’ and ‘zp’ appropriately, that should work.
##### 6 CommentiMostra 4 commenti meno recentiNascondi 4 commenti meno recenti
Waseem Akhtar il 7 Mag 2021
Thank you Star Strider, your first answer was correct.
Star Strider il 7 Mag 2021
As always, my pleasure!

Accedi per commentare.

### Categorie

Scopri di più su Numerical Integration and Differentiation in Help Center e File Exchange

### Community Treasure Hunt

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

Start Hunting!

Translated by