harj53.m

% Lasketaan ympyräkartion tilavuus
function harj53()

    h = 10; % Kartion korkeus
    r0 = 2; % Kartion pohjan säde

    Vanalyyttinen = pi*r0^2*h/3;

    V = 0;
    dx = 0.1;
    dy = 0.1;
    dz = 0.1;
    dV = dx*dy*dz;

    for z = 0:dz:10
        for y = -2:dy:2
            for x = -2:dx:2
                r = sqrt(x^2 + y^2);
                if z < h-h*r/r0          % Ollaanko sisällä?
                    V = V + dV;
                end    
            end
        end
    end

    fprintf( 'Tilavuus = %g\n', V )
    fprintf( 'Analyyttinen = %g\n', Vanalyyttinen )
end