How did I get this crazy Plot?

clc
clear all
close all

A = 1.9;
R = 10;
a = 1;

X = linspace(0.0001,0.1,100);
% Y = a.*(1 + (A^(-R))*(X.^(R-1)));

for i = 1:length(X) 
    Y(i) = a*(1 + (A^(-R))*(X(i)^(R-1)));
end

plot(X, Y)

If you run the above code you will get a weird plot like this (shown below).

1. How come those many ones in the Y-axis?
2. Instead of a straight line, How come the plot it is a curve?

• Could anybody please explain??*

[Dennis]

Y is not really 1 at every point. It is a question of precision, try vpa(Y).

For i=1 Y is 1 + 1.63e-39

[Learner]

Yes, Solved. Thanks a lot...

[Guillaume]

This is the plot I get (R2018a).

I see nothing unusual about it and I don't understand why you expect a straight line when you're effectively plotting y = c1*x.^9 + c2.