clf
[T,X]=meshgrid(t,x);
Z=f(T,X);
%flipud(Z)   % havainnollistukseen, tosi käytössä kommentti tai puolipiste
%plot(T(:),X(:),'x')
hold on
U=ones(size(Z));
normit=sqrt(1+Z.^2);
quiver(T,X,0.5*h*U./normit,0.5*h*Z./normit,0)
axis equal
axis([a-h,b+h,c-h,d+h])
grid
xlabel('t');ylabel('x')
%title('x''=tx')
shg
