急；这个四阶龙伯格求解方程组那里出错了？
1）设：粗加工加工x天，细加工加工y天；方程组（1）5x+15y=150（2）x+y=14 解:x=6,y=8 2）设甲组为x人，乙组为y人；方程组（1）x+y=44;(2)x+6=y-6；解之得：x=16;y=28
3)设甲粮仓为x吨，乙粮仓y豚；方程组（1）14+x=y-10；（2）6（x+8）=y+18；解x=18,y=42
设：快车的速度为x米/秒，慢车的速度为y米/秒
4x+4y=168+184
16x-16y=168+184
把上面的方程组整理可得：x+y=88
x-y=22
解方程组得：x=55 y=33
a是积分下限，b是积分上限，eps是精度，f是被积函数。
double romb(double a,double b,double eps,double (*f)(double))
{ int m,n,i,k;
double y[10],h,ep,p,x,s,q;
h=b-a;
y[0]=h*((*f)(a)+(*f)(b))/2.0;
m=1; n=1; ep=eps+1.0;
while ((ep>=eps)&&(m<=9))
{ p=0.0;
for (i=0;i<=n-1;i++)
{ x=a+(i+0.5)*h;
p=p+(*f)(x);
}
p=(y[0]+h*p)/2.0;
s=1.0;
for (k=1;k<=m;k++)
{ s=4.0*s;
q=(s*p-y[k-1])/(s-1.0);
y[k-1]=p; p=q;
}
ep=fabs(q-y[m-1]);
m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
}
return(q);
}
调用实例(C语言实现)：
//---------------------------------------------------------------------------
#include <stdio.h>
#include <math.h>
double rombf(double x)
{ double y;
y=x/(4.0+x*x);
return(y);
}
double romb(double a,double b,double eps,double (*f)(double))
{ int m,n,i,k;
double y[10],h,ep,p,x,s,q;
h=b-a;
y[0]=h*((*f)(a)+(*f)(b))/2.0;
m=1; n=1; ep=eps+1.0;
while ((ep>=eps)&&(m<=9))
{ p=0.0;
for (i=0;i<=n-1;i++)
{ x=a+(i+0.5)*h;
p=p+(*f)(x);
}
p=(y[0]+h*p)/2.0;
s=1.0;
for (k=1;k<=m;k++)
{ s=4.0*s;
q=(s*p-y[k-1])/(s-1.0);
y[k-1]=p; p=q;
}
ep=fabs(q-y[m-1]);
m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
}
return(q);
}
int main(void)
{
double a,b,eps,t,rombf(double);
a=0.0; b=1.0; eps=0.000001;
t=romb(a,b,eps,rombf);
printf("\n");
printf("t=%e\n",t);
printf("\n");
return 0;
}
//---------------------------------------------------------------------------
这种问题很少能得到正确而及时的解答……无能为力，唉
想看个丛林方面的电影求大家推荐?
无论是《侏罗纪公园》2中的丛林猎杀场景还是《侏罗纪世界》1中迅猛龙小队出击的画面,都给我留下了...第四部分:人与灵(也就是面对自己)。高级阶段 这也是影片后半段大篇幅叙述尤西自与凯文分离之后,面对...金伯格是 探索 频道纪录片《生存大挑战》中的冒险者,他曾与两个朋友、一位向导穿越亚马逊热带丛林,...
求介绍好看的电影
求介绍好看的电影 我来答 4个回答 #热议# 没有文化的年迈农民工退休后干点啥好? batesmotel...终于,本杰明迎来了自己在人世上生存的第11个年头,此时的他看起来也似乎健康年轻一些,但依旧老态龙钟...影片根据2009年美利坚航空公司1549号航班机长切斯利·舒伦伯格的真实英雄事迹改编。讲述了萨利机长在...
请问关于宝藏的电视剧或电影有哪些!!!
吉祥(洪金宝饰〕、刘曦〔吴京饰〕护送“贡佛天珠”来港途中被四名武功高强的人物抢走!最令人感到惊讶的...电影《木乃伊3:龙帝之墓》 主演:布兰登·弗莱瑟 杨紫琼 李连杰 黄圣依 剧情简介: 李连杰扮演一名公元前...2013-10-25 求关于找宝藏的电影!!! 36 2010-10-22 寻宝藏的电影有哪些 5 2010-01-02 请推荐...
如何做一个总经理
总经理绝不能事必躬亲、事事想管，著名企业家路易士·龙伯格认为，总经理履行三大职责：一、为公司的未来设定战略目标和发展方向;二、确定合适的人才是否被安置到合适的位置，考虑目前的同时应考虑未来之需求和变化;三、查证公司各阶层之每一个人，对于预先设定之目标与期望是否确已达成。三、总...
周杰伦成名历程?
周杰伦出生于台湾省新北市，祖籍福建省永春县。周杰伦4岁的周杰伦小时候周杰伦小时候(17张)时候，母亲叶惠美把他送到淡江山叶幼儿音乐班学习钢琴。周杰伦读初中二年级时，他的父母因性格不合离婚，而周杰伦则归母亲叶惠美抚养。中考时，周杰伦没有考上普通高中。同年，他因为擅长钢琴而被淡江中学第一届音乐...
有哪些罕见的姓氏?
罕见的姓氏，这些姓氏可能只在某些特定的地区或族群中出现：邢薇（Xing Wei）：这是一个中国姓氏，据说只有几十人拥有这个姓氏。大谷（Ooya）：这是一个日本姓氏，在日本非常罕见。帕特里奇（Patrick）：这是一个英国姓氏，在英国不算常见。乌姆比加（Umbijak）：这是一个印尼的巴塔克族姓氏。布龙伯格（...
什么是歌舞秀
赫伯特及后来的龙伯格(Sigmund Romberg)、弗雷蒙(Roudolf Friml)等人的小歌剧在艺术特征上依然与维也纳小歌剧一脉相承,较为注重音乐的表现,强调音乐的主导地位,追求曲调的优美动听,然而对戏剧的展开和各种艺术因素的有机结合,则没有考虑为特定的人物性格和特定时间地点写作,只把精彩的布置,好看的戏服,优美的曲调,演员...
三角龙的资料
三角龙（属名：Triceratops）鸟臀目角龙下目角龙科的草食性恐龙的一属，化石发现于北美洲的晚白垩纪晚马斯特里赫特阶地层，约6800万年前到6500万年前。三角龙是最晚出现的恐龙之一，经常被作为晚白垩纪的代表化石。目前已有两个有效种：褶皱三角龙、究竟三角龙，但还有其他属被命名。三角龙也是最著名的恐龙...
恐龙的头上有三个角是什么龙
三角龙。三角龙是一种中等大小的四足恐龙，全长6—10米、高2.4-2.8米、重5—10吨。他们有非常大的头盾，以及三根角状物，令人联想起现代犀牛。传统上这些结构被认为是用来抵抗掠食者的武器，但新的理论认为这些结构可能用在求偶，以及展示支配地位，如同现代驯鹿、山羊、独角仙的角状物。三角龙最显著...
iphone4拳皇2012 我想知道连招的 近P>下S >GAUGE 是什么意思 还有带星...
空中P怎么打出来小白求高手详细解释下... 空中P 怎么打出来 小白 求高手详细解释下 展开 ...特瑞·伯格特殊技步步高:↘+P回冲拳:→+P必杀技★能量波:↓↘→+P★火焰冲拳:↓↙←+P 0r P...龙爪击:空中↓↙←+P★超球弹:↓↙←+P龙扑鼓:→↓↘+P(连打)★龙倒打:↓↘→+P┗★剥鳞...
向燕卡前： 主程序 clc;clear x0=0;y0=1; h=0.1; Nstep=floor(1/h); [x2,y2]=Rung_Kutta(@myfun,x0,y0,h,Nstep); plot(x2,y2) xlabel('x') ylabel('y') box off grid on 函数程序1 function y1=myfun(x,y) %% a self-defined function y1=sqrt(x+y); end 函数程序2 function [x,y]=...
枣阳市18421048825： 急!!!求matlab 用四阶龙格 - 库塔法求解常微分方程 - ？
向燕卡前： 建立.m文件 --------------------------------------------- function theta=danbai(t,X) x=X(1); dx=X(2); ddx=-sin(x); theta=[dx;ddx]; ---------------------------------------------- 命令窗口输入 >> [t,Y]=ode45(@danbai,[0 6],[pi/3 -1/2]); >> plot(t,Y(:,1),'ro-',t,Y(:,2),'bv-'); >> legend...
枣阳市18421048825： 急求MATLAB编程源代码用四阶龙格库塔法解如下微分方程 y'=y - 2x/y(0<x<1),y(0)=1,步长为h=0.2 - ？
向燕卡前： % 以下另存为文件 myrk4.m function [x,y]=myrk4(ufunc,y0,h,a,b)%参数: 函数名称,初始值向量,步长,时间起点,时间终点 n=floor((b-a)/h);%求步数 x(1)=a;%时间起点 y(:,1)=y0;%赋初值%按龙格库塔方法进行求解 for ii=1:nx(ii+1)=x(ii)+h;...
枣阳市18421048825： 在matlab中用四阶龙格 - 库塔法解二阶微分方程怎么做??最好有代码!! - ？
向燕卡前： 例 y'= - y+x+1,y(0) = 1 首先建立M-文件 (weif.m) function f = weif(x,y) f=-y+x+1; 求解:[x,y]=ode45('weif',[0,1],1) 再如:建立文件:function dy = rigid(t,y) dy = zeros(3,1); % a column vector dy(1) = y(2) * y(3); dy(2) = -y(1) * y(3); dy(3) = -0.51 * y(1)...
枣阳市18421048825： 急急急!求助matlab用龙格 - 库塔方法求解方程组 - ？
向燕卡前： function df=ode45_fun(t,xyzuvw) %%注意小写的v和大写的V %常数(请修正) R_0=1; rho_0=1; beta=1; G=6.67e-11; M=1.5e24;x=xyzuvw(1); y=xyzuvw(2); z=xyzuvw(3); u=xyzuvw(4); v=xyzuvw(5); w=xyzuvw(6);R=sqrt(x*x+y*y+z*z); V=sqrt(u*u+v*v+w...
枣阳市18421048825： matlab用四阶龙格库塔法解微分方程组; - ？
向燕卡前： function dy=test(t,y) dy=[-12*cos(y(2))-120*cos(208*2*pi/360-3*y(2));(12*sin(y(2))+120*sin(208*2*pi/360-3*y(2)))/y(1);];[t,y]=ode45('test',[0.01,1],[1,1])plot(t,y(:,1),t,y(:,2)); x=3000*sin(70*2*pi/360)-y(1).*sin(y(2)); z=3000*cos(70*2*pi/360)-...
枣阳市18421048825： Matlab用四阶龙格库塔法求解不可微分方程组的初值？
向燕卡前： %% function dx=aaa(t,x)%定义matlab函数来描述方程 dx=[120-2*x(1)+2*x(2);2*x(1)-5*x(2)] %% x0=[0;0] [t,y]=ode45('aaa',[0,1],x0); plot(t,y)
枣阳市18421048825： 用四阶龙哥库塔法解三阶微分方程y'''+yy''=0 - ？
向燕卡前： 查MATLAB中的ode45的帮助文件即可完成 help ode45
枣阳市18421048825： MATLAB 用四阶龙格库塔解微分方程组 - ？
向燕卡前： 可以用ode45求解.需定义mass函数,用以计算质量矩阵,即化为下述形式求解: M(t,y) * dy/dt = f(t,y).请给出ABCDEF的具体值.
枣阳市18421048825： 用四阶龙格库塔法求解矩阵微分方程？
向燕卡前： global R M U syms t R=[ 0.0247,0,0,0,0,0; 0,0.0247,0,0,0,0; 0,0,0.0247,0,0,0; 0,0,0,0.0193,-0.0193,0; 0,0,0,0,0.0193,-0.0193; 0,0,0,1,1,1 ]; M=[ 0,0,0,-15727/10000*sin(5/12*pi+80*pi*t)*pi,15727/10000*cos(1/4*pi+80*pi*t)*pi,15727/10000*sin(1/12*...
