Hello every one,
I am gettting this error "unable to store HFloat", while solving a system of ODES numerically.
restart:
eq1:=n*(-diff(f(eta),eta$2))^(n-1)*diff(f(eta),eta$3)+(2*n/(n+1))*f(eta)*diff(f(eta),eta$2)
-diff(f(eta),eta$1)^2-M*diff(f(eta),eta$1)+M*epsilon+epsilon^2+lambda*f1(eta)=0;
eq2:=diff(f1(eta),eta$2)+P*(2*n/(n+1))*(f(eta)*diff(f1(eta),eta$1))=0;
bc:=f(0)=0,D(f)(0)=1,D(f1)(0)=1,D(f)(N)=epsilon,f1(N)=0;
N:=10:para:={P=0.7,lambda=1,n=0.4,M=2};
dsolve(subs(para,epsilon=0.6, {bc, eq1, eq2}), numeric, maxmesh = 35000,
output = array([seq(0.1e-1*i, i = 0 .. 100*N)])):
Error, (in dsolve/numeric/BVPSolve) unable to store 'HFloat(0.03902115204890281)-
HFloat(0.12009475726252063)*I' when datatype=float[8]
Please help me out.
Thanks
