Matematické Fórum

Function ObjemE(p#,L#,R#,B#) as Double
' Tato funkce vypočítá objem válcové nádrže s eliptickými čely
' Hodnoty mimo povolený vstupní rozsah se berou jako krajní
' p - zadaná výška hladiny [rozsah 0 až 1], 1 je maximum - odpovídá 100%
' L - je délka válcové části bez eliptických čel
' R - je poloměr válcové části
' B - velikost vyboulení eliptického čela 
    Dim q#, Pi as Double
    Pi = 3.1415926535897932385
    If p<=0 Then ' Nulový objem nádrže [nula je 0%]
        ObjemE = 0
    ElseIf p>=1 Then 'Maximální Objem nádrže [jednička je 100%]
        ObjemE = Pi*R^2*(3*L+4*B)/3
    ElseIf p=1/2 Then ' Poloviční objem nádrže [1/2 je 50%]
        ObjemE = Pi*R^2*(3*L+4*B)/6
    Else ' Všechny ostatní případy budeme počítat ...
        q = p-p^2
        ObjemE = R^2*(L*(atn((p-1/2)/sqr(q))+(4*p-2)*sqr(q)+Pi/2)+4*Pi*(1-2*p/3)*p^2*B)
    End If
End Function

Function pInvQ(V#,L#,R#,B#) as Double
' Tato funkce počítá procentovou výšku z objemu, je inverzní vůči předchozí funkci
' Hodnoty mimo povolený vstupní rozsah se berou jako krajní
' V - zadaný objem nádrže 
' L - je délka válcové části bez eliptických čel
' R - je poloměr válcové části
' B - velikost vyboulení eliptického čela 
    Dim p#, z#, m#, o#, n#, h#, g#, q#, Pi as Double
    Pi = 3.1415926535897932385
    If V<=0 Then ' Nulová výška nádrže je při nulovém objemu
        p = 0
    Elseif V>=Pi*R^2*(3*L+4*B)/3 Then ' Maximální výška nádrže je při maximálním objemu
        p = 1
    Else ' Všechny ostatní případy se budou počítat ...
        p = V/(Pi*R^2*(3*L+4*B)/3)
        z = 1/2
        While abs(p-z)>1e-12 
            z = p
            m = 2*p-1 
            n = sqr(1-m^2)
            q = p-p^2
            o = atn(-m/(2*sqr(q)))
            g = 6*sqr(q)
            h = 6*(o-n-Pi/2)
            p = (p*(6*(p-1)*V+((3*n*(4*p^2-3)+h*(p-1)-g)*L-8*Pi*(p-1)*p^2*(4*p-3)*B)*R^2)) _
                /(((3*(2*m^2-1)*n-g)*L-3*Pi*(m-1)^2*(m+1)^2*B)*R^2)
        Wend
    End If
    pInvQ = p
End Function