Tečná rovina

1
y=var('y')
2
def tecna_rovina(FUN,bod,x_hr,y_hr):
3
    Gh=plot3d(FUN,(x,x_hr[0],x_hr[1]),(y,y_hr[0],y_hr[1]),opacity=0.8,color='blue')
4
    H0=point([bod[0],bod[1],FUN(bod[0],bod[1])],color='yellow',size=10)
5
    hy(x,y)=diff(FUN,y)
6
    hx(x,y)=diff(FUN,x)
7
    z=FUN(bod[0],bod[1])+hx(bod[0],bod[1])*(x-bod[0])+hy(bod[0],bod[1])*(y-bod[1])
8
    rovina=plot3d(z,(x,x_hr[0],x_hr[1]),(y,y_hr[0],y_hr[1]),opacity=0.8,color='grey')
9
    pretty_print("Rovnice tečné roviny:", "z= ",z)
10
    show(Gh+H0+rovina)
11
    
12
#f(x,y)=x^2-y^2
13
#tecna_rovina(FUN=f,bod=[-1,0], x_hr = [-3,3], y_hr=[-3,3])
14