Help with iteration problem

2 views (last 30 days)
mohsen
mohsen on 19 Nov 2012
>> x=0.7; % x is fixed value
>> syms s % s is symbols
>> v=char(x./(s^2+1)); % v is function of (x,s)
% numerical method to compute inversion Laplace transformation
% this function is available through remark at below
>> [t1,ft1]=INVLAP(v,0.0001,5,1000,6,39,89);
>> plot(t1,ft1)
please, help me the previous small code run easily, but i need (x not equal fixed), i.e., i need x- run over the interval (0:0.1:1), and also make all the following steps, then make 3-D figures between mesh(x,t1,ft1)
% remark
% INVLAP Numerical Inversion of Laplace Transforms
  1 Comment
mohsen
mohsen on 20 Nov 2012
thank u, i do the step which u refer to it in your answer but their exist error and the following message is appeared
??? Error using ==> mesh at 80
Z must be a matrix, not a scalar or vector.
code according to your vision
for x=0:0.1:1;
syms s % s is symbols
v=char(x./(s^2+1)); % v is function of (x,s)
% numerical method to compute inversion Laplace transformation
% this function is available through remark at below
[t1,ft1]=INVLAP(v,0.0001,5,1000,6,39,89);
end
mesh(x,t1,ft1)

Sign in to comment.

Answers (1)

bym
bym on 19 Nov 2012
seems like a simple 'for loop' would do the trick
x = 0:.1:1;
for x
v=char(x./(s^2+1)); % v is function of (x,s)
% modify call to INVLAP to store ft1 in a matrix
...
end
mesh(...)
  2 Comments
mohsen
mohsen on 20 Nov 2012
thank u, i do the step which u refer to it in your answer but their exist error and the following message is appeared
??? Error using ==> mesh at 80
Z must be a matrix, not a scalar or vector.
code according to your vision
for x=0:0.1:1;
syms s % s is symbols
v=char(x./(s^2+1)); % v is function of (x,s)
% numerical method to compute inversion Laplace transformation
% this function is available through remark at below
[t1,ft1]=INVLAP(v,0.0001,5,1000,6,39,89);
end
mesh(x,t1,ft1)
bym
bym on 20 Nov 2012
yes, that is why I put in my answer the comment
% modify call to INVLAP to store ft1 in a matrix

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!