x=indgen(200)*0.01
c=[3.0,4.0,0.5]
y=mysin(x,c)+randomn(seek,200)

d=[3.0,2.0,2.0]
p=mpfitfun('mysin',x,y,0.1,d)


plot,x,y
oplot,x,mysin(x,c),color=120
oplot,x,mysin(x,p)