Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Triple Integral for earth potensial

  1. Apr 22, 2013 #1
    I am trying to solve this triple integral for computing the potential: G*ρ*∫∫∫-atan((x*y)/(z*sqrt(x^2+y^2+z^2)))dxdydz with the limits x = 0 to 1000 z = 0 to 1000 and z = 0 to 1. The G is tha gravitational constant and ρ is the density of rock/earth.
    I have tried to use multiple softwares(matlab, maxima, wolframalpha) for solving this, but can't find any solutions. I am not sure of what I am doing wrong, can anyone help me?

  2. jcsd
  3. Apr 22, 2013 #2
    And i want the answar to be in mgal (milli gal)
  4. Apr 22, 2013 #3
    Your integrand choice seems a bit weird for the potential, are you sure this is the correct setup? Is this the potential of a slab of earth? or for the whole earth? (your bounds of integration are for the first case)

    No one will know what is wrong with your attempts to use software unless you show what you input into the program.
  5. Apr 22, 2013 #4


    User Avatar
    Science Advisor
    Homework Helper

    If x and y are lengths, "atan(x*y)" doesn't make any sense. You can't take the inverse tangent of an area.

    Maybe you meant atan(x/y)? Did you make the same mistake using Wolfram, etc?
  6. Apr 22, 2013 #5
    First, the formula for V_zz(potensial) has been given to me by my professor, the limits I use is for testing that I can compute the integral (but I can't...). I am supposed to solve a problem about how the masses around the gravimeter is working on it. I have measured the gradient in two different places, in the basement and one level up. I am supposed to get the same gradient and potensial when I move the gravimeter up from the basement (making a model). The formula is for computing the potensial of a prism.
    The ρ-value is supposed to be 2670.

    Okey, so this is how I try to do it in Maxima (just the inetgral, without the constants):
    But this will only compute the dx.

    In matlab I am trying to use the integral3 function, like this:
    Code (Text):
    f = @(x,y,z) (-atan((x.*y)./(z.*sqrt(x.^2+y.^2+z.^2))))

    xmin = 0
    xmax = 1000
    ymin = @ (x) 0
    ymax = @ (x) 1000
    zmax = @ (x,y) 1
    zmin = @ (x,y) 0

    Q = integral3(f,xmin,xmax,ymin,ymax,zmin,zmax)
    But this will have some dimensions porblems, so matlab wont execute.

    I am not sure if I am supposed to make a bougerplate or if this is some kind of bougerplate?
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook