Цитата
у меня получается что div(a)=0, такое возможно или я что то нитак вычислила.

Всё верно, дивергенция равна нулю. Это из задачи с соседней ветки? Там требуется найти поток векторного поля только через одну поверхность, а не всю пирамиду...