farahnaz

Ask a Question

Create a Post

Name: Dr. farahnaz farahnaz

2 Badges

Member for: 9 years, 321 days

Contact Dr. farahnaz farahnaz

MaplePrimes Activity

These are questions asked by farahnaz

HOW I CAN DSOLVE DIFFERENTIAL EQUATION W...

Asked: farahnaz 15 Product: Maple 2015

July 02 2016

0 5

HI.please help me for solve differenrtial equation with finite difference method not dsolve numeric solver in maple

thanks

FDM2.mw

L := 1/50000000; -1; eta := 1; -1; PDE[111] := 7.65692307692309*10^(-8)*(diff(f1(x), x, x, x, x))-1.56784615384616*10^12*(diff(f1(x), x, x))+220.592307692308*(diff(f2(x), x, x, x))-3.52947692307693*10^21*(diff(f2(x), x))+43.7538461538462*(diff(f3(x), x, x, x))+4.81292307692309*10^20*(diff(f3(x), x))+6.50473846153848*10^30*f1(x)-7.90000000000000*10^(-8)*eta*f1(x)

0.7656923077e-7*(diff(diff(diff(diff(f1(x), x), x), x), x))-0.1567846154e13*(diff(diff(f1(x), x), x))+220.592307692308*(diff(diff(diff(f2(x), x), x), x))-0.3529476923e22*(diff(f2(x), x))+43.7538461538462*(diff(diff(diff(f3(x), x), x), x))+0.4812923077e21*(diff(f3(x), x))+0.6504738462e31*f1(x)

(1)

PDE[222] := 2.14211538461539*10^(-8)*(diff(f2(x), x, x, x, x))-1.64988461538462*10^12*(diff(f2(x), x, x))+7.90486153846156*10^30*f2(x)-220.592307692308*(diff(f1(x), x, x, x))+3.52947692307693*10^21*(diff(f1(x), x))-5.94323076923080*10^11*(diff(f3(x), x, x))+5.13378461538463*10^30*f3(x)-7.90000000000000*10^(-8)*eta*f2(x)

0.2142115385e-7*(diff(diff(diff(diff(f2(x), x), x), x), x))-0.1649884615e13*(diff(diff(f2(x), x), x))+0.7904861538e31*f2(x)-220.592307692308*(diff(diff(diff(f1(x), x), x), x))+0.3529476923e22*(diff(f1(x), x))-0.5943230769e12*(diff(diff(f3(x), x), x))+0.5133784615e31*f3(x)

(2)

-0.6380769231e-30*(diff(diff(diff(diff(diff(diff(f3(x), x), x), x), x), x), x))+0.9665370462e-7*(diff(diff(diff(diff(f3(x), x), x), x), x))-0.3101547535e13*(diff(diff(f3(x), x), x))-43.7538461538462*(diff(diff(diff(f1(x), x), x), x))-0.4812923077e21*(diff(f1(x), x))-0.5943230769e12*(diff(diff(f2(x), x), x))+0.5133784615e31*f2(x)+0.2299890587e32*f3(x)-0.7901053333e-7*omega^2*f3(x)

(3)

bcs := {f1(0) = 0, f1(L) = 0, f2(0) = 0, f2(L) = 0, f3(0) = 0, f3(L) = 0, ((D@@1)(f1))(0) = 0, ((D@@1)(f1))(L) = 0, ((D@@1)(f2))(0) = 0, ((D@@1)(f2))(L) = 0, ((D@@1)(f3))(0) = 0, ((D@@1)(f3))(L) = 0, ((D@@2)(f3))(0) = 0, ((D@@2)(f3))(L) = 0}

{f1(0) = 0, f1(1/50000000) = 0, f2(0) = 0, f2(1/50000000) = 0, f3(0) = 0, f3(1/50000000) = 0, (D(f1))(0) = 0, (D(f1))(1/50000000) = 0, (D(f2))(0) = 0, (D(f2))(1/50000000) = 0, (D(f3))(0) = 0, (D(f3))(1/50000000) = 0, ((D@@2)(f3))(0) = 0, ((D@@2)(f3))(1/50000000) = 0}