fsolve函数怎样定义非线性函数方程。

局部程序如下：

syms y a

f1=-a^2*sinh(a*(-l))*[-sin(a*y),cos(a*y)/a-y*sin(a*y),cos(a*y),sin(a*y)/a+y*cos(a*y),0,0,0,0,0,0,0,0]*W

F1=int(f1,y,0,h1)

function f=myfun(a) %在这里就报错啦，如图，说是函数不能在此情况中定义，请问是什么情况？？

f=F1-q*l

a0=1;

options=optimset('Display','iter');

[x,fval]=fsolve(@myfun,a0,options)

W是一个矩阵，

f里面只含有参数a，可是f式子很长，超出了Matlab显示屏幕范围，所以我想调用fsolve函数求解参数a，望大神指导！