Thank you for this! I did not know it was this straight foward.
Although i am trying to understand explicitly what this part of the code does. maybe you can break it down for me.
getVals:= (V, space, time)-> rhs(pdsA[V]:-value( w(x,t))(space,time)):
for k in Vval do
Matrix(11, 11, (i,j)-> getVals( k, (i-1)/10, (j-1)/10));
pdsA[k]:-plot3d( w(x,t), x=0..1, t=0..1);
You did state that the *best* method to do this probably depends on what I subsequently want to do with the resulting data. Let me explain what I am having in mind. So i intend to use the values of the displacment w(x,t) to solve for frequency using the formula w(x,t)/ (A*cos2tPi). where A is the amplitude (maximum displacement attained in w) and t is the time.
This would allow me get values for frequency F at the varying times. of which i would use to plot those values against the velocity V range.
I hope this gives a picture of the task i am trying to accomplish. Thank you