[slgu]

这个问题事实上比较困难，试了一些常用的函数似乎都不太好。一种简单的近似可用：
fy=a*log(b*exp((x-c)).^2));