Why is my MATLAB's bode plot wildly off?
I suspect that the reason for the difference is the number of samples used in the plot, since the extreme points are likely going to infinity (or very large numbers) such that if less points are used, those large values are simply missed in the computation. If the OP uses the same number of points in each, then the same result should be achieved- and for all the points shown if they were super-imposed on each plot it would be clear that for the points given, the results are the same.
My first suspicion before the OP added the plots was that it was due to the units of the axis. That is not the case but leaving this part of my answer below as it is another common reason for Bode plot differences:
The units of the frequency axis should be confirmed as they could be in units of:
Continuous domain units of frequency:
Hertz, or cycles/sec
radians/sec
Discrete domain units of (normalized) frequency:
cycles per sample (unique domain is $[-1,1)$)
radians/sample (unique domain is $[-\pi, \pi)$)