Greece and the surrounding area are characterized by high seismicity. Strong earthquakes (M>6.5) have repeatedly occurred in this area, as historical information and instrumental recording reveal. Using as input data 67 strong earthquakes that occurred since the beginning of the 20th century, the coseismic stress changes are calculated, in order to make a probabilistic earthquake forecast in the study area under the influence of past events. For this purpose, the calculated coseismic stress changes are translated into earthquake probability using an earthquake nucleation constitutive relation, which includes permanent and transient effect of the sudden stress changes. According to this method, a sudden change in stress seems to modify earthquake rate, moving other faults toward or away from failure, changing the probability of potential earthquake on these faults. Earthquake probability on a fault is lowest after the last event but as tectonic stress grows the odds of another earthquake increase. For all needed calculations a pdf (probability density function) for the time of failure for an earthquake of defined magnitude on the fault of interest must be taken into account along with the calculated stress changes on the fault. Specifically, the estimated probability values concern the probability in each part of a given fault or fault segment, and the probability distribution is illustrated across the specific fault. All calculations were performed at 10 km depth but it was necessary to check whether the estimated probability values vary with depth. Therefore, all estimations were performed for each fault or fault segment at the depth of 8, 12 and 15 km. The probability calculations were carried out and given for the whole study area during the next 30 years in the form of tables and maps.