Помогите реализивать в среде matlab программу

[katoraxa]

Проблема в реализации в среде matlab программы по расчету несобственного интеграла методом Гаусса. На Delphi 7 кое как получилось реализовать, а вот в Matlab всю голову уже сломал по реализации данного метода. Или я что то не понимаю или лыжи совсем не едут.
Вот вроде как работает на Delphi такая реализация но это с уже заданной функцией, а не с вводимой(
Function Fu(x:double):double;
Begin
Fu:=0.265*(4.774-((x*x)))*exp(-101.4*(1-(x))*(1-(x)))/(exp(h*c/k/T/x/L0)-1);
End;
Procedure TForm1.Button1Click(Sender: TObject);
Var
I:integer;
Begin
Series1.Clear;
H:=6.63*intpower(10,-34);
C:=3*intpower(10,8);
K:=1.38*intpower(10,-23);
T:=400;
Tk:=1000;
L0:=558*intpower(10,-9);
A:= 2*h*c*c/L0/L0/L0/L0/L0;
While T<Tk do begin
T:=T+1;
S:=0;
for i:=0 to 4 do begin
S:=S+Ck[i]*Fu(Xk[i]);
end;
S:=S*a;
Series1.ADDXY (T, S);
end;
end;
End.

[Taus]

Если вы хотите работать с вводимой функцией, то вам придётся писать парсер математических выражений или использовать готовые варианты (не могу сказать).