A muscle model is described that uses electromyogram (EMG), muscle length and speed of contraction to predict muscle force. Physiological parameters are the Hill constants and the shape of the twitch response to a single stimulus. The model was incorporated in a jaw model of the rabbit and tested by predicting the bite force produced by the jaw muscles during mastication. The time course of the calculated force appeared to match the bite force, measured in vivo by a strain gauge, applied to the bone below the teeth. The variation in peak strain amplitude from cycle to cycle correlated with the variation predicted by the model. The peak amplitude of the integrated EMGs of individual jaw muscles showed an average correlation with peak strain of 0.41. Use of the sum of the available peak amplitudes, weighted according to their effect upon the bite force increased the correlation to 0.46; the model predicted bite forces showed a correlation of 0.57 with the strain. The increase in correlation was statistically significant. The muscle forces were calculated using a minimum number of easily obtainable constants.