Explicit integral could not be found

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))))

More Answers (0)

Categories

Find more on Physics 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!