Hello, I am try to solve the simultaneous equations with integrand recently, but I am stuck now, below is the detail:
The simultaneous equations that I want to solve:
and I write some code like this:
kursi = 1; b1 = 3*1.33+1; b2 = 3*1+1; alpha = 0.05; beta = 0.10;
fun1 = @(t,c0,n) chi2cdf( (n-1)*(b1*sqrt(n)-t).^2./(9*n*c0.^2),n-1 )...
*( normpdf(t+kursi*sqrt(n))+normpdf(t-kursi*sqrt(n)) );
fun2 = @(t,c0,n) chi2cdf( (n-1)*(b2*sqrt(n)-t).^2./(9*n*c0.^2),n-1 )...
*( normpdf(t+kursi*sqrt(n))+normpdf(t-kursi*sqrt(n)) );
eq = @(c0,n) [ integral( @(t)fun1(t,c0,n),0,b1*sqrt(n) )-(1-alpha);...
integral( @(t)fun2(t,c0,n),0,b2*sqrt(n) )-beta ];
result = fsolve(eq,[0,100]);
the feedback from Matlab:
Not enough input arguments.
Error in
@(c0,n)[integral(@(t)fun1(t,c0,n),0,b1*sqrt(n))-(1-alpha);integral(@(t)fun2(t,c0,n),0,b2*sqrt(n))-beta]
Error in fsolve (line 230)
fuser = feval(funfcn{3},x,varargin{:});
Caused by:
Failure in initial objective function evaluation. FSOLVE cannot continue.
===================================================================== After comparing to the doc fsolve, I guess the problem are function handle and integral. Does anyone can give me some advices? Thanks!
