MATLAB: How to solve a system of equations involving distribution

equationssolve

 gamcdf(3.01,a,b)=0.4
 gamcdf(4.63,a,b)=0.6

how to solve for a and b?

I tried solver & evalin(symengine)

Best Answer

this worked:

 F = [gamcdf(6.63,x(1),x(2))-0.4;
       gamcdf(8.89,x(1),x(2))-0.6];
 x0 = [2; 0.8];            
 options=optimset('Display','iter');   
 [x,fval] = fsolve(@SPI05,x0,options)

spasiba Andrei !

Related Solutions

MATLAB: Error in using fsolve

try this cod:

b = [4;6];
F = @(x,b)[gamcdf(b(1),x(1),x(2))-0.4;gamcdf(b(2),x(1),x(2))-0.6];
[x,fval] = fsolve(@(x)F(x,b),[2;.8])

MATLAB: How to solve singularity problem while using fsolve

Try splitting the fsolve approach with real() and imag() because with real() i get something:

..
F = @(x)[real(((x-a)/(x-c))^(m))+((x-a)/(x-b))];
..
[xd, fval, exitflag, output]= fsolve(F, x0, options)
                                           Norm of      First-order   Trust-region
   Iteration  Func-count     f(x)          step         optimality    radius
       0          2        0.366931                          1.77               1
       1          4     0.000271336       0.207172         0.0421               1
       2          6     2.03148e-07     0.00645204        0.00109               1
       3          8     1.38596e-13    0.000186725       8.97e-07               1
       4         10     7.81816e-26    1.54486e-07       6.74e-13               1
       5         12     3.08149e-33    1.16029e-13       1.34e-16               1
Equation solved, inaccuracy possible.
The vector of function values is near zero, as measured by the selected value
of the function tolerance. However, the last step was ineffective.
<stopping criteria details>
xd =
          0.30
fval =
         -0.00
exitflag =
          3.00
output = 
       iterations: 5.00
        funcCount: 12.00
        algorithm: 'trust-region-dogleg'
    firstorderopt: 0.00
          message: 'Equation solved, inaccuracy possible.…'

repeat for the imag() part.

Or alternatively use abs() and arg()

all together, you should catch at least 2 poles on b and c.

If you find this answer of any help solving your question,

please click on the thumbs-up vote link,

thanks in advance

John

jgb2012@sky.com

Best Answer

Related Solutions

MATLAB: Error in using fsolve

MATLAB: How to solve singularity problem while using fsolve

Related Question