Hi everyone,

I am trying to visualize the integral of x³ over the interval x=−1 to x=1. I tried using:

with(plots):
display(
   plot(x^3, x = -1 .. 1, color = black, thickness = 2),
   shadebetween(x^3, 0, x = -1 .. 1, color = cyan)
);

This works, but I wondered if there’s a better or more elegant way to visualize definite integrals. For example:

• Can I add transparency to the shaded region?
• Is there a built-in function that directly plots definite integrals with shading?
• Any tips for improving the aesthetics of such plots?

Thanks in advance for any help!

Download εμβαδόν_χωρίου.mw

I am working on visualizing the results of an SIR model in Maple. I’ve created multiple plots for different parameter values using dsolve and odeplot, and I’d like to arrange them in a grid (e.g., 2 rows x 3 columns) for easier comparison. Here’s what I’ve tried:

Download grid.mw

Unfortunately, the insequence=true and gird=[2, 4] options didn’t work as I expected. Instead of arranging the plots in a grid, they appear like an animation

How can I correctly arrange multiple plots into a grid layout in Maple? Are there specific options or steps I’m missing? I would appreciate any guidance or examples.

I'm working on solving a system of ODEs in Maple that models an epidemic scenario, tracking the number of susceptible, infected, and recovered individuals over time. Here's what I've done so far:
Download sir_model.mw

Created a table of results at daily intervals with the following code:

results := [seq([t = tval, s = round(evalf(sol(tval)[2][1])), i = round(evalf(sol(tval)[2][2])),r = round(evalf(sol(tval)[2][2]))], tval = 0 .. 50)];

printf("%-10s %-15s %-15s %-15s\n", "Day", "Infected", "Recovered");
printf("---------------------------------------------\n");
for entry in results do
    printf("%-10d %-15d %-15d %-15d\n", entry[t],entry[s], entry[i], entry[r]);
end do;

but I encountered an error (in fprintf) integer expected for integer format

Hello,

I am working on a problem involving Newton's Law of Cooling in Maple, I want to:

1. Solve for the cooling constant k,
2. Find how long it takes for the pizza to cool to 300 F,

Here is the code I am currently using:

with(plots):

T__s := 75;  # Surrounding temperature
ode := diff(T(t), t) = k*(T(t) - T__s);

# General solution
general_solution := dsolve(ode, T(t));

# Initial condition: T(0) = 350
T := t -> 75 + exp(k*t)*c__1;
C := solve(350 = 75 + exp(0)*C__1, C__1);

# Solve for k using T(5) = 340
k := solve(340 = 75 + 275*exp(5*k), k);
k := evalf(k);

# Solve for time to reach T = 300
time_to_reach := solve(300 = 75 + 275*exp(k*t), t);
time_to_reach := evalf(time_to_reach);

# Define T(t) with solved constants
T := t -> 75 + 275*exp(k*t);

# Plot temperature over time
temp := plot(T, t = 0 .. 280, 75 .. 360, title = "Exponential Decay of Pizza Temperature", labels = ["t (minutes)", "T(t) (°F)"], color = blue);

# Plot T__s as a horizontal line
T__sur := plot(T__s, t = 0 .. 280, color = red, linestyle = dash);

# Combine and display plots
display([temp, T__sur]);

The code works, but I am wondering if there is a more efficient or genuine way.

I would greatly appreciate any advice or alternative approaches to solving and visualizing this in Maple.

Thank you in advance for your insights!

PDF ready for download.

Error occurred during PDF generation. Please refresh the page and try again

Is there another way to change the x-axis values other than tickmarks? For example on an interval [0,4*Pi]

Download directly_integrable.mw