Question: How can i obtain numerical data instead of of plot using bvp[midrich]?

hi all..

i m solving a system of four coupled differential equations using bvp[midrich] technique on maple..

i have two questions..

1). is it possible to calculate the numerical data using bvp[midrich] instead of graphs?

2). there is a physical constant "A" in my second equation..can i compute its numerical values while fixing all other physical parameters of my equations as fixed?,,, i m writting my system of equations as under:

> eq1 := diff(f(eta),eta,eta,eta)+(1+lambda)*(f(eta)*diff(f(eta),eta,eta)-(diff(f(eta),eta))^2+(B)^2)+k*(((diff(f(eta),eta,eta))^2-f(eta)*diff(f(eta),eta,eta,eta,eta)));

> eq2 := diff(h(eta),eta,eta)+(1+lambda)*(f(eta)*diff(h(eta),eta)-diff(f(eta),eta)*h(eta)-A)+k*((diff(f(eta),eta,eta)*diff(h(eta),eta)-diff(f(eta),eta)*diff(h(eta),eta,eta)+diff(f(eta),eta,eta,eta)*h(eta)-f(eta)*diff(h(eta),eta,eta,eta)));

> eq3 := diff(theta(eta),eta,eta)+Pr*(f(eta)*diff(theta(eta),eta)+`N[b]`*diff(phi(eta),eta)*diff(theta(eta),eta)+`N[t]`*(diff(theta(eta),eta))^2);

> eq4 := diff(phi(eta),eta,eta)+S[c]*f(eta)*diff(phi(eta),eta)+`N[t]`/`N[b]`*(diff(theta(eta),eta,eta));

> bc:=f(0)=0,D(f)(0)=1,D(f)(N)=B,D(D(f))(N)=0,h(0)=0,D(h)(N)=1,D(D(h))(N)=0,theta(0)=1,theta(N)=0,phi(0)=1,phi(N)=0;

(2)
bc := f(0) = 0, D(f)(0) = 1, D(f)(10) = 0.1, (D )(f)(10) = 0,

(2)
h(0) = 0, D(h)(10) = 1, (D )(h)(10) = 0, theta(0) = 1,

theta(10) = 0, phi(0) = 1, phi(10) = 0

> N:=10;lambda:=0.1;A:=0.84;B:=0.1;k:=0.1;`N[b]`:=0.1;`N[t]`:=0.1;Pr:=10;S[c]:=1;

N := 10

lambda := 0.1

A := 0.84

B := 0.1

k := 0.1

N[b] := 0.1

N[t] := 0.1

Pr := 10

S[c] := 1

> lambda := .1;

lambda := 0.1

> k := .1;

k := 0.1

> A := .84;

A := 0.84

> A1:=dsolve({eq1,eq2,eq3,eq4,bc},numeric,method=bvp[midrich]);

A1 := proc(x_bvp) ... end proc

> with(plots):
> p1:=odeplot(A1, [[eta,diff(f(eta),eta), color=red,style=line,thickness=2]],0..N):
> display({p1}, axes=boxed, title=`Velocity f'`);

> q1:=odeplot(A1, [[eta,diff(h(eta),eta), color=blue,style=line,thickness=2]],0..N):
> display({q1}, axes=boxed, title=`Velocity h'`);

> r1:=odeplot(A1, [[eta,theta(eta), color=green,style=line,thickness=2]],0..N):
> display({r1}, axes=boxed, title=`temp`);

> s1:=odeplot(A1, [[eta,phi(eta), color=orange,style=line,thickness=2]],0..N):
> display({s1}, axes=boxed, title=`concentration`);

waiting for a good response