I suggest the following script :
clearvars; clc; close all;
syms x y a b h k
% ellipse axis
a=2; b=1;
%ellipse center
h=1; k=0;
%ellipse equation
ellipse= (((x-h)^2)/(a^2))+(((y-k)^2)/(b^2))==1;
%plot the ellipse
plotZoom=(max(a,b)+max(abs(h), abs(k)))+1;
fimplicit(ellipse, [-plotZoom plotZoom]);hold on;
plot([-plotZoom plotZoom], [0 0], '-k');
plot([0 0 ], [-plotZoom plotZoom], '-k');
plot(h, k, 'b*');
xlabel('x'); ylabel('y');
grid on; zoom on;
In case you execute it , you will get this:
