Advanced Search

Question: How to plot residual error

Question:How to plot residual error

666 jvbasha 130 Maple 18
ode dsolve numeric differential-equations bvp duplicate-question
May 03 2020
1

Dear maple users,
Greetings.
How to plot residual error for BVP.
Here I have enclosed the file.rerror.mw
 

> 

restart; with(plots)
> 

fcns := {f(eta), g(eta)}:
> 

bet := 0.:
> 

eq1 := diff(f(eta), `$`(eta, 3))+(diff(f(eta), `$`(eta, 2)))*f(eta)+be*((diff(g(eta), `$`(eta, 1)))^2-(diff(g(eta), `$`(eta, 2)))*g(eta)) = 0:
> 

eq2 := pr*lam*(diff(g(eta), `$`(eta, 3)))+(diff(g(eta), `$`(eta, 2)))*f(eta)-(diff(f(eta), `$`(eta, 2)))*g(eta) = 0:
> 

bc := f(0) = 0, (D(f))(0) = 0, (D(f))(N) = 1, g(0) = 0, ((D@@2)(g))(0) = 0, (D(g))(N) = 1;

f(0) = 0, (D(f))(0) = 0, (D(f))(5) = 1, g(0) = 0, ((D@@2)(g))(0) = 0, (D(g))(5) = 1

 (1)
> 

R := dsolve(eval({bc, eq1, eq2}), fcns, type = numeric, method = bvp[midrich], output = listprocedure):
> 

Sol_exact := dsolve(eval({bc, eq1, eq2}), fcns):
> 

``


 

Download rerror.mw

Have a good day

 
Please Wait...