here is the paper: www.amleto.pwp.blueyonder.co.uk/paper.pdf
I am tryng to reproduce figures 2 and 4.
my scripts they are quite simple and easy to compare with the paper. They are at the bottom of this post :)
I believe the author has put the inital condition sin and cos the wrong way around...
dot product on the rows of A:
z.*A(i,:)
dot product on the columns of A:
z.*A(:,i)
If you want to keep z as a row/column vectorI think you will still need to loop, Although you could make a temporary variable that was padded with lots of z's...
ztemp=zeros(size(A))
ztemp(1:size(A,1),:)...
hi, thanks a lot for the help.
dr(5) is a universal quantity and doesn't depend on the bubbles so maybe I should make this dr(1), and then the start of the bubble quantities would be 2,8,12 etc?
Would this require all r(x:5:z) to be changed to r(x:4:z)?
I would have thought...
ok, the deal is - I have a few bubbles. They pulsate due to an externally changing applied pressure. The pulsations and the external pressure make the bubbles move, and the following code attempts to model that for two bubbles.
function dr = trans1func(t,r)
global sPeq sPatm sP0 sPv Re We dk...
I am currently using ode45 but am a bit stuck.
What I need to do is basically progress the solution of two instances of the same equation, with different initial conditions, at the same time.
I could just write out the equation again in the function but eventually I may be dealing with 50 of...
