Direktivitet i Matlab

Hej,

Jeg skal beregne direktiviteten for en antennegruppe, hvor punkterne findes i en matrix med varierende theta vandret og varierende phi lodret.
Spørgsmål går på, hvordan man regner et dobbelt-integrale ud fra disse punkter. Skal jeg virkelig oprette en funktion, og bruge "dblquad"? Det passer nemlig meget dårligt sammen med resten af koden...

Håber at der er nogen, der har prøvet noget lignende.