janhardo

910 Reputation

13 Badges

11 years, 312 days
B. Ed math

MaplePrimes Activity


These are replies submitted by janhardo

Improved FastPursuit2 procedure

FastPursuit2(5, 1.5, 1.0, 0, 2.0, 0., 1200, 0.05) example: the wolf is running in a almost concentric circle at almost fixed distant (not this plot)

achtervolgings_kromme_wolf_en_geit_mprimes_DEF_14-9-2025_.mw

achtervolgings_kromme_wolf_en_geit_mprimes_14-9-2025_.mw

Making now a FastPursuit2  procedure using dsolve :
sol := dsolve(sys union ics, {xw(t), yw(t)}, numeric, method = rkf45)

@Alfred_F 
 

Starting position and speed are decisive?

The goat continues to move in a circle

The wolf remains focused on the goat

Was it not with a pursuit curve?

Bring back  the number of sliders ? 

@Alfred_F 
Yes, that is a proof by contradiction, looks like on the proof voor the irrationality of square root two number
aside:
The ancient Greeks were confronted with the existence of the square root of 2 and believe that it came as a shock to their world view at that time.

periodiciteit_cos_achtergronden_mprimes7-9-2025.mw

@vv 
In my opinion, this is a mathematical proof sketch and the correct details should be filled in.
Getting this expression for cos(2pi.n) seems to be straight forwards ?, but now further..?

Was @Alfred not reffering to the function  exp(x) ?
bewijs_perodiciteit_mprimes_4-9-2025.mw
y=cos(x) is not bijective 

I can enlarge screen in windows 11 pro  with Ctrl + Middle mouse wheel ?

@Alfred_F 
This is userfriendlier by using y(x) explicit 

NULL

with(Student:-ODEs)

ode1 := y(x)*(diff(y(x), x)) = a

y(x)*(diff(y(x), x)) = a

(1)

ODESteps(ode1)

"[[,,"Let's solve"],[,,y(x) ((ⅆ)/(ⅆx) y(x))=a],["•",,"Highest derivative means the order of the ODE is" 1],[,,(ⅆ)/(ⅆx) y(x)],["•",,"Integrate both sides with respect to" x],[,,∫y(x) ((ⅆ)/(ⅆx) y(x)) ⅆx=∫a ⅆx+`c__1`],["•",,"Evaluate integral"],[,,((y(x))^2)/2=a x+`c__1`],["•",,"Solve for" y(x)],[,,{y(x)=sqrt(2 a x+2 `c__1`),y(x)=-sqrt(2 a x+2 `c__1`)}],["•",,"Redefine the integration constant(s)"],[,,{y(x)=sqrt(2 a x+`c__1`),y(x)=-sqrt(2 a x+`c__1`)}]]"

(2)
 

restart

"maple.ini in user"

(3)

with(Student:-ODEs)

ODE2 := {y(x)*(diff(y(x), x)) = a, y(0) = b}

{y(x)*(diff(y(x), x)) = a, y(0) = b}

(4)

ODESteps(ODE2)

"[[,,"Let's solve"],[,,{y(x) ((ⅆ)/(ⅆx) y(x))=a,y(0)=b}],["•",,"Highest derivative means the order of the ODE is" 1],[,,(ⅆ)/(ⅆx) y(x)],["•",,"Integrate both sides with respect to" x],[,,∫y(x) ((ⅆ)/(ⅆx) y(x)) ⅆx=∫a ⅆx+`c__1`],["•",,"Evaluate integral"],[,,((y(x))^2)/2=a x+`c__1`],["•",,"Solve for" y(x)],[,,{y(x)=sqrt(2 a x+2 `c__1`),y(x)=-sqrt(2 a x+2 `c__1`)}],["•",,"Redefine the integration constant(s)"],[,,{y(x)=sqrt(2 a x+`c__1`),y(x)=-sqrt(2 a x+`c__1`)}],["•",,"Use initial condition" y(0)=b],[,,b=sqrt(`c__1`)],["•",,"Solve for" `c__1`],[,,`c__1`=b^2],["•",,"Substitute" `c__1`=b^2 "into general solution and simplify"],[,,y(x)=sqrt(2 a x+b^2)],["•",,"Use initial condition" y(0)=b],[,,b=-sqrt(`c__1`)],["•",,"Solve for" `c__1`],[,,`c__1`=b^2],["•",,"Substitute" `c__1`=b^2 "into general solution and simplify"],[,,y(x)=-sqrt(2 a x+b^2)],["•",,"Solutions to the IVP"],[,,[y(x)=sqrt(2 a x+b^2),y(x)=-sqrt(2 a x+b^2)]]]"

(5)
 

with(plots); p1 := plot(sqrt(2*x), x = 0 .. 5, color = red, thickness = 2, legend = "y = √(2x), a=1, b=0"); p2 := plot(sqrt(x), x = 0 .. 5, color = green, thickness = 2, legend = "y = √(x), a=0.5, b=0"); p3 := plot(sqrt(-2*x), x = -5 .. 0, color = blue, thickness = 2, legend = "y = √(-2x), a=-1, b=0"); display([p1, p2, p3], title = "Solutions of y(x) = sqrt(2*a*x + b²) with b = 0 (positive branch)", labels = ["x", "y(x)"], axes = boxed)

 
 

``

Download ode_via_odesteps_mprimes24-8-2025.mw

The dependent variable cannot be easily isolated. ?

If b= 0 in ic y(0)=0

y1 := sqrt(2*x);
y2 := sqrt(2*0.5*x);
y3 := sqrt(2*(-1)*x);
p1 := plot(y1, x = 0 .. 5, color = red, thickness = 2, legend = "a = 1");
p2 := plot(y2, x = 0 .. 5, color = green, thickness = 2, legend = "a = 0.5");
p3 := plot(y3, x = -5 .. 0, color = blue, thickness = 2, legend = "a = -1");
display([p1, p2, p3], title = "Solutions of y(x) = sqrt(2*a*x) with Initial Condition y(0) = 0", labels = ["x", "y(x)"], axes = boxed);


@salim-barzani 
I don't think I fully understand you.

@salim-barzani 

Do you want to calculate all this manually?

Then you need to extract the procedures from the wkptest module.

The procedures must be reduced from general to a single calculation.

First 9 10 11 12 13 14 15 Last Page 11 of 85