The sunspot drawings of Johann Staudacher of 1749--1799 were used to determine the solar differential rotation in that period. These drawings of the full disk lack any indication of their orientation. We used a Bayesian estimator to obtain the position angles of the drawings, the corresponding heliographic spot positions, a time offset between the drawings and the differential rotation parameter \delta\Omega, assuming the equatorial rotation period is the same as today. The drawings are grouped in pairs, and the resulting marginal distributions for \delta\Omega were multiplied. We obtain \delta\Omega=-0.048 \pm 0.025 d^-1 (-2.75^o/d) for the entire period. There is no significant difference to the value of the present Sun. We find an (insignificant) indication for a change of \delta\Omega throughout the observing period from strong differential rotation, \delta\Omega\approx -0.07 d^-1, to weaker differential rotation, \delta\Omega\approx-0.04 d^-1.

From high-precision MOST photometry spanning 35 days the existence of two spots rotating with slightly differing periods is confirmed. From the marginal probability distribution of the derived differential rotation parameter k its expectation value as well as confidence limits are computed directly from the data. The result depends on the assumed range in inclination i, not on the shape of the prior distributions. Two cases have been considered: (a) The priors for angles, inclination i of the star and spot latitudes beta_1,2, are assumed to be constant over i, beta_1, and beta_2. (b) The priors are assumed to be constant over cosi, sin beta_1, and sin beta_2. In both cases the full range of inclination is considered: 0^o< = i< = 90^o. Scale-free parameters, i.e. periods and spot areas (in case of small spots) are taken logarithmically. Irrespective of the shape of the prior, k is restricted to 0.03 < = k < = 0.10 (one-sigma limits). The inclination i of the star is photometrically ill-defined.