I am trying to compute the Laplace transform of a known function. I have replicated the code shown below in Mathematica and it works fine, but I am unsure as to why it is not working in Matlab. The result I get is just laplace (some expression).
syms G w0 wc wk tau s C
a(tau) = int(wk^3*exp(-1i*(wk-w0)*tau), wk, 0, wc);
corr(tau) = (G/(w0^3))*(a(tau));
F(s) = laplace(corr(tau),tau,s);
C(s) = solve(s*C-1==-F(tau)*C,C);