Explicit integral could not be found
1 view (last 30 days)
Show older comments
Hello There, I have a problem on integrating a function. The function is important to describe band bending on the surface of a semiconductor.. The function already contains MTFA factor and potential dependance Fermi-Dirac function. I want to integrate w.r.t E by keeping V and z as constant from 0 to inf. When I use usual procedure I get an error Warning: Explicit integral could not be found. I was able to integrate for a small energy(E) interval by finding the area under the curve. But I need an efficient way to solve this problem. Is there any better way? Thanks for your help.
Here is the function..
Syms E V z;
f= -(160822346162939158374951306007842683360580644730317667175195738112000000000000000*E*((315656633149031*sin((75557863725914323419136*pi*z*((10000*E)/3543 + 1)^(1/2)*((5534527050034529*E)/140737488355328)^(1/2))/315656633149031))/(75557863725914323419136*pi*z*((10000*E)/3543 + 1)^(1/2)*((5534527050034529*E)/140737488355328)^(1/2)) - 1))/(618501041962717063047950574462370983128147*(exp((5534527050034529*E)/140737488355328 + (5534527050034529*V)/140737488355328 - 7742803342998306071/14073748835532800) + 1)*((98722420244898788167299958078358473932800*((6877444662696557*E)/618970019642690137449562112)^(1/2))/(100093691394021*((34081300774431287738506648835156271556395069066039782704939008000000*E)/1113194117431279288697172049 + (((34081300774431287738506648835156271556395069066039782704939008000000*E)/1113194117431279288697172049 - 15026719/25000000000)^2 - ((278101189509843229200380651449663578591555031138635517788160000000000*E)/3339582352293837866091516147 + 13522409/300000000)^3)^(1/2) - 15026719/25000000000)^(1/3)) - (((98722420244898788167299958078358473932800*((6877444662696557*E)/618970019642690137449562112)^(1/2)*((278101189509843229200380651449663578591555031138635517788160000000000*E)/3339582352293837866091516147 + 13522409/300000000)^2)/33364563798007 - (604921630050617324495130493125141549023232*((34081300774431287738506648835156271556395069066039782704939008000000*E)/1113194117431279288697172049 - 15026719/25000000000)*((6877444662696557*E)/618970019642690137449562112)^(1/2))/834114094950175)/(2*(((34081300774431287738506648835156271556395069066039782704939008000000*E)/1113194117431279288697172049 - 15026719/25000000000)^2 - ((278101189509843229200380651449663578591555031138635517788160000000000*E)/3339582352293837866091516147 + 13522409/300000000)^3)^(1/2)) - (302460815025308662247565246562570774511616*((6877444662696557*E)/618970019642690137449562112)^(1/2))/834114094950175)/(3*((34081300774431287738506648835156271556395069066039782704939008000000*E)/1113194117431279288697172049 + (((34081300774431287738506648835156271556395069066039782704939008000000*E)/1113194117431279288697172049 - 15026719/25000000000)^2 - ((278101189509843229200380651449663578591555031138635517788160000000000*E)/3339582352293837866091516147 + 13522409/300000000)^3)^(1/2) - 15026719/25000000000)^(2/3)) + (3952027511271527982106121913080217600000*((6877444662696557*E)/618970019642690137449562112)^(1/2))/33364563798007 + ((((98722420244898788167299958078358473932800*((6877444662696557*E)/618970019642690137449562112)^(1/2)*((278101189509843229200380651449663578591555031138635517788160000000000*E)/3339582352293837866091516147 + 13522409/300000000)^2)/33364563798007 - (604921630050617324495130493125141549023232*((34081300774431287738506648835156271556395069066039782704939008000000*E)/1113194117431279288697172049 - 15026719/25000000000)*((6877444662696557*E)/618970019642690137449562112)^(1/2))/834114094950175)/(2*(((34081300774431287738506648835156271556395069066039782704939008000000*E)/1113194117431279288697172049 - 15026719/25000000000)^2 - ((278101189509843229200380651449663578591555031138635517788160000000000*E)/3339582352293837866091516147 + 13522409/300000000)^3)^(1/2)) - (302460815025308662247565246562570774511616*((6877444662696557*E)/618970019642690137449562112)^(1/2))/834114094950175)*((278101189509843229200380651449663578591555031138635517788160000000000*E)/3339582352293837866091516147 + 13522409/300000000))/(3*((34081300774431287738506648835156271556395069066039782704939008000000*E)/1113194117431279288697172049 + (((34081300774431287738506648835156271556395069066039782704939008000000*E)/1113194117431279288697172049 - 15026719/25000000000)^2 - ((278101189509843229200380651449663578591555031138635517788160000000000*E)/3339582352293837866091516147 + 13522409/300000000)^3)^(1/2) - 15026719/25000000000)^(4/3))))
0 Comments
Accepted Answer
More Answers (0)
See Also
Categories
Find more on Calculus 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!