%TFY4108 Fysikk hoesten 2014
%Oving 6, oppgave 5: 
%Plotting av potensialfunksjonen
%U(x,y) = x^2*y^2
%og det tilhoerende kraftfeltet F(x,y)
%med komponentene
%Fx = -2*x*y^2 og Fy = -2*x^2*y
%Her har vi satt U0=1 og R=1
clear all;
figure(1);
domain = [-2,2,-2,2];  %Syntaks: [xmin,xmax,ymin,ymax]
ezsurf('x^2*y^2',domain);
figure(2);
[x,y] = meshgrid(-2:0.5:2);
Fx = -2.*x.*y.*y;
Fy = -2.*x.*x.*y;
quiver(x,y,Fx,Fy);