double integration of torsional function

Hello,

I am new to matlab, I want to calculate this double integral of torsional equations. I am struggling to find an answer for this

Answer 1

[David Hill]

syms x h y;
eqn=x^3-h*x^2-3*x*y^2+4*h^3/27-h*y^2;
I=4*int(int(eqn,y,0,(2*sqrt(3)*h-3*sqrt(3)*x)/9),x,0,h/3);