Plotting problem of function in one variable and four known parameters

Hello I have a problem with the following code. From previous calculations I have obtained the equation of T (don't be alarmed by the length of the equation, it comes from previous codes which I will not explain) which expressed as a function of the variable "radius" and the known parameters B, v, p, r, R.

After assigning numerical values to the known parameters and values ranging from 0.25 to 0.50 of the variable "radius", I used the following code to plot the equation:

B = 20230000000;
p = 6000000;
R = 0.50;
r = 0.25;
v = 0.3;
raggio = linspace(0.25,0.50,26);
T = arrayfun(@(raggio)B*(((R^6*p - 3*p*r^6 - p*r^6*v - 2*R^6*p*log(R) - R^2*p*r^4 + 3*R^4*p*r^2 + 10*p*r^6*log(r) - R^6*p*v - 2*R^2*p*r^4*log(R) - 12*R^4*p*r^2*log(R) - 12*R^2*p*r^4*log(r) + 18*R^4*p*r^2*log(r) + 24*R^2*p*r^4*log(R)^2 + 8*R^4*p*r^2*log(r)^2 + 2*R^6*p*v*log(R) + R^2*p*r^4*v + R^4*p*r^2*v + 2*p*r^6*v*log(r) + 8*R^4*p*r^2*v*log(r)^2 + 2*R^2*p*r^4*v*log(R) - 4*R^4*p*r^2*v*log(R) - 4*R^2*p*r^4*v*log(r) + 2*R^4*p*r^2*v*log(r) + 8*R^2*p*r^4*v*log(R)^2 - 24*R^2*p*r^4*log(R)*log(r) - 8*R^4*p*r^2*log(R)*log(r) - 8*R^2*p*r^4*v*log(R)*log(r) - 8*R^4*p*r^2*v*log(R)*log(r))/(32*B*(R^4 - r^4*v - R^4*v - 3*r^4 + 2*R^2*r^2 + 4*R^2*r^2*log(R)^2 + 4*R^2*r^2*log(r)^2 + 2*R^2*r^2*v - 8*R^2*r^2*log(R) + 8*R^2*r^2*log(r) + 4*R^2*r^2*v*log(R)^2 - 8*R^2*r^2*log(R)*log(r) + 4*R^2*r^2*v*log(r)^2 - 8*R^2*r^2*v*log(R)*log(r))) - (3*p*raggio^2)/(16*B) + (3*R^6*p - 15*p*r^6 - 3*p*r^6*v + 21*R^2*p*r^4 - 9*R^4*p*r^2 - 3*R^6*p*v - 36*R^2*p*r^4*log(R) + 12*R^4*p*r^2*log(R) + 36*R^2*p*r^4*log(r) - 12*R^4*p*r^2*log(r) + 3*R^2*p*r^4*v + 3*R^4*p*r^2*v - 12*R^2*p*r^4*v*log(R) + 12*R^4*p*r^2*v*log(R) + 12*R^2*p*r^4*v*log(r) - 12*R^4*p*r^2*v*log(r))/(32*B*(R^4 - r^4*v - R^4*v - 3*r^4 + 2*R^2*r^2 + 4*R^2*r^2*log(R)^2 + 4*R^2*r^2*log(r)^2 + 2*R^2*r^2*v - 8*R^2*r^2*log(R) + 8*R^2*r^2*log(r) + 4*R^2*r^2*v*log(R)^2 - 8*R^2*r^2*log(R)*log(r) + 4*R^2*r^2*v*log(r)^2 - 8*R^2*r^2*v*log(R)*log(r))) + (log(raggio)*(R^6*p - 5*p*r^6 - p*r^6*v + 7*R^2*p*r^4 - 3*R^4*p*r^2 - R^6*p*v - 12*R^2*p*r^4*log(R) + 4*R^4*p*r^2*log(R) + 12*R^2*p*r^4*log(r) - 4*R^4*p*r^2*log(r) + R^2*p*r^4*v + R^4*p*r^2*v - 4*R^2*p*r^4*v*log(R) + 4*R^4*p*r^2*v*log(R) + 4*R^2*p*r^4*v*log(r) - 4*R^4*p*r^2*v*log(r)))/(16*B*(R^4 - r^4*v - R^4*v - 3*r^4 + 2*R^2*r^2 + 4*R^2*r^2*log(R)^2 + 4*R^2*r^2*log(r)^2 + 2*R^2*r^2*v - 8*R^2*r^2*log(R) + 8*R^2*r^2*log(r) + 4*R^2*r^2*v*log(R)^2 - 8*R^2*r^2*log(R)*log(r) + 4*R^2*r^2*v*log(r)^2 - 8*R^2*r^2*v*log(R)*log(r))) - (R^2*r^2*(2*R^4*p + 4*p*r^4 + p*r^4*v - R^4*p*log(R) - 6*R^2*p*r^2 + 5*p*r^4*log(R) + R^4*p*log(r) - 5*p*r^4*log(r) + R^4*p*v - R^4*p*v*log(R) - 2*R^2*p*r^2*v + p*r^4*v*log(R) + R^4*p*v*log(r) - p*r^4*v*log(r)))/(16*B*raggio^2*(R^4 - r^4*v - R^4*v - 3*r^4 + 2*R^2*r^2 + 4*R^2*r^2*log(R)^2 + 4*R^2*r^2*log(r)^2 + 2*R^2*r^2*v - 8*R^2*r^2*log(R) + 8*R^2*r^2*log(r) + 4*R^2*r^2*v*log(R)^2 - 8*R^2*r^2*log(R)*log(r) + 4*R^2*r^2*v*log(r)^2 - 8*R^2*r^2*v*log(R)*log(r))))/raggio - ((raggio*(R^6*p - 3*p*r^6 - p*r^6*v - 2*R^6*p*log(R) - R^2*p*r^4 + 3*R^4*p*r^2 + 10*p*r^6*log(r) - R^6*p*v - 2*R^2*p*r^4*log(R) - 12*R^4*p*r^2*log(R) - 12*R^2*p*r^4*log(r) + 18*R^4*p*r^2*log(r) + 24*R^2*p*r^4*log(R)^2 + 8*R^4*p*r^2*log(r)^2 + 2*R^6*p*v*log(R) + R^2*p*r^4*v + R^4*p*r^2*v + 2*p*r^6*v*log(r) + 8*R^4*p*r^2*v*log(r)^2 + 2*R^2*p*r^4*v*log(R) - 4*R^4*p*r^2*v*log(R) - 4*R^2*p*r^4*v*log(r) + 2*R^4*p*r^2*v*log(r) + 8*R^2*p*r^4*v*log(R)^2 - 24*R^2*p*r^4*log(R)*log(r) - 8*R^4*p*r^2*log(R)*log(r) - 8*R^2*p*r^4*v*log(R)*log(r) - 8*R^4*p*r^2*v*log(R)*log(r)))/(32*B*(R^4 - r^4*v - R^4*v - 3*r^4 + 2*R^2*r^2 + 4*R^2*r^2*log(R)^2 + 4*R^2*r^2*log(r)^2 + 2*R^2*r^2*v - 8*R^2*r^2*log(R) + 8*R^2*r^2*log(r) + 4*R^2*r^2*v*log(R)^2 - 8*R^2*r^2*log(R)*log(r) + 4*R^2*r^2*v*log(r)^2 - 8*R^2*r^2*v*log(R)*log(r))) - (p*raggio^3)/(16*B) + (raggio*(R^6*p - 5*p*r^6 - p*r^6*v + 7*R^2*p*r^4 - 3*R^4*p*r^2 - R^6*p*v - 12*R^2*p*r^4*log(R) + 4*R^4*p*r^2*log(R) + 12*R^2*p*r^4*log(r) - 4*R^4*p*r^2*log(r) + R^2*p*r^4*v + R^4*p*r^2*v - 4*R^2*p*r^4*v*log(R) + 4*R^4*p*r^2*v*log(R) + 4*R^2*p*r^4*v*log(r) - 4*R^4*p*r^2*v*log(r)))/(32*B*(R^4 - r^4*v - R^4*v - 3*r^4 + 2*R^2*r^2 + 4*R^2*r^2*log(R)^2 + 4*R^2*r^2*log(r)^2 + 2*R^2*r^2*v - 8*R^2*r^2*log(R) + 8*R^2*r^2*log(r) + 4*R^2*r^2*v*log(R)^2 - 8*R^2*r^2*log(R)*log(r) + 4*R^2*r^2*v*log(r)^2 - 8*R^2*r^2*v*log(R)*log(r))) + (raggio*log(raggio)*(R^6*p - 5*p*r^6 - p*r^6*v + 7*R^2*p*r^4 - 3*R^4*p*r^2 - R^6*p*v - 12*R^2*p*r^4*log(R) + 4*R^4*p*r^2*log(R) + 12*R^2*p*r^4*log(r) - 4*R^4*p*r^2*log(r) + R^2*p*r^4*v + R^4*p*r^2*v - 4*R^2*p*r^4*v*log(R) + 4*R^4*p*r^2*v*log(R) + 4*R^2*p*r^4*v*log(r) - 4*R^4*p*r^2*v*log(r)))/(16*B*(R^4 - r^4*v - R^4*v - 3*r^4 + 2*R^2*r^2 + 4*R^2*r^2*log(R)^2 + 4*R^2*r^2*log(r)^2 + 2*R^2*r^2*v - 8*R^2*r^2*log(R) + 8*R^2*r^2*log(r) + 4*R^2*r^2*v*log(R)^2 - 8*R^2*r^2*log(R)*log(r) + 4*R^2*r^2*v*log(r)^2 - 8*R^2*r^2*v*log(R)*log(r))) + (R^2*r^2*(2*R^4*p + 4*p*r^4 + p*r^4*v - R^4*p*log(R) - 6*R^2*p*r^2 + 5*p*r^4*log(R) + R^4*p*log(r) - 5*p*r^4*log(r) + R^4*p*v - R^4*p*v*log(R) - 2*R^2*p*r^2*v + p*r^4*v*log(R) + R^4*p*v*log(r) - p*r^4*v*log(r)))/(16*B*raggio*(R^4 - r^4*v - R^4*v - 3*r^4 + 2*R^2*r^2 + 4*R^2*r^2*log(R)^2 + 4*R^2*r^2*log(r)^2 + 2*R^2*r^2*v - 8*R^2*r^2*log(R) + 8*R^2*r^2*log(r) + 4*R^2*r^2*v*log(R)^2 - 8*R^2*r^2*log(R)*log(r) + 4*R^2*r^2*v*log(r)^2 - 8*R^2*r^2*v*log(R)*log(r))))/raggio^2 - (3*p*raggio)/(8*B) + (R^6*p - 5*p*r^6 - p*r^6*v + 7*R^2*p*r^4 - 3*R^4*p*r^2 - R^6*p*v - 12*R^2*p*r^4*log(R) + 4*R^4*p*r^2*log(R) + 12*R^2*p*r^4*log(r) - 4*R^4*p*r^2*log(r) + R^2*p*r^4*v + R^4*p*r^2*v - 4*R^2*p*r^4*v*log(R) + 4*R^4*p*r^2*v*log(R) + 4*R^2*p*r^4*v*log(r) - 4*R^4*p*r^2*v*log(r))/(16*B*raggio*(R^4 - r^4*v - R^4*v - 3*r^4 + 2*R^2*r^2 + 4*R^2*r^2*log(R)^2 + 4*R^2*r^2*log(r)^2 + 2*R^2*r^2*v - 8*R^2*r^2*log(R) + 8*R^2*r^2*log(r) + 4*R^2*r^2*v*log(R)^2 - 8*R^2*r^2*log(R)*log(r) + 4*R^2*r^2*v*log(r)^2 - 8*R^2*r^2*v*log(R)*log(r))) + (R^2*r^2*(2*R^4*p + 4*p*r^4 + p*r^4*v - R^4*p*log(R) - 6*R^2*p*r^2 + 5*p*r^4*log(R) + R^4*p*log(r) - 5*p*r^4*log(r) + R^4*p*v - R^4*p*v*log(R) - 2*R^2*p*r^2*v + p*r^4*v*log(R) + R^4*p*v*log(r) - p*r^4*v*log(r)))/(8*B*raggio^3*(R^4 - r^4*v - R^4*v - 3*r^4 + 2*R^2*r^2 + 4*R^2*r^2*log(R)^2 + 4*R^2*r^2*log(r)^2 + 2*R^2*r^2*v - 8*R^2*r^2*log(R) + 8*R^2*r^2*log(r) + 4*R^2*r^2*v*log(R)^2 - 8*R^2*r^2*log(R)*log(r) + 4*R^2*r^2*v*log(r)^2 - 8*R^2*r^2*v*log(R)*log(r)))),raggio);
figure();
plot(raggio,T)

the error matlab exposes to me is as follows:

Error using ^

Incorrect dimensions for raising a matrix to a power. Check that the matrix is square and the power is a

scalar. To operate on each element of the matrix individually, use POWER (.^) for elementwise power.

Error in ProvaGrafico (line 7)

Please help me. Thank you.

[Torsten]

Corrected (see above).

Answer 1

[John D'Errico]

Does it not tell you the problem?

x = 1:5;

For example, what do these two operations do in MATLAB for the vector x?

x.^2

ans = 1×5

1 4 9 16 25

x^2

Error using ^
Incorrect dimensions for raising a matrix to a power. Check that the matrix is square and the power is a scalar. To operate on each element of the matrix individually, use POWER (.^) for elementwise
power.

As you can see, the first one works, allowing you to raise each element of the vector to the indicated power. The second fails, because the VECTOR x cannot be raised to a power. What is the square of a vector? It would be x^2 = x*x, correct? But that would fail, since linear algebra does not allow you to multiply a 1x5 array by a 1x5 array. The dimensions do not conform for multiplication.

Remember that * and / and ^ are all operators that apply to matrices and vectors. They just happen to work nicely for scalars too, but that gets us feeling lazy, so we try to use them for all operations. When you want to work on the elements of matrices and vectors, you need to use the .* and ./ and .^ operators.