ciao ragazzi ho un problemino con matlab!
devo:
Sviluppare una function Matlab per determinare la soluzione del sistema di equazioni alle derivate parziali, su un segmento, di avvezione (1D) +diffusione (1D) + predatore – preda di due specie biologiche la cui concentrazione è c1(x,t) e c2(x,t):


dc1\dt +v*dc1\dx=K*d^2c1\dx^2 -ac1+ bc1c2

dc2\dt +v*dc2\dx=K*d^2c2\dx^2 -pc2- rc1c2




Utilizzare la tecnica del metodo delle linee senza splitting, con il metodo upwind per l’avvezione, centrale esplicito per la diffusione ed Eulero esplicito per preda predatore.

Utilizzare la function per risolvere un problema in cui le condizioni al contorno di inflow sono nulle, la condizione iniziale delle prede e dei predatori sono degli impulsi di valore massimo 10 e 5 (rispettivamente) centrati vicino all’estremo sinistro del dominio (estremo di inflow), la velocità è 0.2, il coefficiente di diffusione è 0.2, a=0.2, b=0.3, p=0.4 r=0.1.
Fare vari esperimenti variando alcuni parametri biologici.


avete qualche dritta???
come affrontare il problema delle differenze finite???
grazie!!!![/img]