How to make color magnitude by quiver plot?

Hello

I have a problem with quiver plot. My code makes the plot perfectly but I would like to have also color magnitude on it. I mean the electric field can be observed by color change when arrow rotates. I appreiate all the helps in this regard. My code and sample date is below:

clc
clear all

Ex=csvread('position2_Ex_new.csv');
Ey=csvread('position2_Ey_new.csv');
theta = Ex(2:end, 1);
phi = Ex(1, 2:end);
data_x = Ex(2:end,2:end);
data_x = transpose(data_x);
% data_x = max(data,0);
data_y = Ey(2:end,2:end);
data_y = transpose(data_y);
% data_y = max(data,0);

[theta_,phi_]= meshgrid(theta,phi);
r=1e-3;
x=r*sin(theta_).*cos(phi_);
y=r*sin(theta_).*sin(phi_);
z=r*cos(phi_);
quiver(x,y,data_x,data_y);
% xlabel('x');
% ylabel('y');
set(gca,'XDir','reverse');
set(gca,'DataAspectRatio',[10 10 10])
camroll(-90)

Answer 1

[KSSV]

You can use some fileexchange functions:

https://in.mathworks.com/matlabcentral/fileexchange/3225-quiverc

https://in.mathworks.com/matlabcentral/fileexchange/31480-quiverc-m-updated-with-new-features

https://in.mathworks.com/matlabcentral/fileexchange/7180-quiverc2wcmap

[Masoud Taleb]

I tried all of the above functions. They do not work! The show different errors...

[KSSV]

What errors you getting?

[Masoud Taleb]

First link stays in play mode and does not give any plot at all (I stayed 20 min without getting response)

Second link shows following error. I tried to fix it but was not successful.

Unrecognized function or variable 'samesize'.

Error in quiverc (line 129)
if ~isempty(samesize)

Error in Untitled (line 20)
quiverc(x,y,data_x,data_y);

Third link brings a ghraph which is totally something else.