If you know what is the feeding function of the array element and only interested to see the radiation behaviour, then the suggested resources by other experts are usefull. You can also write your own program to find the radiation. it is not that much difficult.
However, if you have a desired radiation pattern in mind and eager to find the appropriate feeding of the array element (array synthesis), then it is another stroy.
Normally the classic ways of feeding would not be very effecient (Dolfchebychef, fourier, ...). In this case you should look at the global optimization techniuqes.