  Results of analysis of 23 events of the 2014-2015 mutual event series from
the Vainu Bappu Observatory are presented. Our intensity distribution model for
the eclipsed/occulted satellite is based on the criterion that it simulates a
rotational light curve that matches the ground based light curve. Dichotomy in
the scattering characteristics of the leading and trailing sides explain the
basic shape of the rotational light curves of Europa, Ganymede and Callisto. In
case of Io the albedo map from USGS along with global values of scattering
parameters work well. Mean values of residuals in $(O-C)$ along and
perpendicular to the track are found to be -3.3 and -3.4 mas respectively
compared to "L2" theory for the seven 2E1/2O1 events. The corresponding R.M.S
values are 8.7 and 7.8 mas respectively. For the five 1E3/1O3 events, the along
and perpendicular to the track mean residuals are 5.6 and 3.2 mas respectively.
The corresponding R.M.S. residuals are 6.8 and 10.5 mas respectively. We
compare the results using the chosen model (Model 1) with a uniform but limb
darkened disk (Model 2). The residuals with Model 2 of the 2E1/2O1 and 1E3/1O3
events indicate a bias along the satellite track. The extent and direction of
bias is consistent with the shift of the light center from the geometric
center. Results using Model 1, which intrinsically takes into account the
intensity distribution show no such bias.
