Simplified 3D Moffat

(%i1) f3s3(x3,x2,x1,alpha,beta) := ( 1 + (x1^2 + x2^2 +x3^2)/ alpha ^2 )^(-beta);

Result

(%i2) f3s2(x2,x1,alpha,beta) := defint( f3s3(x3,x2,x1,alpha,beta), x3, -inf, inf );

Result

(%i3) f3s1(x1,alpha,beta) := defint( f3s2(x2,x1,alpha,beta), x2, -inf, inf );

Result

(%i4) f3s0(alpha,beta) := defint( f3s1(x1,alpha,beta), x1, -inf, inf );

Result

(%i5) norm(alpha, beta) := factor(makegamma(ratsimp(f3s0(alpha,beta))));

Result

(%i29) killcontext(initial); assume( alpha > 0, beta > 0, beta < 0.5 ); norm(alpha,beta);

Result

(%i32) killcontext(initial); assume( alpha > 0, equal(beta,0.5)); norm(alpha,beta);

Result

(%i35) killcontext(initial); assume( alpha > 0, beta > 0.5, beta < 1 ); norm(alpha,beta);

Result

(%i40) killcontext(initial); assume( alpha > 0, equal(beta,1)); norm(alpha,beta);

Result

(%i43) killcontext(initial); assume( alpha > 0, beta > 1, beta < 3/2 ); norm(alpha,beta);

Result

(%i46) killcontext(initial); assume( alpha > 0, equal(beta,3/2) ); norm(alpha,beta);

Result

(%i49) killcontext(initial); assume( alpha > 0, beta > 3/2, beta < 3 ); norm(alpha,beta);

Result

(%i70) killcontext(initial); assume( alpha > 0, equal(beta,3)); norm(alpha,beta);

Result

(%i55) killcontext(initial); assume( alpha > 0, beta > 3, beta < 7/2 ); norm(alpha,beta);

Result

(%i64) killcontext(initial); assume( alpha > 0, equal(beta,7/2) ); norm(alpha,beta);

Result

(%i67) killcontext(initial); assume( alpha > 0, beta > 7/2 );norm(alpha,beta);

Result


Created with wxMaxima.