thanks again this exactly what I want.... thank u

thanks again this exactly what I want.... thank u

can I ask you another question?

now I want the solution for specific X values so I defined An array, and solved for it. it gives me the values for (x,f,f',f'') is there any way to print out just the value of (x,f')?

my code is:

> restart; interface(rtablesize = 45);
> ode := diff(f(x), x, x, x)+f(x)*(diff(f(x), x, x))+B*(1-(diff(f(x), x))^2) = 0;
> bcs := f(0) = 0, (D(f))(0) = 0, (D(f))(5) = 1;
> sys := {bcs, ode};
> A := Array(Vector[row](28, {(1) = 0, (2) = .1, (3) = .2, (4) = .3, (5) = .4, (6) = .5, (7) = .6, (8) = .7, (9) = .8, (10) = .9, (11) = 1, (12) = 1.2, (13) = 1.4, (14) = 1.6, (15) = 1.8, (16) = 2, (17) = 2.2, (18) = 2.4, (19) = 2.6, (20) = 2.8, (21) = 3, (22) = 3.2, (23) = 3.4, (24) = 3.6, (25) = 3.8, (26) = 4, (27) = 4.5, (28) = 5}));
> for b in [-.199, -.18, 0, .3, 1, 2, 10] do 

sys0 := eval(sys, B = b);

 dsn := dsolve(sys0, numeric, output = A);

 print(dsn) end do;

Thanks again

can I ask you another question?

now I want the solution for specific X values so I defined An array, and solved for it. it gives me the values for (x,f,f',f'') is there any way to print out just the value of (x,f')?

my code is:

> restart; interface(rtablesize = 45);
> ode := diff(f(x), x, x, x)+f(x)*(diff(f(x), x, x))+B*(1-(diff(f(x), x))^2) = 0;
> bcs := f(0) = 0, (D(f))(0) = 0, (D(f))(5) = 1;
> sys := {bcs, ode};
> A := Array(Vector[row](28, {(1) = 0, (2) = .1, (3) = .2, (4) = .3, (5) = .4, (6) = .5, (7) = .6, (8) = .7, (9) = .8, (10) = .9, (11) = 1, (12) = 1.2, (13) = 1.4, (14) = 1.6, (15) = 1.8, (16) = 2, (17) = 2.2, (18) = 2.4, (19) = 2.6, (20) = 2.8, (21) = 3, (22) = 3.2, (23) = 3.4, (24) = 3.6, (25) = 3.8, (26) = 4, (27) = 4.5, (28) = 5}));
> for b in [-.199, -.18, 0, .3, 1, 2, 10] do 

sys0 := eval(sys, B = b);

 dsn := dsolve(sys0, numeric, output = A);

 print(dsn) end do;

Thanks again