janhardo

715 Reputation

12 Badges

11 years, 101 days

MaplePrimes Activity


These are replies submitted by janhardo

A procedure , numeric, symbolic sphere equation for given sphere centre and 3 circles 

A phytagoras circles example

res := SphereFromSections([8, 4, 0], 4, 3, 5, symbolic = true, plot = true);
res[equation];
res[center], res[radius];

"Sphere equation: ", (x - 3.000000000)^2 + (y - 4.000000000)^2 + z^2 = 25.00000000

opgave_rouben_bol_mprimes_procedure_nmeriek_en_symbolisch_11--8-2025DEF.mw

@Rouben Rostamian  
Ai says  , its ill posed?
At first, the Ai said that the task was impossible, but later AI said that it was possible.( again question) 



.

Download opgac_rouben_bol_mprimes_11--8-2025.mw

 

@sand15 
The ai  is improved lately, and seems to be much better , making this is not a achievement
bolvergelijking_opstellen_mprimes_10-8-2025.mw

@AHSAN 

odeplot chooses its own (adaptive or fixed) sample points to draw a nice curve.

Difficulties with interpolation in dsol

Maple stores the numerical solution on an adaptive grid (mesh) — not on my fixed 9 y-values that I want.

That is why I get → Float(undefined) when I try calculate points.

There is no explicit function that is so easily known for the 5 graphs of the fraction coefficient to obtain the plot values.

@AHSAN 
 

That code does not work because the two tables are not filled.

The plot with 5 curves is drawn, so there are known y values, but the current code fails to place function values in the two tables shown.

@AHSAN 

It seems more practical to first create a function value table for 1 plot.

Does Maple have a procedure for this already ?

@AHSAN 
Maple 2025.1  gives right output table structure , but still not complete

@AHSAN 

I enclosed the same file again yesterday ( I must delete post ?) , so i go further with this one (intented for yesterday.). Must be further improved.
 tabel_skin_fraction_coeifficent_mprimes_A_7-8-2025.mw

@AHSAN 
Is the intention to view the function values in a table format?

@AHSAN 
Good day , the Excle export is now only in one column, must be further  improved ..
If you are on Windows computer , then chance username ( here it is : jando ) in your own username. 
nusselt-brinkmann_getallen_en_k_waarde_in_excel_5-8-2025.mw 

@lucaud 
What do you want to achieve now?

The system of odes is now for first order solved and it seems be better for handling?

 

restart; with(plots); with(DEtools); M := 2; m := 1; R := 1; r := .3; g := 9.81; sys := {diff(x1(t), t) = x2(t), diff(x2(t), t) = (-g*(M+m)*sin(x1(t))-m*r*x4(t)^2*sin(x1(t)-x3(t))+m*(g*sin(x3(t))-x2(t)^2*(R-r)*sin(x1(t)-x3(t)))*cos(x1(t)-x3(t)))/(2*M*R-2*M*r-R*m*cos(x1(t)-x3(t))^2+R*m+m*r*cos(x1(t)-x3(t))^2-m*r), diff(x3(t), t) = x4(t), diff(x4(t), t) = (-(2*M+m)*(g*sin(x3(t))-x2(t)^2*(R-r)*sin(x1(t)-x3(t)))+(g*(M+m)*sin(x1(t))+m*r*x4(t)^2*sin(x1(t)-x3(t)))*cos(x1(t)-x3(t)))/(r*(2*M-m*cos(x1(t)-x3(t))^2+m))}; ics := x1(0) = (1/6)*Pi, x2(0) = 0, x3(0) = (1/2)*Pi, x4(0) = 0; sol := dsolve(`union`(sys, {ics}), numeric, output = listprocedure, abserr = 0.1e-7, relerr = 0.1e-7); theta_plot := odeplot(sol, [t, x1(t)], 0 .. 10, color = blue, title = "θ(t)", thickness = 2); phi_plot := odeplot(sol, [t, x3(t)], 0 .. 10, color = red, title = "φ(t)", thickness = 2); display([theta_plot, phi_plot], title = "Hoekbeweging in de tijd", legend = ["θ(t)", "φ(t)"])

{diff(x1(t), t) = x2(t), diff(x2(t), t) = (-29.43*sin(x1(t))-.3*x4(t)^2*sin(x1(t)-x3(t))+(9.81*sin(x3(t))-.7*x2(t)^2*sin(x1(t)-x3(t)))*cos(x1(t)-x3(t)))/(3.5-.7*cos(x1(t)-x3(t))^2), diff(x3(t), t) = x4(t), diff(x4(t), t) = 3.333333333*(-49.05*sin(x3(t))+3.5*x2(t)^2*sin(x1(t)-x3(t))+(29.43*sin(x1(t))+.3*x4(t)^2*sin(x1(t)-x3(t)))*cos(x1(t)-x3(t)))/(5-cos(x1(t)-x3(t))^2)}

 

 

NULL


 

Download lagragne_popgave_wiel_omgezet_naar_1e_orde_systeem_ode_30-7-2025.mw

2 3 4 5 6 7 8 Last Page 4 of 77