%**************************************************************** % File: FLUSSO.M % % Scopo: Creazione linee di campo di flusso dell' equazione % y'=f(x,y) % % Input: ax,bx estremi intervallo x % ay,by " " y % nx,ny numero intervalli in x, y % f stringa contenente la funzione f(x,y). % % Output: % Richiami: % Commenti: % % Autori: mf %***************************************************************** % hx=(bx-ax)/(nx-1); hy=(by-ay)/(ny-1); h=min(hx,hy); % lunghezza frecce xx=[ax:hx:bx]; yy=[ay:hy:by]; axis([ax bx ay by]) plot(ax,bx,'x',bx,by,'x') hold on for i=1:nx x=xx(i); xp(1)=x; for j=1:ny y=yy(j); ff=eval(f); t=atan(ff); yp(1)=y; xp(2)=x+h*cos(t); yp(2)=y+h*sin(t); plot(xp,yp,xp(2),yp(2),'+') end end title(' Campo di flusso') xlabel(['f(x,y) = ' f]) hold off