Buongiorno a tutti,
sto cercando di implementare in Matlab il metodo delle differenze finite centrato con le condizioni di Neumann della seguente equazione:
$\frac{\partial u(x,t)}{\partial t}-B \frac{\partial u(x,t)}{\partial t^2}=f$ dove $f(x,t)=-x^2+1-2B+2Bt$.
Le condizioni iniziali e di Neumann sono le seguenti:
$u(x,0)=x^2-1$ e $\frac{\partial u}{\partial x}(x_0,t_n)=\frac{u_1^n-u_0^n}{h}$.
Non riesco a capire come devo fare ad implementare il codice su Matlab, in particolare come scrivere le condizioni di Neumann
Grazie a tutti coloro che mi aiuteranno!!