Jakubi

The suggestions work perfectly.

Many thanks!

Alec and Jakubi,

I apologize for not defining "l" in the parameters. Thanks for catching that error.

Jakubi,

Thanks for the suggestions. I am sure that will do the trick. The code generation is for generating a Fortran form of the eventual equation. So it is kinda needed but not necessarily at this moment.

I will try these suggestions as soon as I go back to the lab.

Appreciate your help!

Alec,

Its the letter "L" albeit in small caps in the definition of "b".

Thanks guys!

Thanks Jakubi!

You nailed the issue correctly. The parameters were becoming function calls and as such Maple was unable to figure it out. Unfortunately now when I run the problem after about 500-700 secs. when it says kernel lost.

Here are the command lines that I run (basically I changed the boundary conditions and updated the dsolve command in the end):

> with(CodeGeneration);

> with(PDEtools);

> declare(prime = x);

> b := ap/(l+si)^1.5;

> NULL;

> q := proc (x) options operator, arrow; ap/(x+si)^1.5-b end proc;

> bc1 := w(0) = 0;

> E*MI*(diff(w(x), `$`(x, 3)))/(A*G)+(diff(q(x), x))/(A*G)^2+diff(w(x), x); convert(%, D); eval(%, x = 0);

> -E*MI*(diff(w(x), `$`(x, 3)))/(A*G)-E*MI*(diff(q(x), x))/(A*G)^2; convert(%, D); eval(%, x = 1);

> diff(w(x), `$`(x, 2))+q(x)/(G*A); convert(%, D); eval(%, x = 1);

> myode := E*MI*(diff(w(x), `$`(x, 4)))+E*MI*(diff(q(x), `$`(x, 2)))/(G*A) = q(x);

> mysol1 := rhs(dsolve({bc1, myode, E*MI*((D@@3)(w))(0)/(A*G)+(D(q))(0)/(A^2*G^2)+(D(w))(0) = 0, -E*MI*((D@@3)(w))(1)/(A*G)-E*MI*(D(q))(1)/(G^2*A^2) = 0, ((D@@2)(w))(1)+q(1)/(G*A) = 0}, w(x)));

I figure there might be inefficiencies or redundancy issues with my method of writing the commands. Do you have any suggestions to improve it? Is there a way to increase memory allocation for the problem? I am running it on a machine which has 2 GB RAM and 2.13 Ghz processor.

with regards,

lawkc.