i need to solve this hideous equation.
2*x3+(2^(1/2)*exp(-x3^2/2) + 2^(1/2))/(pi^(1/2)*erf((2^(1/2)*x3)/2))+(2^(1/2)*exp(-x3^2/2))/(pi^(1/2)*(erf((2^(1/2)*x3)/2) - 1))=0
my_function = @(x3) 2*x3+(2^(1/2)*exp(-x3^2/2) + 2^(1/2))/(pi^(1/2)*erf((2^(1/2)*x3)/2))+(2^(1/2)*exp(-x3^2/2))/(pi^(1/2)*(erf((2^(1/2)*x3)/2) - 1))==0;
z = solve(my_function)
when i compile this, it says "Empty sym: 0-by-1" i don't think the answer can be found. I tried the fpsolve, fzero functions too but none were effective. How do i find the root for that equation? Thank you

 채택된 답변

dpb
dpb 2018년 6월 2일

0 개 추천

Always try to plot() something to see what you have...
>> fnf = @(x3) 2*x3+(2^(1/2)*exp(-x3^2/2) + 2^(1/2))/(pi^(1/2)*erf((2^(1/2)*x3)/2))+(2^(1/2)*exp(-x3^2/2))/(pi^(1/2)*(erf((2^(1/2)*x3)/2) - 1));
>> ezplot(fnf)
Warning: Function failed to evaluate on array inputs; vectorizing the function may speed up its evaluation and avoid
the need to loop over array elements.
> In ezplot>ezplot1 (line 498)
In ezplot (line 154)
>> ylim([-8 8])
>>
yields
which clearly shows there is no intersection of the zero line

댓글 수: 2

Walter Roberson
Walter Roberson 2018년 6월 2일
... shows there is no real-valued solution.
dpb
dpb 2018년 6월 2일
Yes...

댓글을 달려면 로그인하십시오.

추가 답변 (1개)

Walter Roberson
Walter Roberson 2018년 6월 2일

0 개 추천

That equation has no real-valued solutions. It goes to negative infinity as x3 approaches 0 from below, and to positive infinity as x3 approaches 0 from above.
You might find it easier to rewrite it as a system of two equations, one involving the real part and the other involving the imaginary part, after having rewritten x3 to x3r+1i*x3i where x3r and x3i are both real-valued.

댓글 수: 1

Walter Roberson
Walter Roberson 2018년 6월 2일
My experimentation suggests there are two complex-valued roots, at approximately .3382099048573570581422493 +/- 1.359657273974301366628947*1i

댓글을 달려면 로그인하십시오.

카테고리

도움말 센터File Exchange에서 Programming에 대해 자세히 알아보기

제품

질문:

2018년 6월 2일

댓글:

dpb
2018년 6월 2일

Community Treasure Hunt

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

Start Hunting!

Translated by