m_magooda

Ask a Question

Create a Post

One Badge

Member for: 7 years, 92 days

Contact m_magooda

MaplePrimes Activity

These are replies submitted by m_magooda

@mmcdara yes, it's a heat equ...

Answered: m_magooda 20

April 20 2018

0 0

@mmcdara

yes, it's a heat equation using finite difference method

I wanted to show the role of many parameters, but it's not pedagogical so i guess fsolve would do just fine

now I realize how mistaken i was regarding the order of equations, it would take like forever and no added value would be achieved

thank you again for your assistance, much appreciated

@mmcdara Thank you for the assist...

Answered: m_magooda 20

April 19 2018

0 0

@mmcdara

Thank you for the assistance, that was helpful

the 53 file worked just fine, i had problem with defining eqs as you mentioned

regarding the 54 file, I tried plugging values to anything that isn't T[1..9] so now I have a system of 9 equations with 9 T values

yet maple freezes? is it that hard with only couple of power 4 terms?

Fsolve works good, but I'm having a problem that may sound stupid ... how do i know which T is which in the array, does it use the numbering i used initially in defining the equations?

thank you

$eqs := {T[4]+T[4]+T[8]+T[8]-4*T[5] = 0, (1/2)*hi*`Δx`*(Ti-T[1])+(1/2)*k*`Δy`*(T[2]-T[1])/`Δx`+(1/2)*k*`Δx`*(T[3]-T[1])/`Δy` = 0, (1/2)*k*`Δy`*(T[1]-T[2])/`Δx`+(1/2)*hi*`Δx`*(Ti-T[2])+k*`Δx`*(T[4]-T[2])/`Δy` = 0, (1/2)*k*`Δy`*(T[8]-T[9])/`Δx`+(1/2)*ho*`Δx`*(To-T[9])+(1/2)*epsilon*sigma*`Δx`*(Ts^4-T[9]^4) = 0, T[4]+T[1]+T[4]+T[6]-4*T[3] = 0, (1/2)*k*`Δx`*(T[3]-T[6])/`Δy`+(1/2)*k*`Δy`*(T[7]-T[6])/`Δx`+(1/2)*ho*`Δx`*(To-T[6])+(1/2)*epsilon*sigma*`Δx`*(Ts^4-T[6]^4) = 0, T[3]+T[2]+T[5]+T[7]-4*T[4] = 0, (1/2)*k*`Δy`*(T[6]-T[7])/`Δx`+k*`Δx`*(T[4]-T[7])/`Δy`+(1/2)*k*`Δy`*(T[8]-T[7])/`Δx`+ho*`Δx`*(To-T[7])+epsilon*sigma*`Δx`*(Ts^4-T[7]^4) = 0, (1/2)*k*`Δy`*(T[7]-T[8])/`Δx`+k*`Δx`*(T[5]-T[8])/`Δy`+(1/2)*k*`Δy`*(T[9]-T[8])/`Δx`+ho*`Δx`*(To-T[8])+epsilon*sigma*`Δx`*(Ts^4-T[8]^4) = 0}$

${2*T[4]+2*T[8]-4*T[5] = 0, 2006.025000-4.900000000*T[1]+.7000000000*T[2]+.7000000000*T[3] = 0, .7000000000*T[1]-5.600000000*T[2]+2006.025000+1.400000000*T[4] = 0, 2*T[4]+T[1]+T[6]-4*T[3] = 0, .7000000000*T[8]-15.35750000*T[9]+304.8097426-0.2551500000e-8*T[9]^4 = 0, T[3]+T[2]+T[5]+T[7]-4*T[4] = 0, .7000000000*T[3]-16.05750000*T[6]+.7000000000*T[7]+304.8097426-0.2551500000e-8*T[6]^4 = 0, .7000000000*T[6]-32.11500000*T[7]+1.400000000*T[4]+.7000000000*T[8]+609.6194853-0.5103000000e-8*T[7]^4 = 0, .7000000000*T[7]-32.11500000*T[8]+1.400000000*T[5]+.7000000000*T[9]+609.6194853-0.5103000000e-8*T[8]^4 = 0}$