Dear All,
I'm trying to substitute a series ansatz into a coupled couple of equations to find coefficients. The equations are quite cumbersome, but the principle shouldn't be too difficult. The equations are
f3A(alpha,mu,Omega,delta,m) =0
f4A(alpha,mu,Omega,m) = 0
where f3A and f4A are complicated functions of the several variables. The region I'm working in is
delta=0, m=1
Now, I have reason to believe that the following ansatz will work:
alpha = A0 + A2/Omega^2 + A4/Omega^4 +....
-mu=M0 + M2/Omega^2 + M4/Omega^4 +.....

Dear all,
I'm trying to solve a difficult metric in maple. I did the following:
restart;
with(tensor):
coords:=[t, r, theta,phi]:
g:=array(symmetric,sparse,1..4,1..4):
g[1,1]:=1+mu*a(r);
g[2,2]:=-(1+mu*b(r)):
g[3,3]:=-r^2*(1+mu*b(r)):
g[4,4]:=-r^2*(sin(theta))^2*(1+mu*b(r)):
g[4,1]:=mu*omega*f(r)*r^2*sin(theta)^2*(1+mu*b(r))/c:
metric:=create([-1,-1], eval(g));
tensorsGR(coords,metric,contra_metric,det_met, C1, C2, Rm, Rc, R, G, C):
display_allGR (coords,metric,contra_metric, det_met, C1, C2, Rm, Rc, R, G, C);

Dear All,
I'm trying to make an asymptotic expansion of the following integral:
A := m-> (x-sqrt(x^2+2*x*m+m^2-1))/(sqrt(x)*sqrt(x^2+2*x*m+m^2-1)) ;
IntegralA := m-> evalf(Int(A(x, m), x = 0 .. infinity)) ;
and similarly for two others
P3 := m-> -(-1+4*x^2*m^2+8*x^3*m+2*x*m-3*x^2+m^2)/((2*x^2+1+2*sqrt((x+m)^2/((x+m+1)*(x+m-1)))*x^2)*(x+m+1)*(x+m-1)*sqrt(x)) ;
IntegralP3:=m->evalf(Int(P3(x,m)/2,x=0..infinity));
and lastly
P1:=(x,m)->sqrt(x)/(sqrt((x+m)^2-1)*(x+m+sqrt((x+m)^2-1)));
IntegralP1:=m->evalf(Int(P1(x,1),x=0..infinity));

Hello,
Does anybody know if there's any way to display a 2-D plot & a 3-D plot on the same axes? I tried using display and I got the message:
Error, (in plots:-display) cannot display 2-D and 3-D plots together
Regards,
Drew

Hello,
This is related to my recent posts at
http://www.mapleprimes.com/forum/integral-equations#comment-8339
I have a function of three variables one of which is inside a (difficult) integral. I want to calculate & store the array of points which satisfy the function, then use transform to change those points and then plot them. The function looks like this:
P:=(x,m)->sqrt(x)/(sqrt((x+m)^2-1)*(x+m+sqrt((x+m)^2-1)));
IntegralP:=(m)->Int(P(x,m),x=0..infinity);
f4A:=(beta,Omega,m)->beta^2+(3*(beta*Omega)^(3/2)*evalf(IntegralP(m)))/(4*2^(3/4))-1;
In my earlier case, m=1 (and there was another variable, but it wasn't under the integral) so the integral was not a problem. I tried the same method;