error in solving numerical integral

Question

0 个投票

Hi,

I want t slve a Numerical Integral fOr Leonard-johnes Potential

How can I Solve this error ?

Fun=@(x,t)(x^2*(exp(-((1/(x^12))-(1/x^6))/t))-1));
Q=integral(@(x)Fun(x,10),0,2)

Error In integral

I Also need to plot Q in different t

Do You Have Any Comments ?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

Answer 1

madhan ravi 2019-1-17

编辑：madhan ravi 2019-1-17

0 个投票

syms t x % at the very beginning
%remove @(x,t) and replace it with matlabFunction()

Answer 2

Walter Roberson 2019-1-17

1 个投票

Fun=@(x,t) x.^2.*(exp(-((1./(x.^12))-(1./x.^6))./t))-1

You had too many ) and you needed to vectorize. integral() always invokes the function with a vector of values rather than just a single value.