Question: Cut of Errorfunction

I'm currently wondering about the cut I'm looking for in the following worksheet.

I evaluate it in 2 ways but get different answers. Any idea what the problem here is?

Thanks


 

restart; dIs := sqrt(Pi/(I*s))*exp(I*s*t-I*s*omega0^2); Is1 := `assuming`([simplify(int(dIs, s))], [s > 0]); dIs := `assuming`([int(exp(-I*(omega^2+omega0^2-t)*s), omega = -infinity .. infinity)], [s > 0]); Is2 := int(%, s); plot3d(Im(eval(Is1, [t = x+I*y, s = 1, omega0 = 1])), x = -3 .. 3, y = -3 .. 3)

(-I*Pi/s)^(1/2)*exp(I*s*t-I*s*omega0^2)

 

(1/2-(1/2)*I)*Pi*2^(1/2)*erf(s^(1/2)*(I*(omega0^2-t))^(1/2))/(I*(omega0^2-t))^(1/2)

 

exp(I*s*t-I*s*omega0^2)*Pi^(1/2)/(I*s)^(1/2)

 

-I*Pi*erf((omega0^2-t)^(1/2)*(I*s)^(1/2))/(omega0^2-t)^(1/2)

 

 

``


 

Download CutErrorFunction.mw

Please Wait...