MATLAB: I am not able to get a line to appear on the graph

graph

clear all;

radius = 6378:50000; force = (3.986e5/(6378.^2));

figure plot(radius,force); axis([0 50e3 0 10]) xlabel('r(km)'); ylabel('F(m/s)');

Best Answer

You had force just as one single number, not a list of forces for each radius. Try this:

radius = 6378:50000;
force = 3.986e5 ./ (radius .^ 2);
figure
plot(radius,force, 'b-', 'LineWidth', 2);
% axis([0 50e3 0 10])
xlabel('r(km)');
ylabel('F(m/s)');
grid on;
title('force vs. radius', 'FontSize', 20);

Related Solutions

MATLAB: Colour of legend doesn’t change

Your a looks to have multiple columns, so S1 would have multiple columns.

plot(f,S1,'b');

would generate one line object per column.

MATLAB: Get x and y coordinates of motion from ode45

You've got a 4th order system (2 equations for x and y, each 2nd order --> 2 x 2 = 4), so your state vector will need to be 4 elements, not 2 elements like you currently have. Define the state vector as follows and then work from there:

y(1) = the x position
y(2) = the y position
y(3) = the x velocity
y(4) = the y velocity

Then in your derivative function you will have something like this:

dy=zeros(4,1); % 4th order system so 4 elements
dy(1) = y(3);
dy(2) = y(4);
% the -GM *rvec/r^3 calculation goes here (this will be a vector result)
dy(3) = _____; % you fill this in, based on the vector result above

dy(4) = _____; % you fill this in, based on the vector result above

Side Note: The DEQ in your Question is missing a minus sign (even though you do have it in your code)

Best Answer

Related Solutions

MATLAB: Colour of legend doesn’t change

MATLAB: Get x and y coordinates of motion from ode45

Related Question