Numerical integration for a vector.

Can anybody help me with this? I do have problems to be integrating, when z is also defined as a vector, so I can't be using syms.

5 Comments

How are z and d related?
What is the expected result?
This is the exppected results, where z is a vector and d = length(z)/2
And n should be a very large value, example 0:1:1e4.
And at t=0, then the higest value of u_e is equal u_i, where u_i is a constant.
Try either of these, after supplying , , and t
z = ...;
d = numel(z)/2;
n = linspace(0, 1E4, 2*d);
ftrpz = (ui,cv,t,z) trapz(ui.*sin(pi*n.*z/(2*d))) .* sin(pi*n.*z/(2*d)) .* exp(-n.^2.*pi^2.*cv.*t/(4*d^2));
fintg = (ui,cv,t,z) integral(@(z)ui.*sin(pi*n.*z/(2*d)), 0, inf, 'ArrayValued',1) .* sin(pi*n.*z/(2*d)) .* exp(-n.^2.*pi^2.*cv.*t/(4*d^2))
See the documentation on trapz and integral for details. Then decide which version of ‘f’ works best in your application. Use the sum function complete the expression.
I cannot test this, so I am posting it as UNTESTED CODE. It will be necessary to experiment with it with your data to produce the desired result. The variable does not appear anywhere, so I have no idea what to do with it.
There is not enough information to proceed further.

Sign in to comment.

Answers (0)

Asked:

on 17 Apr 2022

Commented:

on 17 Apr 2022

Community Treasure Hunt

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

Start Hunting!