1fx(x):=erfc(sin(x));
2dfx(x):=diff(fx(x),x);
3df2x(x):=diff(dfx(x),x);
4df3x(x):=diff(df2x(x),x);
5Fortran(expr):= fortran(ev(expr,numer));
6with_stdout("submaxima_erfc.F",
7  print("      subroutine submaxima_erfc(x,f,df,df2,df3)"),
8  print("      implicit none"),
9  print("      double precision x,f,df,df2,df3"),
10  Fortran(f = fx(x)),
11  Fortran(df = dfx(x)),
12  Fortran(df2 = df2x(x)),
13  Fortran(df3 = df3x(x)),
14  print("      end"));
15