MA 105 - KODER FOR MAPLE-FIGURER BRUKT I FORELESNINGENE

> with(plots):

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

> implicitplot3d({f(x,y,z) = 0, f(x,y,z) = 1, f(x,y,z) = -1}, x=-2..2, y=-2..2, z=-1..4);

>

> with(plots):

> f:= y^3 - x^2;

> contourplot(f, x=-2..2, y=-1.5..1.5);

> plot3d(f, x=-2..2, y=-1.5..1.5);

>

31/01

> with(plots);

> f:= 2*x^2*y/(x^4+y^2);

> plot3d(f(x,y), x=-0.5..0.5, y=-0.5..0.5, grid=[70,70]);

> contourplot (f(x,y), x=-0.5..0.5, y=-0.5..0.5, grid=[30,30]);

>

02/02

> with(plots):

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

> plot3d(f(x,y), x=-2..2, y=-2..2);

> contourplot(f(x,y), x=-2..2, y=-2..2);

> gradplot(f(x,y), x=-2..2, y=-2..2);

Error, (in plots/fieldplot/arrowsf2d) division by zero

> fieldplot([-2*x*exp(-x^2-y^2), -2*y*exp(-x^2-y^2)], x=-2..2, y=-2..2);

>

07-09/02

• with(plots):
• spacecurve([cos(t), sin(t), sin(2*t)], t=0..2*Pi);

• with(plots):
• f:= x^2 + y^2;
• contourplot(f(x,y), x=-2..2, y=-2..2);
• plot3d(f(x,y), x=-2..2, y=-2..2);
• fieldplot([2*x, 2*y], x=-2..2, y=-2..2);

• with(plots):
• f:= (-x^2 - y^2)*cos(x*y);
• plot3d(f(x,y), x=-1.5..1.5, y=-1.5..1.5);

16/02

• with(plots):
• f:= x*y;
• plot3d(f(x,y), x=-1.5..1.5, y=-1.5..1.5);

21/02

> with(plots):

> A:= plot3d(sqrt(x*x + y*y), x=-3..2, y=-4..2):

> B:= plot3d(-sqrt(x*x + y*y), x=-3..2, y=-4..2):

> C:= plot3d( (3-x-2*y)/3, x=-3..2, y=-4..2):

> display(A, B, C);

>

28/02

• with(plots):
• fieldplot([x,y], x=-1..1, y=-1..1);
• fieldplot([-y,x], x=-1..1, y=-1..1);
• fieldplot([-y/(x^2+y^2), x/(x^2+y^2)], x=-1..1, y=-1..1);
• fieldplot([ 1-y^2, 0], x=-1..1, y=-1..1);
• fieldplot([-x*(x^2+y^2)^(-3/2), -y*(x^2+y^2)^(-3/2)], x=-1..1, y=-1..1, grid=[10,10]);

• with(plots):
• fieldplot([x+y, sin(x)], x=-10..10, y=-7.5..7.5, grid=[10,10]);
• #strømlinjene til dette vektorfeltet ble plottet med McMac, programmet DiffEq, Phase Plane.

01/03

> with(plots):

> plot3d(1 + x^2 + y^2, x=0..1, y=0..1);

> with(linalg):

Warning, new definition for norm

Warning, new definition for trace

> A:= matrix(4, 4, (i,j) -> 1 + ((2*i - 1)/8)^2 + ((2*j - 1)/8)^2);

> matrixplot(A, heights=histogram, style=patch, shading=zgrayscale);

> B:= matrix(8, 8, (i,j) -> 1 + ((2*i - 1)/16)^2 + ((2*j - 1)/16)^2);

> matrixplot(B, heights=histogram, style=patch, shading=zgrayscale);

>

06/03

> with(plots):

> A:=plot3d([r*cos(t), r*sin(t), r*r*(1+sin(t)^2)], r=0..2, t=0..2*Pi):

> B:=plot3d([r*cos(t), r*sin(t), 12 - r*r*(1+cos(t)^2)], r=0..2, t=0..2*Pi):

> display(A,B);

>

08/03

> with(plots):

> A:=plot3d([cos(t), sin(t), u], t=0..2*Pi, u=-1.5..1.5):

> B:=plot3d([u, cos(t), sin(t)], t=0..2*Pi, u=-1.5..1.5):

> display(A,B);

>

> with(plots):

> A:=plot3d([u,v,2-u*u], u=-1..1, v=0..3):

> B:=plot3d([u,v,u*u], u=-1..1, v=0..3):

> display(A,B);

>

13/03

> with(plots):

> plot([(1+cos(t))*sin(t), (1+cos(t))*cos(t), t=0..Pi]);

> plot3d([1+cos(phi), theta, phi], theta=0..2*Pi, phi=0..Pi, coords=spherical);

>

20-22/03

> with(plots):

> plot3d([cosh(u)*cos(t),cosh(u)*sin(t)/2,sinh(u)], u=-1..1, t=0..2*Pi);

>

27/03

> with(plots):

> plot3d([r*cos(t), r*sin(t), r*r*sin(2*t)/2], r=0..1, t=0..2*Pi);

>

> with(plots):

> plot3d([4*cos(th)-t*cos(th)*sin(th/2), 4*sin(th)-t*sin(th)*sin(th/2), t*cos(th/2)], th=0..2*Pi, t=-1..1);

> #Møbius-flaten