%%%%%%%%%
% pot.m %
%%%%%%%%%

clf;

% x = [0.51603557035733 4.08888282515727 0.97703114535312]; %
% x = [0.92231006075817 2.21570837209438 1.94740518295180]; %
x = [0.37742291628244 3.57823943232922 0.80848585672605];
B = 2 * x(3);
A = x(2) + log(x(1)) / B;
D = x(3);
C = x(2) + log(2 * x(1)) / D;
[A B C D]

A = 3.62229716811456;
B = 1.56475000320749;
C = 0.39558139377561;
D = 0.70741137104679;
E = 0.5;
RC = 6.10699597371401;

r = 1:0.01:RC;

y = A * ( ( B ./ r) .^ C - D ) .* exp(E ./ (r - RC));
plot(r,y);

axis ([1 RC -0.5 1]);

