I need to graph this function C(R) and I have problems with the vector R. The code is:

clear all

close all

clc

C0= 1500;

Dab= 1.3896e-7;

Mw= 152;

p= 1200000;

k= 36;

r=0.0025;

R= 0:0.1:10;

alpha= r*sqrt(k/Dab);

C=(C0*r/R)*exp(alpha*(1-R/r));

plot(R,C)