Question: How to plot this Pde

restart;

PDEtools[declare]((f, g)(x), prime = x);

de1 := diff(f(x), x, x, x, x)-(H*H)*(diff(f(x), x, x))-R*(diff(f(x), x, x))*(diff(f(x), x))+R*(diff(f(x), x, x, x))*f(x);

de2 := diff(g(x), x, x)-(H*H)*g(x)-R*(diff(f(x), x))*g(x)+R*(diff(g(x), x))*f(x); R := 1; H := 1;

dd1 := {de1 = 0, de2 = 0, f(0) = 0, f(1) = 1, g(0) = 1, g(1) = 1, (D(f))(0) = 0, (D(f))(1) = 0};

r1 := dsolve(dd1, numeric, output = Array([0., 0.5e-1, .10, .15, .20, .25, .30, .35, .40, .45, .50, .55, .60, .65, .70, .75, .80, .85, .90, .95, 1.00]));

restart;

PDEtools[declare]((f, g)(x), prime = x);

de1 := diff(f(x), x, x, x, x)-(H*H)*(diff(f(x), x, x))-R*(diff(f(x), x, x))*(diff(f(x), x))+R*(diff(f(x), x, x, x))*f(x);

de2 := diff(g(x), x, x)-(H*H)*g(x)-R*(diff(f(x), x))*g(x)+R*(diff(g(x), x))*f(x); R := 5; H := 5;

dd1 := {de1 = 0, de2 = 0, f(0) = 0, f(1) = 1, g(0) = 1, g(1) = 1, (D(f))(0) = 0, (D(f))(1) = 0};

r2 := dsolve(dd1, numeric, output = Array([0., 0.5e-1, .10, .15, .20, .25, .30, .35, .40, .45, .50, .55, .60, .65, .70, .75, .80, .85, .90, .95, 1.00]));

odeplot({r1, r2}, [x, f(x)]);
print(`output redirected...`); # input placeholder
                f(x) will now be displayed as f
 

Please Wait...