closed contour complex integral

2 次查看(过去 30 天)
salah zetreni
salah zetreni 2015-3-13
回答: Mike Hosea 2015-3-16
hello.. please help me ∮(sinπz^2+cosπz^2)/((z-1)(z-2)) dz around z=2
answer this integral by using matlab ?
My attempt at the bottom but it was wrong .. please help me to detection the error ..thank u
fun=@(z)(sin(pi*z.^2)+cos(pi*z.^2))/((z-1)*(z-2));
>> g=@(theta)2*cos(theta)+1i*2*sin(theta);
>> gprime=@(theta)-2*sin(theta)+1i*2*cos(theta);
>> q1=quad(@(t) fun(g(t)).*gprime(t),0,2*pi)
??? Error using ==> mtimes
Inner matrix dimensions must agree.

回答(1 个)

Mike Hosea
Mike Hosea 2015-3-16
MTIMES is matrix multiplication. The functional notation is mtimes(A,B), but the operator notation for matrix multiplication is A*B. So A*B assumes A and B have consistent matrix dimensions, i.e. size(A,2) must be the same as size(B,1), unless either is a scalar. You used .^ instead of ^, which is good. That's presumably because you meant element-wise powers of array inputs, but you did not always use .* instead of *, nor did you always use ./ instead of /. Go back and fix all of those operators.
The QUAD function (which is obsolete, BTW--use INTEGRAL or QUADGK instead), makes evaluations of your function in "batches", so your function must accept an array and return an array. It is not enough that it works only when you pass in a scalar. In other words f([1,2]) must work and return the same result as [f(1),f(2)].

类别

Help CenterFile Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息

标签

尚未输入任何标签。

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by