Green's Function Solution in Matlab

I wrote a code for stationary temperature distribution which is;

f = 0.09;
b = 0.0044;
q = 3.73e-9;
L = 1;
Tw = 250;
Tam = 27;
syms c x g
T = 2*c*cosh(x*((f-b*g)/q)^0.5)+g/(f-b*q);
c = solve(subs(T,'x',L/2)==0,c);
z = simplify(int(subs(T,'c',c),x,-L/2,L/2));
g = solve(z==L*(Tw-Tam),g)

T

Now I try to use perturbation method and find T1 temperature distribution. My equation is;

q*T1''(x)-T1(x)*(f-b*g+i*w*p)=T*(f1-b*g1)-g1

And my prof. say that can be solved by Green's function G(x,y), where the G(x,y) is soluton of;

q*G''(x,y)-G(x,y)*(f-b*g+i*w*p)=Diract(x-y)

And

G(L/2)=0 , G(-L/2)=0

if I can find the G(x,y), I will get the solution of T1.

Can anyone help me please?

Thank you for everything

Yusuf

Answer 1

[Ted Shultz]

The person who asked this question has some (slow) sample code that will solve Green's function.

Answer 2

[Sousheel reddy Lachagari]

clc
clear all
syms x y 
N = input('the value of N = ');
M = input('M=')
My=diff(M,y)
Nx=diff(N,x)
x1=input('limit of x');
x2=input(' x u');
y1=input('limit of y');
y2=input('limit of y up');
I=int(int(Nx-My,y,y1,y2),x,x1,x2);

[Walter Roberson]

Please add comments explaining how this implements Green's function.