The method is always a visible argument, never a hidden default. The same plant discretised by zero-order hold and by Tustin has different discrete poles, so a silent choice would make two correct-looking answers disagree with no way to see why.
Tustin, the bilinear transforms = (2/Ts)*(z-1)/(z+1) — what a filter designer wants. Maps the left half-plane onto the unit disc exactly, so stability is preserved.
Tustin, the bilinear transforms = (2/Ts)*(z-1)/(z+1) — what a filter designer wants. Maps the left half-plane onto the unit disc exactly, so stability is preserved.