Vrstevnice

> restart;

> with(plots):

> f:=(x,y)->x^2+y^2;

zobrazeni fce f(x,y)

> P1:=plot3d(f(x,y), x=-3..3, y=-sqrt(9-x^2)..sqrt(9-x^2), axes=framed, tickmarks=[7,7,5],labels=[x,y,z],style=patchcontour):

zobrazeni roviny z=0 ..8

> p_2:=seq(plot3d(k,x=-3..3,y=-3..3, style=patchnogrid),k=0..8):

> p_22:=display(p_2, axes=framed, tickmarks=[7,7,5],labels=[x,y,z],insequence=true):

zobrazeni vrstevnic v rovine z=0..8

> P_3:=seq(spacecurve([sqrt(k)*cos(t), sqrt(k)*sin(t), k], t=0..2*Pi, color=red, thickness=3),k=0..8):

> p_33:=display(P_3,insequence=true):

> for i from 0 to 8 do

> for j from 0 to 8 do

> if j<i then

> v[i][j]:=spacecurve([sqrt(0)*cos(t), sqrt(0)*sin(t), 0], t=0..2*Pi, color=black, thickness=2);

> else

> v[i][j]:=spacecurve([sqrt(i)*cos(t), sqrt(i)*sin(t), 0], t=0..2*Pi, color=black, thickness=2);

> end if:

> end do:

> vr[i]:=seq((v[i][l]),l=0..8);

> disp[i]:=display(vr[i],insequence=true);

> end do:

>

> P_4:=seq(spacecurve([sqrt(l)*cos(t), sqrt(l)*sin(t), 0], t=0..2*Pi, color=black, thickness=3),k=0..8):

> p_44:=display(P_4,insequence=true):

>

zobrazeni textu rovina c=0..8

> tx:=seq(textplot3d([[3,-3,-2,`c=`],[2.5,-2.5,-2,k]],color=blue,font=[TIMES,BOLD,20]),k=0..8):

> txt:=display(tx,insequence=true):

>

> display(P1,p_22,p_33,p_44,txt,disp[0],disp[1],disp[2],disp[3],disp[4],disp[5],disp[6],disp[7],disp[8], axes=normal,orientation=[55,70]);

> display(P1,p_22,p_33,axes=normal,orientation=[0,0],scaling=constrained);

> display(P1,p_22,p_33,axes=normal,orientation=[90,180],scaling=constrained);