Your comments

??? Error using ==> mp>mp.pinv at 3009

The algorithm failed to compute a singular value.
The update process of divide and conquer failed.
Please send us the input matrix to support@advanpix.com
Error in ==> LM_Xc at 149
s=pinv(A)*den;

matrix sent to pavel@advanpix.com

clear all

clc
addpath('C:\........\Multiprecision Computing Toolbox\')

digito=2000;

mp.Digits(digito);

load variables.mat

s=pinv(a)*den; %a\den;


??? Error using ==> mp>mp.pinv at 2964

DGESSD failed with exit code = 3
Error in ==> productmatrix at 11

s=pinv(a)*den; %a\den;

could create length mantissa and exponent precision real, arbitrary?

would political project create real length and precision, arbitrary ?

Using the code that I shared with you, was a way to find the orders of magnitude of R and Rtc

mp.Digits(N);

ExpMax=fix(log(mp.realmax)/log(10))

ExpMin=fix(log(mp.realmin)/log(10))

___________________________________

N=100

ExpMax = 323228495.0000

ExpMin =-323228496.0000


N=1000

ExpMax = 323228495.0000

ExpMin =-323228496.0000

function Num=Random1(Max,Min)

r=mp(rand(1,1));

Num=abs(Min+(Max-Min)*r);

---------------------------------------------------

% main program


clear all

clc

% Estimator order parameters

% I need to estimate the order of the parameters, because it was out range

mp.Digits(34)

GMinR=mp('0');

GMaxR=mp('1e+10000');

GMinRtc=mp('0');

GMaxRtc=mp('1e+10000');

GMinc=mp('0');

GMaxc=mp('1');

fre= mp('0.13');

DRe= mp('8241.3');
DIm= mp('-14145');

for j=1:30

j

for k=1:1000

R=Random1(GMinR,GMaxR);

Rtc=Random1(GMinRtc,GMaxRtc);

c=Random1(GMinc,GMaxc);


Fun=mp(funcion(fre,R,Rtc,c,DRe,DIm));

L(k,:)=[R Rtc c Fun];

end

L=sortrows(L,4);


GMaxR=L(1,1);

GMaxRtc=L(1,2);

GMaxc=L(1,3);

end


GMaxR

GMaxRtc

GMaxc


Fun=L(1,4)
----------------------------------------
OUT
GMaxR = NaN
GMaxRtc = NaN
GMaxc = 0.7807950309700372182319938474392984

Fun = NaN

function Num=Random1(Max,Min)

r=mp(rand(1,1));

Num=Min+(Max-Min)*r;

---------------------------------

main program


mp.Digits(1000000)

GMin=mp(0);

GMax=mp(1e+5000)

Number=abs(Random1(GMin,GMax))

---------------------------------

OUT

GMax=inf

Number = NaN


I thought you mean the next version of Multiprecision Computing Toolbox for MATLAB ADVANPIX

si sera en la próxima versión y en muy pocos días puedo esperar

if it will be in the next version and in a few days I can wait