Um nicht ewig um den Brei zu labern - eine Lösung könnte so aussehen:
A = Erdradius^2 * (Pi/36) * |sin(alpha1) - sin(alpha2)|
wobei die alpha jeweils die Breitengrade im Bogenmaß (Radiant) sind, die das Gebiet begrenzen. Die Betragsbildung ist sinnvoll, denn ansonsten kann bei sorgloser Eingabe der Winkelwerte das Ergebnis negativ werden.
Die Genauigkeit ist zu steigern, wenn du den Erdradius an das jeweilige Gebiet anpaßt. Für eine Kugel wäre die Lösung exakt.