Earthquakes are a phenomenon that is still a mystery in terms of predicting events, one of which is the magnitude. As technology develops, there are many algorithms that can be used as approaches in earthquake forecasting. In the context of magnitude forecasting, the application of GaussianNB, Random Forest, and SVM has the potential to reveal these patterns and relationships in the data. With the six main phases of this research, namely data acquisition, data pre-processing, feature selection, model training, forecast result evaluation, and performance analysis, this study is expected to contribute to the development of more accurate and effective earthquake forecasting methods. From these results we first obtain the result that the GaussianNB model has a relatively simple and fast method in training its model. However, the weakness lies in the assumption of a Gaussian distribution, which may not always suit the complex and diverse characteristics of earthquake data. Second, Random Forest, this method can increase accuracy and overcome the overfitting problem that occurs when forecasting magnitudes. In contrast to GaussianNB, it tends to result in models with greater complexity and requires more time to compute. The third option is SVM, which has both benefits and drawbacks that must be taken into account. The capacity of SVM to separate data that has both linear and nonlinear separation is one of its key advantages; nevertheless, the main drawback is that it is sensitive to hyperparameter adjustments.