3Â¥: Originally posted by תʯǧØð at 2013-05-21 14:24:09
MÎļþ£º
function =ji(m)
a=quad('power(power(r,5)-power(m,5),-1./2)',m,973.15)
£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­
function =ji(m)
a=quad('power(power(r,5)-power(50,5),-1./2)' ...

function a = myji(m)
myfun = @(r,m)power(power(r,5)-power(m,5),-1./2);
a = quadl(@(r)myfun(r,m),m,973.15);
end

2Â¥: Originally posted by change0618 at 2013-05-20 14:29:33
±ßÖµÎÊÌ⣬¶øÇÒÊǶþ½×µÄ£¬Ïȱ仯³ÉÒ»½×·½³Ì×飬ȻºóÓÃbvp4c½â¾ö

3Â¥: Originally posted by תʯǧØð at 2013-05-21 14:24:09
MÎļþ£º
function =ji(m)
a=quad('power(power(r,5)-power(m,5),-1./2)',m,973.15)
£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­£­
function =ji(m)
a=quad('power(power(r,5)-power(50,5),-1./2)' ...

function MyInteger
clear all
clc
m = 10:10:50;
for i = 1:length(m)
    a(i) = quadl(@(t)myfun(t,m(i)),m(i),973.15);
end
a
end
function f = myfun(t,m)
f = power(power(t,5)-power(m,5),-1./2);
end