We have investigated the solar activity variation with period shorter than 1000 days, through Fourier transformation of solar cycle 21 and 22 data. And real time predictions of the flare maximum intensity have been made by multilinear regression method to allow the use of multivariate vectors of sunspot groups or active region characteristics. In addition, we have examined the evolution of magnetic field and current density in active regions at times before and after flare occurrence, to check short term variability of solar activity. According to our results of calculation, solar activity changes with periods of 27.1, 28.0, 52.1, 156.3, 333.3 days for solar cycle 21 and of 26.5, 27.1, 28.9, 54.1, 154, 176.7, 384.6 days for solar cycle 22. Periodic components of about 27, 28, 53, 155 days are found simultaneously at all of two solar cycles. Finally, from our intensive analysis of solar activity data for three different terms of 1977~1982,1975~1998,and1978~1982, we find out that our predictions coincide with observations at hit rate of 76%,63% 76%,63% , 59 respectively.