Plotting function in interval

Question

I have the following function :

ySol2=(2*(x^3 + 1)^(1/2))/cos(x) +2/cos(x)

my question is how can I draw it in the interval 1000;5000 with the constrains cos(x) != 0 and x> -1?

Answer 1

Something like this perhaps?

% your function. Note element-wise operations (.^, ./)
ySol2 = @(x) (2*(x.^3 + 1).^(1/2))./cos(x) +2./cos(x);
% your x interval, with step 0.1. Note also that the argument for 'cos' is in [rad].  
x = 1000:0.1:5000;
% calculate function values
y = ySol2(x);
% keep only non-constrained values
ind = (cos(x)~=0) | (x>-1);
x = x(ind);
y = y(ind);
% plot
figure;
plot(x,y);