function derTwoSided input Real u[:]; input Real y[size(u, 1)]; output Real der_y[size(u, 1)] "dy/du"; end derTwoSided;