Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

Manipulating and posing hypothesis on functions

With the nice help of this community, I have been able to write this code replicating what is in the main books dedicated to theoretical insurance (I am sorry for the lengh of the code but it is necessary).

First I need this function

def use_prime(expr):
    r"""
    Return the expression with one-variable function derivatives "primed".
    """
    op = expr.operator()
    if op:
        args = expr.operands()
        aargs = (use_prime(a) for a in args)
        opp = op
        if (isinstance(op, sage.symbolic.operators.FDerivativeOperator)
                and len(args) == 1):
            name = op.function().name()
            primes = "'" * len(op.parameter_set())
            opp = function(f"{name}{primes}")
        return opp(*aargs)
    else:
        return expr

then the code

var('x y p w_0 θ I')
assume(p>=0)
assume(p<=1)
assume(θ>=1)
U = function('U')#déclare U comme le nom d'une fonction
EU(x,y,p) = p * U(x) + (1-p) * U(y)
w_1 = function('w_1')
w_1(w_0, π) = (w_0-π)
w_2 = function('w_2')
w_2(w_0, π, x, I, p) = (w_0 - π) - x + I 
π = function('π')
π(p,I,θ)=θ*p*I
U = function('U')
EU(x, y, p) = p*U(x) + (1 - p)*U(y)
x, y, p, w_0, θ, I, = SR.var('x y p w_0 θ I')
f = EU(w_1(w_0, π(p,I,θ)),w_2(w_0, π, x, I, π(p,I,θ)) )
h = diff(f, I)
hh = use_prime(h)
k = diff(f, I,2)
kk = use_prime(k)
show(LatexExpr(r'''\text{La condition d'optimalité du premier ordre est } dEU = 0 \text{, soit : } ''' ))
show(LatexExpr(r'\,\,\,\,\,\,\,\,\,\,\,\,\,\, dEU = '),hh, ' = 0')
d2EU=diff(EU(w_1(w_0, π),w_2(w_0, π, x, I, p),p),I,2)
show(LatexExpr(r'''\text{La condition d'optimalité du second ordre est } dEU = 0 \text{, soit : } ''' ))
show(LatexExpr(r'\,\,\,\,\,\,\,\,\,\,\,\,\,\,'),kk, ' < 0')

This should gives this

image description

And now come my questions :

1) Can I assume $U^\prime>0$,$U^{\prime\prime}<0$ and verify with the other asumptions that the second order is verified, without hand manipulation. 2) How to rewrite the first order condition in such a way that all the terms in $U^\prime$ be on right side of the equation (a ratio) and the other terms on the left side ? 3) But, I think to know the answer, how to substitute $x$ to $I$ in the final stage ?

Manipulating and posing hypothesis on functions

With the nice help of this community, I have been able to write this code replicating what is in the main books dedicated to theoretical insurance (I am sorry for the lengh of the code but it is necessary).

First I need this function

def use_prime(expr):
    r"""
    Return the expression with one-variable function derivatives "primed".
    """
    op = expr.operator()
    if op:
        args = expr.operands()
        aargs = (use_prime(a) for a in args)
        opp = op
        if (isinstance(op, sage.symbolic.operators.FDerivativeOperator)
                and len(args) == 1):
            name = op.function().name()
            primes = "'" * len(op.parameter_set())
            opp = function(f"{name}{primes}")
        return opp(*aargs)
    else:
        return expr

then the code

var('x y p w_0 θ I')
assume(p>=0)
assume(p<=1)
assume(θ>=1)
U = function('U')#déclare U comme le nom d'une fonction
EU(x,y,p) = p * U(x) + (1-p) * U(y)
w_1 = function('w_1')
w_1(w_0, π) = (w_0-π)
w_2 = function('w_2')
w_2(w_0, π, x, I, p) = (w_0 - π) - x + I 
π = function('π')
π(p,I,θ)=θ*p*I
U = function('U')
EU(x, y, p) = p*U(x) + (1 - p)*U(y)
x, y, p, w_0, θ, I, = SR.var('x y p w_0 θ I')
f = EU(w_1(w_0, π(p,I,θ)),w_2(w_0, π, x, I, π(p,I,θ)) )
h = diff(f, I)
hh = use_prime(h)
k = diff(f, I,2)
kk = use_prime(k)
show(LatexExpr(r'''\text{La condition d'optimalité du premier ordre est } dEU = 0 \text{, soit : } ''' ))
show(LatexExpr(r'\,\,\,\,\,\,\,\,\,\,\,\,\,\, dEU = '),hh, ' = 0')
d2EU=diff(EU(w_1(w_0, π),w_2(w_0, π, x, I, p),p),I,2)
show(LatexExpr(r'''\text{La condition d'optimalité du second ordre est } dEU = 0 \text{, soit : } ''' ))
show(LatexExpr(r'\,\,\,\,\,\,\,\,\,\,\,\,\,\,'),kk, ' < 0')

This should gives this

image description

And now come my questions :

1) Can I assume $U^\prime>0$,$U^{\prime\prime}<0$ and verify with the other asumptions that the second order is verified, without hand manipulation. 2) How to rewrite the first order condition in such a way that all the terms in $U^\prime$ be on right side of the equation (a ratio) and the other terms on the left side ? 3) But, I think to know the answer, how to substitute $x$ to $I$ in the final stage ?