function vectorField2( f ,g, It, Ix, Iy  )
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
t=It(1):0.25:It(2);
x=Ix(1):0.25:Ix(2);
y=Iy(1):0.25:Iy(2);
[T,X,Y]=meshgrid(t,x,y);
vx=cos(atan(g(T,X,Y)./f(T,X,Y)));
vy=sin(atan(g(T,X,Y)./f(T,X,Y)));
subplot(2,2,1)
quiver(X,Y,vx,vy)
%plano t-x
vt1=cos(atan(f(T,X,Y)));
vx1=sin(atan(f(T,X,Y)));
subplot(2,2,3)
quiver(T(:,:,30),X(:,:,30),vt1(:,:,30),vx1(:,:,30))
%plano t-y
y(30)
vt2=cos(atan(g(T,X,Y)));
vy2=sin(atan(g(T,X,Y)));
subplot(2,2,4)
quiver(T,Y,vt2,vy2)
end