Bonjour mon ami !
Ta demande date d'il y a une semaine mais au cas où tu n'aurais toujours pas de réponse, voici une manière de faire, probablement pas la meilleure mais peut-être cela te suffira-t-il ;-)
L'astuce qui me vient à l'esprit est de forcer ta fonction à zéro sur le triangle supérieur (0<y<x<30) sur lequel tu n'intègre pas. Cela te permet d'utiliser sans crainte la fonction dblquad sur le carré [0,30]².
Exemple :
f=@(x,y)(x+y);
intf=dblquad(@(x,y)((x<=30-y).*f(x,y)),0,30,0,30);
Je te souhaite une bien bonne journée !
Tous les animaux criaient bien haut
Qu’il était le crapaud le plus beau, quand il jouait du banjo