## 55 Reputation

6 years, 260 days

## Get result from pdsolve...

Maple

hello ,

how i can exract value from pdsolve ,i need to use dU(x,R)/dR

thank you

 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 >
 (6)
 >
 >
 >
 >

## Fourth order differential equation...

Maple

Hello evrey one , I need help for solve these equation with boundary conditions

Boundary Conditions

My COde + equation

 >
 >
 >
 >
 (1)
 >
 >
 >
 >
 >

Thank you

## plot simpson integral...

hello i need plot integrale of siampson thank you

> Simpson := proc(f, a, b, n)
> local h, S1, S2, S, i;
> h := (b-a)/n;
> S1 := 0.0;
> for i from 0 to n-1 do
> S1 := S1 + f(a + (2*i+1)*h);
> end do;
> S2 := 0.0;
> for i from 1 to n-1 do
> S2 := S2 + f(a + (2*i)*h);
> end do;
> S := (h/3) * ( f(a)+f(b) + 4*S1 + 2*S2 );
> return S;
> end proc:
> Digits := 5;

Digits := 5

> f := x -> 1/sqrt(39.24*x-44.65*(x*arccos(x)-sqrt(1-x^2)-13.88*(1-x^2)^1.5));

f := x -> 1/sqrt(39.24 x

2                2 1.5
- 44.65 (x arccos(x) - sqrt(1 - x ) - 13.88 (1 - x )   ))

> Simpson(f, 0, 1, 100):
>
> p:=int(f(x), x=0..0.1);

p := 0.0038931

> w:=int(f(x), x=0.1..0.2);

w := 0.0039570

> m:=int(f(x), x=0.2..0.3);

m := 0.0040826

> l:=int(f(x), x=0.3..0.4);

l := 0.0042836

> kohv:=int(f(x), x=0.4..0.5);

kohv := 0.0045860

> q:=int(f(x), x=0.5..0.6);

q := 0.0050373

> s:=int(f(x), x=0.6..0.7);

s := 0.0057306

> d:=int(f(x), x=0.8..0.9);

d := 0.0089874

> f:=int(f(x), x=0.9..1);

f := 0.013349

## fsolve b*x*ln(x)-x*ln(a)+a=0...

hello, i went solve these equation ,with a & b take any value

b*x*ln(x)-x*ln(a)+a=0

thank you

## Warning, cannot evaluate the solution fu...

Maple

Pleaz i nees help i have probleme withe singularity

 >

Paramétres

 >
 (1.1)

Equation suivant x :

 >
 (2.1)

Equation suivant z :

 >
 (3.1)

Equation suivant y :

 >
 (4.1)

Equation suivant y

 >
 (5.1)

Résolution :

 > CI:= x(0)=0,z(0)=0,theta(0)=0,alpha(0)=0,D(x)(0)=0,D(alpha)(0)=0,D(z)(0)=0,D(theta)(0)=0;
 (6.1)
 > if theta(t) <> 0 then  solution:=dsolve([eq1,eq2,eq3,eq4,CI],numeric,maxfun=0):  odeplot(solution, [[t, x(t)]], t = 0 .. 100, thickness = 2);  odeplot(solution, [[t, z(t)]], t = 0 .. 100, thickness = 2);  odeplot(solution, [[t, theta(t)]], t = 0 .. 100, thickness = 2);  odeplot(solution, [[t, alpha(t)]], t = 0 .. 100, thickness = 2);  #odeplot(solution,[[t,x(t)],[t,alpha(t)],[t,z(t)],[t,theta(t)]], t=0..100, thickness=2);  end ;

thank you !