Why does:
diffpressure = gradient(pressure)./gradient(density);
and
diffpressure = diff(pressure)./diff(density);
give drastically different results? The derivative using each method are completely different. See attached figure
differential equationsplot
diffpressure = gradient(pressure)./gradient(density);
diffpressure = diff(pressure)./diff(density);
Best Answer