Minimize a Hexp function with respect to beta?

Questions
,
1

Hello Colleagues,

Am trying to solve for Hexp in the attached notebook in order to minimize this function with respect to beta but it seems mathematica is having trouble while evaluating some of the integrals in this function. Do you have any suggestion on how I can guide it evaluate these integrals?

Thanks

Fred

2

This is a bit confused because ro is used both as a parameter and a variable of integration. So first thing is to change the latter to a unique dummy variable. Then you can define

Hexp[\[Alpha]_?NumberQ, \[Beta]_?NumberQ, \[Rho]_?NumberQ, zi_?NumberQ] :=...

where … is similar to what you have, but with Integrate replaced by NIntegrate. Also the nsqrd should be changed into a function of similar form

nsqrd[\[Alpha]_?NumberQ, \[Beta]_?NumberQ, \[Rho]_?NumberQ, zi_?NumberQ] :=...
3

Hi Daniel,

Thanks for this suggestion. I have implemented the suggested changes but still mathematica can’t solve it !