Solving Equation Sinus: Best Software & Arguments

  • Thread starter Thread starter Deuskull
  • Start date Start date
Click For Summary

Discussion Overview

The discussion revolves around solving a specific trigonometric equation involving sine functions, where the goal is to find either an analytical solution or multiple numerical solutions. The context includes the application of this equation to diffusion in porous materials.

Discussion Character

  • Exploratory
  • Technical explanation
  • Mathematical reasoning

Main Points Raised

  • One participant presents the equation a.sin(x(1+k)) + b.sin(x(1-k)) - c.sin(x) = 0 and expresses difficulty in finding solutions manually or using software like Mathematica and Matlab.
  • Another participant suggests that there may not be an analytic solution and implies that the equation might already be in its simplest form.
  • A third participant shares a Mathematica implementation that involves plotting the function and using FindRoot to identify approximate roots, providing several numerical solutions.
  • A later reply acknowledges a misunderstanding of the original question and confirms the validity of the Mathematica solution using Mathcad.
  • The original poster reflects on the importance of plotting the curve and expresses intent to use the identified roots in subsequent calculations.

Areas of Agreement / Disagreement

There is no consensus on the existence of an analytical solution, with some participants suggesting it may not exist while others focus on numerical methods. The discussion remains unresolved regarding the best approach to solving the equation.

Contextual Notes

Participants note the potential complexity of the equation and the challenges in finding solutions, highlighting the dependence on software capabilities and the need for appropriate initial conditions for numerical methods.

Deuskull
Messages
2
Reaction score
0
Good morning everyone,

I'm working on an equation like :
a.sin(x(1+k))+b.sin(x(1-k))-c.sin(x) = 0
where x is the variable and the others one (a,b,c,k) are constant.

I tried to solve it manually but I didn't find any simple solutions.
I tried to solve it using some softwares (Mathematica and Matlab) using solve, nsolve and dsolve but without success.

What's the best software and/or the best argument (xsolve) to use ?
I want to find the analytical solution or, if it's not possible, I would like to have 10 positive numerical solutions.

Thank you.

PS : For anyone interested, this equations is a prior to use an equation of diffusion in porous material.

a = 5,95E-01
b = 5,22E-02
c = -3,53E-01
k = 0,14
 
Physics news on Phys.org
I'm not sure that there is an analytic solution ... I suspect that it's in it's simplest form already.

I may not have understood the problem, though. It seem as though any of the mathematical applications should handle giving you values as it looks as though you have, implicitly, a straightforward function of x. For example, implementing it in Mathcad gives the following ...

attachment.php?attachmentid=52273&stc=1&d=1351106971.jpg
 

Attachments

  • phys - 12 10 24 diffusion equation 01.jpg
    phys - 12 10 24 diffusion equation 01.jpg
    36.9 KB · Views: 787
In[1]:= a=5.95*10^-1;
b=5.22*10^-2;
c= -3.53*10^-1;
k=0.14;
f=a*Sin[x(1+k)]+b*Sin[x(1-k)]-c*Sin[x];
Plot[f,{x,0,25}]

Out[6]= <plot snipped>

By inspection pick starting points "near" each of the roots

In[7]:= Map[FindRoot[f==0,{x,#}]&, {0,3,6,9,12,14,17,20,22,24}]

Out[7]= {{x -> 0.}, {x -> 2.9177220113446776}, {x -> 5.822871346236126}, {x -> 8.701118634689632}, {x -> 11.534799430135319}, {x -> 14.302537173735368}, {x -> 16.98315627393366}, {x -> 19.56859265877026}, {x -> 22.08324005185527}, {x -> 24.587711402920295}}
 
Last edited:
Silly me ... didn't read the question properly. A Mathcad solution validating the Mathematica one:

attachment.php?attachmentid=52276&stc=1&d=1351118674.jpg
 

Attachments

  • phys - 12 10 24 diffusion equation roots 01.jpg
    phys - 12 10 24 diffusion equation roots 01.jpg
    26.1 KB · Views: 733
Thank you for your answers.

I was too much trying to find a solution that I didn't plot the curve.
I will try to take all the roots and implement them in my next equation.

Have a good day !
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 1 ·
Replies
1
Views
932
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 14 ·
Replies
14
Views
4K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 5 ·
Replies
5
Views
3K