Help with iteration problem

1 view (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.

Categories

Find more on Mathematics in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!