24小时热门版块排行榜    

查看: 827  |  回复: 8
当前主题已经存档。

yongbinz

金虫 (著名写手)

[交流] 【讨论】一个小程序,过来讨论讨论!

我需要一个小程序,求解一个方程!
X/L=Sin(θ/2)÷( θ/2) 公式中的θ为弧度,已知X与L,输入X与L后,要能求出θ的大小。
这个程序怎么写呢,大家能否给点意见?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gjliu

银虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
幻影无痕(金币+0,VIP+0):呵呵,亲自上阵!鼓励一下! 6-3 11:04
这个是一个超越方程求解的问题,可以看一些《计算方法》或者《数值分析》的书,比如牛顿迭代什么的方法都可以找到数值解的。如果是Matlab或者Mathematic等软件直接有求解方程的函数,相对比较简单。

就是不知道你要用什么语言去写。给你一个简单的例子吧


用牛顿迭代法求方程X3+4X2-10=0 的根

#include
main()
{float x,x0,f,f1;
x=1.5;
do
{x0=x;
f=x0*x0*(x0+4)-10;
f1=x0*(3*x0+8); /* f的导数,即它的切线 */
x=x0-f/f1; /* 切线与x轴的交点 */
}
while(fabs(x-x0)>=1e-5);
printf("The root is %8.7f\n",x);
}

[ Last edited by gjliu on 2009-6-3 at 11:21 ]
2楼2009-06-03 11:00:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongbinz

金虫 (著名写手)

谢谢楼上,不过,本人就是不会C语言呀!谢谢支持!
3楼2009-06-03 11:40:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

levone

木虫 (正式写手)

那会什么语言?
莫愁前途无知己!
4楼2009-06-03 12:31:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuzhs

金虫 (正式写手)

maple,mathmatics
当所有一切都已看平淡 是否有一种精致还留在心田
5楼2009-06-06 00:14:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duzhuqi

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
suton(金币+1,VIP+0):谢谢! 6-7 09:10
MATLAB程序如下:
function y=f6(x)
y=Sin(x)/( x)-X/L;
end
----------------------------------------------------
a=-1;b=1.5;
y1=f6(a);y2=f6(b);
while abs(b-a)>0.001
x0=0.5*(a+b);y0=f6(x0);
  if abs(y0)< 0.0001
    x0
    elseif y0*y1<0
        b=x0;y2=y0;
    else
        a=x0;y1=y0;
  end
end
x0
6楼2009-06-06 22:28:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongbinz

金虫 (著名写手)

真的很感谢呀!谢谢!
7楼2009-06-06 23:42:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

levone

木虫 (正式写手)

f6(x)是什么意思?
莫愁前途无知己!
8楼2009-06-07 08:40:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sitonyl

铁杆木虫 (文坛精英)


小木虫(金币+0.5):给个红包,谢谢回帖交流
程序中自己定义的函数。
zz
9楼2009-06-09 09:12:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yongbinz 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-29 17:22 by 4FFAWE8HcgUD
[基金申请] 面上意见出来了 +6 黄鸟于飞Chao 2026-08-29 12/600 2026-08-29 17:20 by 黄鸟于飞Chao
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 16:31 by 时光旅行者
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:51 by jCd0dEvKHShX
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 11:49 by jCd0dEvKHShX
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 11:18 by jCd0dEvKHShX
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
信息提示
请填处理意见