In [26]:

    

var('p,q,r')


    

    
Out[26]:





(p, q, r)

In [27]:

    

def SurfArea(a,b,c):
    return (1/4)*sqrt(2*a*a*b*b+2*b*b*c*c+2*c*c*a*a-a*a*a*a-b*b*b*b-c*c*c*c)

F=SurfArea(1,1,1)


    

In [28]:

    

implicit_plot3d(max_symbolic(min_symbolic(max_symbolic(F*(q*q+r*r-p*p),SurfArea(p,q,r)*(1+1-1))+
                            max_symbolic(F*(-q*q+r*r+p*p),SurfArea(p,q,r)*(1+1-1))-2*F*1*r,
                            max_symbolic(F*(r*r+p*p-q*q),SurfArea(p,q,r)*(1+1-1))+
                            max_symbolic(F*(-r*r+p*p+q*q),SurfArea(p,q,r)*(1+1-1))-2*F*1*p,
                            max_symbolic(F*(p*p+q*q-r*r),SurfArea(p,q,r)*(1+1-1))+
                            max_symbolic(F*(-p*p+q*q+r*r),SurfArea(p,q,r)*(1+1-1))-2*F*1*q,
#                           max_symbolic(F*(p*p+r*r-q*q),SurfArea(x,y,z)*(1+1-1))+
#                           max_symbolic(F*(-p*p+r*r+q*q),SurfArea(x,y,z)*(1+1-1))-2*F*1*r,
#                           max_symbolic(F*(r*r+q*q-p*p),SurfArea(x,y,z)*(1+1-1))+
#                           max_symbolic(F*(-r*r+q*q+p*p),SurfArea(x,y,z)*(1+1-1))-2*F*1*q,
#                           max_symbolic(F*(q*q+p*p-r*r),SurfArea(x,y,z)*(1+1-1))+
#                           max_symbolic(F*(-q*q+p*p+r*r),SurfArea(x,y,z)*(1+1-1))-2*F*1*p,                             
                            ),
                            r-p-q,p-r-q,q-r-p
                            ),
                (p, 0, 1.2), (q, 0, 1.2), (r, 0, 1.2),plot_points=100)


    

    
Out[28]:

In [ ]: