I want to use integral&fzero to solve the equation but I could't debug by myself

2 views (last 30 days)
Hello, I have written some code below but I could't debug by myself. Does anyone can give me some advice? Thanks!
n = 100; Cpk_hat = 1.5; r = 0.95; zzz = 1.00;
fun = @(t,C)chi2cdf(((n-1)*((3*C+abs(zzz)*sqrt(n))-t).^2/(9*n*Cpk_hat.^2)),n-1)*...
(normpdf(t+zzz*sqrt(n))+normpdf(t-zzz*sqrt(n)));
C_guess = 1;
fzero(@(C) integral(@(t) fun(t,C),0,((3*C+abs(zzz))*sqrt(n)))-(1-r),C_guess);
error message:
Error using fzero (line 306)
FZERO cannot continue because user-supplied function_handle ==>
@(C)integral(@(t)fun(t,C),0,((3*C+abs(strange_hat))*sqrt(n)))-(1-r) failed with the error
below.
Inner matrix dimensions must agree.
  2 Comments
Adam
Adam on 25 Aug 2017
Simplest debugging option would be to move your complicated function into its own file, split into variables across multiple lines rather than all on one line, then just create a simple function handle to that to pass in. Then you will be able to debug it with the 'stop on errors' option from the breakpoints menu.
At a glance all your inputs are scalar so I can't see simply where your error would be.

Sign in to comment.

Accepted Answer

Matt J
Matt J on 25 Aug 2017
Define fun() with element-wise operators
fun=@(t,C)chi2cdf(((n-1).*((3.*C+abs(zzz).*sqrt(n))-t).^2./(9.*n.*Cpk_hat.^2)),n-1).*(normpdf(t+zzz.*sqrt(n))+normpdf(t-zzz.*sqrt(n)))
  4 Comments
Louis Liu
Louis Liu on 26 Aug 2017
@Matt J Hello, I have found the mistake and get the right value 1.3147 now. Thank for your help!

Sign in to comment.

More Answers (0)

Categories

Find more on 数值积分和微分 in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!