## 2220 Reputation

4 years, 275 days

## I assume that you are familar with Maple...

I assume that you are familar with Maple styles and in your question you want, for example, to change only the font of all Maple input in an existing document.

What you have to do is:

• hide all other screen content using view-> show/hide contents
• select all remaining content with the mouse
• change the font
• show the hinden content using view-> show/hide contents

## Units in (4) are correct if (4) is evalu...

Units in (4) are correct if t in (4) is evaluated with units. With units standard int must be mapped over a column vector (I don't know why that is).

 >
 >
 >
 >
 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 > restart; with(Units:-Standard):
 > `#mover(mi("a"),mo("→"))` := <0*Unit(('m')/'s'^2), g, 0*Unit(('m')/'s'^2)>
 (6)
 > map(int,(6),t)
 (7)
 > int~((6),t)
 (8)

## display...

I assume that you are able to generate plots for the exact and the approximate solution using the plot command.

Each of these plots have to be assigned to a name with the assignment operator :=

Then you can plot boths solutions in one plot using plots:- display command.

If you have these solutions in a worksheet and you can't manage to compare them, you could upload them using the green arrow.

## "Complex" integrand...

To your question: How to calculate this integral?

My preliminary answer is: Don't use Maples sophisticated algorithms.

Zoomed in it looks like this

This could explain why Maples sophisticated integration methods have difficulties. They will try to follow all the oscillations.

Annother reason could be: The function you want to integrate is not always real valued also (see the attachment).

PS:: I have just seen Roubens answer. Zoomed it looks also jagged.

## Gaps...

Do you want the new constant to have the highest index? If filling of gaps is acceptable

```n:=1;
do n++ until not(has(_C||n,myconstants)):
new_constant := _C ||(n)
```

Otherwise, in your example myconstants is a set and the element with the highest index is to the right.

```n:=1;
do n++ until is(_C||n=myconstants[-1]):
new_constant := _C ||(n)```

I guess there are shorter ways

## Technically possible...

but you will get three bulky solutions

For clarity I assume in the following that you combine the coefficients into new ones.
Your equations will look like this

 > eq1:=cos(3*phi)-a*cos(phi)=b
 (1)
 > eq2:=sin(3*phi)-a*sin(phi)=c
 (2)

Solve for phi

 > expand(eq1,trig); sols:=solve(%,[phi])
 (3)

You get 3 solutions for for phi that you can substitute one after each other into eq2

 > simplify(subs(sols[1],eq2))
 (4)
 >

Depending on the initial parameters you might be able to simplify further and select an equation that describes your problem

`evalf(Int(abs(sin(x^4))/(sqrt(x) + x^2), x = 0 .. 5))`