write me a method to plot the distance from the earth to enceladus over a year in matlab

main.m
function plot_enceladus_distance()
% Plot the distance from the Earth to Enceladus over a year

% Set the JPL HORIZONS query parameters
target = 'Enceladus';
observer = 'Earth';
time_span = '1y'; % time span is 1 year
time_step = '1d'; % time step is 1 day
quantity = '3'; % quantity is distance

% Construct the JPL HORIZONS query URL
url = sprintf('https://ssd.jpl.nasa.gov/horizons_batch.cgi?batch=1&COMMAND="%s;%s"&MAKE_EPHEM=YES&TABLE_TYPE=OBSERVER&START_TIME=2022-01-01&STOP_TIME=2023-01-01&STEP_SIZE=%s&QUANTITIES=%s&CENTER=%s', target, observer, time_step, quantity, observer);

% Send the JPL HORIZONS query and parse the results
options = weboptions('ContentType', 'text');
data = webread(url, options);
lines = splitlines(data);
data_start = find(contains(lines, 'Date__(UT)__HR:MN:SC.fff'));
data_end = find(contains(lines, '$$EOE'));
data_lines = lines(data_start+2:data_end-1);
data_cells = cellfun(@(x) strsplit(x, ','), data_lines, 'UniformOutput', false); % parse the CSV data
data_matrix = [data_cells{:}];
time_strings = data_matrix(:,1);
distances = str2double(data_matrix(:,2));
times = datetime(time_strings, 'InputFormat', 'yyyy-MM-dd HH:mm:ss.SSS', 'TimeZone', 'UTC');

% Plot the distance data
plot(times, distances)
title('Distance from Earth to Enceladus over a year')
xlabel('Time (UTC)')
ylabel('Distance (km)')
end
1345 chars
33 lines