Celem zadania jest obliczenie całki z funkcji, która jest zadana na punktach w układzie współrzędnych biegunowych (r,θ). 
Całkę mozna obliczyć np. znajdując płaszczyznę przechodzącą przez 3 punkty i całkując wynikową funkcję symbolicznie. 
Poniżej przykłady siatek.

NT=4 (* ilość ,,promieni'' w zakresie 0<=θ≤Pi/2, minimum 3 *)

NR=4 (* ilość odcinków na które jest dzielony każdy promień, minimum 2 *)

R =Table[RandomInteger[{1,3}],{i,1,NT}] (* maksymalne wartości promienia, dla przykładu losowe *)

r=Table[j/NR*R[[i]],{i,1,NT},{j,1,NR}] (* wsp. radialne *)

θ=Table[Pi/2*i/(NT-1),{i,0,NT-1}] (* wsp. katowe *)

(* Przypominam, ze x= r cos(θ), y = r sin(θ) *)

points=Join[{{0,0}},Flatten[Table[{r[[i,j]]*Sin[θ[[i]]],r[[i,j]]*Cos[θ[[i]]]},{i,1,NT},{j,1,NR}],1]] (* punkty, razem ze środkiem ukł. współrzędnych *)

grid=Join[{{0,0,1}},Flatten[Table[{r[[i,j]]*Sin[θ[[i]]],r[[i,j]]*Cos[θ[[i]]],1},{i,1,NT},{j,1,NR}],1]] (* siatka *)

ListPlot[points,PlotStyle->PointSize[Large]] (* rozkład punktów na płaszczyznie XY *)

ListContourPlot[grid,Mesh->All] (* automatyczna triangulacja *)