How do I plot fractional orders of the Bessel function of the first kind (e.g. J_{1/3}(x))? Can I still use the besselj function?

clc;clear all;clf;

x=-20:0.1:20;

J1/3=besselj(1/3,x);

figure(1)

plot(x,J1/3,'k');

grid on

legend('J_{1/3}')

xlabel('x');ylabel('J_1/3(x)')

Answer 1

[VBBV]

clc;clear all;clf;
x=-20:0.1:20;
J13=besselj(1/3,x);  % change the variable befoe use in besselj function 
figure(1)
plot(x,real(J13),'k');
grid on
legend('J_{1/3}')
xlabel('x');ylabel('J_1/3(x)')

Use the correct variable naming convention