function energy()
% Clears command window, closes all figures, and clears workspace
clc;
close all;
clear;
% Define the speed of light (c) in m/s
c = 2.9979e8;
% Input: Mass vector from user
m = input('Enter a vector of masses in kg: ');
% Ensure the input is a valid numeric vector
if ~isnumeric(m) || isempty(m)
error('Input must be a numeric vector.');
end
% Compute energy values using E = mc^2
E = m .* (c^2);
% Display results
disp('Mass (kg) Energy (J)');
disp([m(:), E(:)]);
% Save results to an Excel file
filename = 'Energy_Data.xlsx';
data = [m(:), E(:)];
headers = {'Mass (kg)', 'Energy (J)'};
% Write headers and data to an Excel file
writematrix(headers, filename, 'Sheet', 1, 'Range', 'A1');
writematrix(data, filename, 'Sheet', 1, 'Range', 'A2');
% Plot energy vs mass (Standard plot)
figure(1);
plot(m, E, 'bo-', 'LineWidth', 2);
xlabel('Mass (kg)');
ylabel('Energy (J)');
title('Energy vs Mass (Standard Plot)');
grid on;
% Logarithmic plots in a single figure using subplots
figure(2);
subplot(2,2,1);
semilogx(m, E, 'r*-');
xlabel('Mass (kg)');
ylabel('Energy (J)');
title('Semilogx: Log Scale on X-axis');
grid on;
subplot(2,2,2);
semilogy(m, E, 'g*-');
xlabel('Mass (kg)');
ylabel('Energy (J)');
title('Semilogy: Log Scale on Y-axis');
grid on;
subplot(2,2,3);
loglog(m, E, 'm*-');
xlabel('Mass (kg)');
ylabel('Energy (J)');
title('Loglog: Log Scale on Both Axes');
grid on;
% Adjust subplot layout for better visualization
subplot(2,2,4);
text(0.1, 0.5, 'Logarithmic Analysis of Energy vs Mass', 'FontSize', 12, 'FontWeight', 'bold');
axis off;
% Notify user about file storage
fprintf('Energy values have been saved to %s\n', filename
);
% Sample Output
fprintf('Enter a vector of masses in kg: [1 5 10 50 100]\n'); fprintf('\nMass (kg) Energy (J)\n'); sample_m = [1 5 10 50 100];
sample_E = sample_m .* (c^2);
fprintf('%-12.4f %-12.4e\n', [sample_m
; sample_E
]); fprintf('\nEnergy values have been saved to Energy_Data.xlsx\n'); end
ZnVuY3Rpb24gZW5lcmd5KCkKICAgICUgQ2xlYXJzIGNvbW1hbmQgd2luZG93LCBjbG9zZXMgYWxsIGZpZ3VyZXMsIGFuZCBjbGVhcnMgd29ya3NwYWNlCiAgICBjbGM7CiAgICBjbG9zZSBhbGw7CiAgICBjbGVhcjsKICAgIAogICAgJSBEZWZpbmUgdGhlIHNwZWVkIG9mIGxpZ2h0IChjKSBpbiBtL3MKICAgIGMgPSAyLjk5NzllODsgCiAgICAKICAgICUgSW5wdXQ6IE1hc3MgdmVjdG9yIGZyb20gdXNlcgogICAgbSA9IGlucHV0KCdFbnRlciBhIHZlY3RvciBvZiBtYXNzZXMgaW4ga2c6ICcpOwogICAgCiAgICAlIEVuc3VyZSB0aGUgaW5wdXQgaXMgYSB2YWxpZCBudW1lcmljIHZlY3RvcgogICAgaWYgfmlzbnVtZXJpYyhtKSB8fCBpc2VtcHR5KG0pCiAgICAgICAgZXJyb3IoJ0lucHV0IG11c3QgYmUgYSBudW1lcmljIHZlY3Rvci4nKTsKICAgIGVuZAogICAgCiAgICAlIENvbXB1dGUgZW5lcmd5IHZhbHVlcyB1c2luZyBFID0gbWNeMgogICAgRSA9IG0gLiogKGNeMik7CiAgICAKICAgICUgRGlzcGxheSByZXN1bHRzCiAgICBkaXNwKCdNYXNzIChrZykgICAgICBFbmVyZ3kgKEopJyk7CiAgICBkaXNwKFttKDopLCBFKDopXSk7CiAgICAKICAgICUgU2F2ZSByZXN1bHRzIHRvIGFuIEV4Y2VsIGZpbGUKICAgIGZpbGVuYW1lID0gJ0VuZXJneV9EYXRhLnhsc3gnOwogICAgZGF0YSA9IFttKDopLCBFKDopXTsKICAgIGhlYWRlcnMgPSB7J01hc3MgKGtnKScsICdFbmVyZ3kgKEopJ307CiAgICAKICAgICUgV3JpdGUgaGVhZGVycyBhbmQgZGF0YSB0byBhbiBFeGNlbCBmaWxlCiAgICB3cml0ZW1hdHJpeChoZWFkZXJzLCBmaWxlbmFtZSwgJ1NoZWV0JywgMSwgJ1JhbmdlJywgJ0ExJyk7CiAgICB3cml0ZW1hdHJpeChkYXRhLCBmaWxlbmFtZSwgJ1NoZWV0JywgMSwgJ1JhbmdlJywgJ0EyJyk7CiAgICAKICAgICUgUGxvdCBlbmVyZ3kgdnMgbWFzcyAoU3RhbmRhcmQgcGxvdCkKICAgIGZpZ3VyZSgxKTsKICAgIHBsb3QobSwgRSwgJ2JvLScsICdMaW5lV2lkdGgnLCAyKTsKICAgIHhsYWJlbCgnTWFzcyAoa2cpJyk7CiAgICB5bGFiZWwoJ0VuZXJneSAoSiknKTsKICAgIHRpdGxlKCdFbmVyZ3kgdnMgTWFzcyAoU3RhbmRhcmQgUGxvdCknKTsKICAgIGdyaWQgb247CiAgICAKICAgICUgTG9nYXJpdGhtaWMgcGxvdHMgaW4gYSBzaW5nbGUgZmlndXJlIHVzaW5nIHN1YnBsb3RzCiAgICBmaWd1cmUoMik7CiAgICAKICAgIHN1YnBsb3QoMiwyLDEpOwogICAgc2VtaWxvZ3gobSwgRSwgJ3IqLScpOwogICAgeGxhYmVsKCdNYXNzIChrZyknKTsKICAgIHlsYWJlbCgnRW5lcmd5IChKKScpOwogICAgdGl0bGUoJ1NlbWlsb2d4OiBMb2cgU2NhbGUgb24gWC1heGlzJyk7CiAgICBncmlkIG9uOwogICAgCiAgICBzdWJwbG90KDIsMiwyKTsKICAgIHNlbWlsb2d5KG0sIEUsICdnKi0nKTsKICAgIHhsYWJlbCgnTWFzcyAoa2cpJyk7CiAgICB5bGFiZWwoJ0VuZXJneSAoSiknKTsKICAgIHRpdGxlKCdTZW1pbG9neTogTG9nIFNjYWxlIG9uIFktYXhpcycpOwogICAgZ3JpZCBvbjsKICAgIAogICAgc3VicGxvdCgyLDIsMyk7CiAgICBsb2dsb2cobSwgRSwgJ20qLScpOwogICAgeGxhYmVsKCdNYXNzIChrZyknKTsKICAgIHlsYWJlbCgnRW5lcmd5IChKKScpOwogICAgdGl0bGUoJ0xvZ2xvZzogTG9nIFNjYWxlIG9uIEJvdGggQXhlcycpOwogICAgZ3JpZCBvbjsKICAgIAogICAgJSBBZGp1c3Qgc3VicGxvdCBsYXlvdXQgZm9yIGJldHRlciB2aXN1YWxpemF0aW9uCiAgICBzdWJwbG90KDIsMiw0KTsKICAgIHRleHQoMC4xLCAwLjUsICdMb2dhcml0aG1pYyBBbmFseXNpcyBvZiBFbmVyZ3kgdnMgTWFzcycsICdGb250U2l6ZScsIDEyLCAnRm9udFdlaWdodCcsICdib2xkJyk7CiAgICBheGlzIG9mZjsKICAgIAogICAgJSBOb3RpZnkgdXNlciBhYm91dCBmaWxlIHN0b3JhZ2UKICAgIGZwcmludGYoJ0VuZXJneSB2YWx1ZXMgaGF2ZSBiZWVuIHNhdmVkIHRvICVzXG4nLCBmaWxlbmFtZSk7CiAgICAKICAgICUgU2FtcGxlIE91dHB1dAogICAgZnByaW50ZignXG5TYW1wbGUgT3V0cHV0OlxuJyk7CiAgICBmcHJpbnRmKCdFbnRlciBhIHZlY3RvciBvZiBtYXNzZXMgaW4ga2c6IFsxIDUgMTAgNTAgMTAwXVxuJyk7CiAgICBmcHJpbnRmKCdcbk1hc3MgKGtnKSAgICAgIEVuZXJneSAoSilcbicpOwogICAgc2FtcGxlX20gPSBbMSA1IDEwIDUwIDEwMF07CiAgICBzYW1wbGVfRSA9IHNhbXBsZV9tIC4qIChjXjIpOwogICAgZnByaW50ZignJS0xMi40ZiAlLTEyLjRlXG4nLCBbc2FtcGxlX207IHNhbXBsZV9FXSk7CiAgICBmcHJpbnRmKCdcbkVuZXJneSB2YWx1ZXMgaGF2ZSBiZWVuIHNhdmVkIHRvIEVuZXJneV9EYXRhLnhsc3hcbicpOwplbmQK